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Foreword 


On  January  13-14.  1977  an  Instrunantatlon  Meeting  was  held  In  the 
Silver  Spring.  Maryland  offices  of  the  U.S.  Nuclear  Regulatory 
Coonlsslon. 

In  the  meeting,  the  R/0  progress  of  two-phase  flow  Instrumentation  made 
by  various  NRC  contractors  was  reported.  An  Informal  survey  of  the 
state-of-art  was  made  by  experts  In  the  field.  Since  nary  request  were 
made  for  copies  of  the  proceedings  of  the  meeting.  It  was  decided  that 
the  Meeting  Minutes  and  the  contractors'  reports  shoul  je  published  In 
the  form  of  a  NUREG  report  for  public  dissemination. 

Ylh-yun  Hsu.  Chairman 
Instrumentation  Review  Group 
Systems  Engineering  Branch 
Division  of  Reactor  Safety  Research 
U.S.  Nuclear  Regulatory  Comnlsslon 
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MEETING  MINUTES 

TWO-PHASE  FLOW  INSTRUMENTATION  REVIEW  GROUP 
By 

Y.  Y.  Hsu.  USNRC 

Date;  January  13-14,  1977 

Place:  Room  150,  Wlllste  Building,  USNRC 

Attendees*  See  Attachmnt  1 

Agenda;  See  Attachment  2 


A  full  text  of  proceedings  will  be  published  as  an  NRC  report  In 
the  near  future.  A  sumnary  is  given  below. 

I.  Progress  Reports  by  NRC  Contractors  (January  13,  AH) 

1.  ORNL,  J.  Sheppard  -  subjects  reported  are: 

a.  Effect  of  drag  body  on  flow  pattern. 

b.  Improvements  of  drag  disk  signal  by  the  disperser 
screen. 

c.  Test  results  using  drag  screen. 

d.  Testing  of  ionization  chamber  as  a  detector. 

e.  Modeling  of  turbine  meter  dynamic. 

f.  Adaptive  cross-correlation  algorithm. 

g.  Interpretation  of  flow  measurements  from  spool 
piece. 

2.  RPI.  R.  Lahey  -  subjects  reported  are: 

a.  Steam-air  content  measurements  with  I.R.  technique 
-  analysis  and  test  facilities. 

b.  Novel  techniques 

(1)  test  facilities  and  electronics, 

(2)  x-ray  standard, 

(3)  exploratory  tests  with  I.R.,  R.F.  and  impedence 
probes. 
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3.  BNLt  0.  C.  Jones 

«.  Description  of  test  facility  and  probe. 

b.  Present  status  of  development  of  local  optical  probes. 

4.  SUNY  Stoney  Brook.  R.  Lee  •  reported  development  of  Laser- 
Ooppler  anemonety.  Including: 

a.  Operating  principles  and  difficulties. 

b.  Proposed  scheme  to  determine  number  density  of 
particles  as  function  of  velocity  and  particle  size. 

c.  Description  of  electronic  package. 

d.  Calibration  method. 

e.  Experimental  data  to  demonstrate  the  feasibility 
of  the  method. 

5.  INEL.  EG&G  Staff 

a.  DTT  and  densitometer  In  LOFT  -  R.  Wesley 

(1)  description  of  system. 

(2)  uncertainty  analysis. 

(3)  results  from  LOFT  tests. 

(4)  comparison  of  density  measurements  from  on^ 
and  from  ydensltometer. 

b.  Turbine  model  analysis  •  L.  Goodrich 

(1)  description  of  various  models, 

(2)  voltage  vs  free  field  velocity  curves  from 
WCL. 

(3)  comparison  of  velocities  of  liquid  and  vapor 
based  upon  various  model  vs  the  true  mass  flux. 
Proposed  "volumetric  theory"  model  Is  shown  to 
be  more  accurate  than  the  or  Rouhani  models. 

c.  Upgraded  drag  disc  turbine  separate  effects  tests  - 

R.  Bearden 

(1)  test  results  of  bearing  life, 

(2)  description  of  single  bearing  and  compound 
bearing  turbine, 

(3)  sunmary  of  single-phase  tests  of  prototype  OTT, 

(4)  scope  of  model  development  at  RPI. 

(5)  development  eddy  current  sensor. 

d.  Rake  design  and  pitot-tube  data  -  R.  Bearden  and 

S.  Banerjee 

(1)  description  of  location  and  design  requirement 
of  3-OTT  rake  for  LOFT  application. 
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(2)  pitot  tube  rake  developed  by  McMaster  Ur.lverslty. 

(3)  sample  results  of  using  pitot  tube  rake  to  measure 
tmo-phase  flow. 

e.  Full  flow  drag  screen  •  J.  Colson.  J.  R.  Fincke  (pre¬ 
sented  by  A.  Arave) 

(1)  principle  of  drag  on  a  blunt  body  discussed. 

(2)  present  drag  screen  design  (3-1ever  arm). 

(3)  air-water  test  results  on  round-wire  screen  (over 
a  wide  range  of  void  fraction,  capable  to  detect 
flow  pattern). 

f.  Transit  time  flow  meter  -  G.  Lassahn 

(1)  discussion  of  error  associated  with  transit  time 
flow  meter. 

(2)  presentation  of  calculated  delay  time  as  function 
of  averaging  time. 

(3)  presentation  of  sample  signal  correlation  using 
thermocouples. 

g.  Nuclear  hardened  densitometer  -  A.  Stephens 

(1)  present  setup  for  testing  In  Adv.  Testing  Reactor. 

(2)  methods  of  removing  background  radiation 

(a)  shielding, 

(b)  background  subtraction  (In  both  digital 
and  analog  modes). 

(3)  chordal  overage  density  statistical  error  position 
and  energy  dependent. 

II.  State-of-Art  Review  (January  13,  PM) 

1.  Local  probes  -  J.  M.  Oelhaye.  Grenoble 

a.  Theoretical  prerequisites  for  local  probe  measurements 
were  discussed.  Including 

(1)  speed  of  displacement  of  moving  Interface, 

(2)  phase  density  function, 

(3)  Instantaneous  space-averaging, 

(4)  local,  time-averaging, 

(5)  conmutatlvlty  of  averaging  operators. 

b.  Survey  of  recent  work  on  local  probes 

(1)  electrical  probes  -  recent  work  at  Karlsruhe 
was  discussed, 

(2)  Laser-Ooppler  volocimetry  -  papers  presented  at 
LOA-symposlum  (including  those  by  Durst  and  Zare 
and  by  Oelhaye),  at  Osaka  University  and  In  AIChE 
Journal  were  discussed. 


(3)  hot  film  anemometry  -  Mork  by  RemkOt  Bremhorst 
and  Gilmore,  and  Lyon  were  discussed. 

(4)  film  thickness  measurements  capabilities  and 
limitations  of  impedence  probe  and  x-ray  technique 
were  discussed. 

(5)  optical  probe  -  current  work  at  England  and  at 
Grenoble  (CENG  and  AID)  were  discusser .  At  Grenoble, 
studies  are  being  carried  out  on  threshold  adjustfpent; 
on  in-bundle  void  measurements;  and  on  monofibre 
probes  for  droplet  flow. 

Z.  Radioactive  beams  -  S.  Banerjee.  McHaster  University 

Current  work  on  attenuation  and  scattering  methods  were 
discussed.  Particular  attention  was  given  to  multi -beam 
Y-ray  and  to  neutron  scattering  (by  B&M,  French  and  ScMaster). 
The  neutron  scattering  technique  is  a  promising  technique  for 
density  measureinent,  since  It  Is  Independent  of  flow  fattern. 
Only  difficulty  Is  the  low  counting  rate,  thus  not  very  good 
for  fast  transient. 

3.  Drag  body,  turbine  and  other  Impact  devices  -  L.  Davis, 
Measurement  Inc.  The  presentation  Includes: 

a.  Parameters  to  be  measured  and  sensor  requirements. 

b.  Current  sensors  used. 

c.  Methods  of  ccrbinlnn  sensor  readings  for  flow  parameters. 

d.  Current  development  to  Improve  sensor  -  work  Is  being 
done  at  Semiscale,  LOFT,  Germany  ("True  Flow  Meter"), 
Battelle,  M.I.  and  C-E. 

e.  Check  list  for  considerations  for  sensor  applications. 

4.  Tracer  techniques  -  P.  Kehler,  ANL 
The  presentation  Includes: 

a.  Outline  of  NRC-ANL  tracer  technique  program. 

b.  Discussion  of  basic  principles. 

c.  Lists  of  Y-em1tHng  tracers  for  gas;  for  liquid; 
commercially  available  radioisotope  milking  systems; 
elements  easily  activated  by  thermal  neutrons;  and 
those  activated  by  fast  neutrons. 

5.  Novel  techniques  -  P.  Griffith,  HIT 

a.  Unresolved  thermal -hydraulic  measurements  needs  were 
listed  as: 
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(1)  code  verification  -  subcfiannel  void  and  flow 
rate;  cnannel  or  pipe  void  and  flow  rate, 

(2)  diagnostic  measurements  -  flow  behavior  In  upper 
and  lower  plena  and  subchannel  flow  distribution 
In  core. 

b.  Novel  techniques  to  be  used: 

(1)  flow  measurements 

(a)  nuclear  magnetic  resonance  (no  disturbance 
but  need  flow  modeling), 

(b)  Laser-Ooppler  anemomety  (non-disturbing  and 
quick  response,  need  clear  line  of  sights, 
may  be  good  for  reflood,  but  not  very  SMitable 
for  Integral  tests). 

(2)  void  fraction 

(a)  x-ray  (flow  regime  dependent) 

(b)  Impedence  d<^v1ce  (flow  regime  dependent  but 
with  good  transient  response) 

(c)  acoustic  method  (flow  regime  and  qualUy 
dependent) . 

(3)  diagnostic 

(a)  Storz  lense  (good  viewing  angle,  needs  light, 
not  good  for  high  temperature,  good  diagnostic 
tool ) . 

c.  Conclusion 

(1)  none  of  the  novel  devices  1s  really  readily  useful. 
Impedence  probe  and  Storz  lense  are  worthwhile 
development, 

(2)  meanwhile,  drag  body-turbine,  y-densitometer  Is 
still  a  good  tool  for  pipe  flow  measurement, 

(3)  for  plena  flow  behavior,  photographic  methods  are 
the  only  ones  that  are  useful . 

III.  Open  Discussions  (January  K,  AM) 

1.  Tracer  technique  - 

Participants  were  concerned  about  the  following  difficulties 
Involved  In  tracer  measurements: 

a.  Diffusion  and  spread  of  tracer. 

b.  Injection  method. 

c.  Length  of  transient  time  vs  transit  time. 

1.  Presence  of  tracer  material  affecting  other  instrumentation. 

e.  Number  and  location  of  detectors. 

f.  Need  of  bench  tests. 


P.  Kahler  answe.'ed  sonc  of  the  concerns,  such  as  a,  b,  c, 
and  e,  and  tM>k  note  of  the  various  cautions. 

Pertinent  references  produced  are: 

(1)  P.  Kehler,  ANL-CT-76-17. 

(2)  Gopal,  N.Sc.  Thesis,  Penn  State  Univ.,  19o>  .n 
flow  meter. 

(3)  D.  Rowe,  BNWL-371-Pt.  3,  ’969. 

(4)  Quant's  work. 

(5)  HcHaster's  work. 

2.  Measurements  of  aP  diametrically  across  core  barrel  - 
(brought  up  by  W.  Hodges,  NRR) 

Range  of  Interest  200  psi.  In  1  millisecond.  Thu'  high 
frequency  fluctuations  are  difficult  to  measure.  Piezo¬ 
electric  transducers  were  used  in  Semiscale  but  tne  potting 
deteriorates  at  600°F.  Diaphram  Is  too  thick  for  response. 
Strain  gauge  on  the  core  barrel  was  proposed  as  one  possible 
answer.  It  may  need  cooling  to  eliminate  thermal  effect. 
Brockett  and  Davis  offered  to  help  Hodges. 

3.  Transit  time 

Last  year,  INEL  reported  sw.ae  Inconclusive  results  based  jpon 
cross-correlation  of  thermocouple  signature*-.  Raptis  o.  ANL 
suggested  that  a  more  sophisticated  approac.  is  that  of  Input 
response  function  (Reference,  Raptis  -  ANL-CT-76-39).  The 
consensus  Is  that  while  cross-correlation  method  should  be 
continuously  studied,  people  should  also  look  into  input- 
response  function  method. 

IV.  Two-Pha^e  Flow  Calibratit-.i  (January  14,  late  AM,  all  PM) 

1.  State-of-art  review  -  M.  Stanley,  INEL 

a.  Statement  of  problems 

(1)  definition  of  calibration, 

(2)  parameters  that  can  be  measured, 

(3)  parameters  that  are  difficult  to  n»easure. 

b.  Survey  of  existing  facilities 

(1)  Karlsruhe  loop, 

(2)  else  loop, 

(3)  WCL  loop. 
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(4)  French  loops. 

(5)  INEL  steedy  state  loop* 

(6)  INEL  transient  flow. 

Characteristics  and  measurement  capability  of  each 
facility  was  given. 

c.  Calibration  methods  currently  used  at  INEL  for  Semiscale 
and  LOFT  was  described. 

d.  Recommendations 

(1)  establish  a  universal  standard. 

(2)  establish  repeatable  conditions. 

(3)  for  phase  velocity,  use  tracer  tech'^lque  and 
local  probes. 

(4)  for  flow  regime,  use  multibeam  Y-densItoweters 
and  local  probes. 

(5)  for  density  distribution  use  multibeam  Y*dens1tometer. 


2.  Survey  of  facilities  for  two>phase  Instrumentation  testing  - 
H.  Estrada.  HPR 

Presentation  Includes: 

a.  Categories  of  testing 

(1)  conventlonttallbratlon  tests  >  to  determine  "iqeter 
factors"  In  each  single  phase. 

(2)  prototype  calibration  «  test  In  two-phase  flow  In 
conjunction  with  other  related  Instrumentation  to 
arrive  at  a  set  of  flow  parametenwith  the  help 
of  proper  modeling.  This  Is  the  "calibration" 
generally  referred  In  two-phase  flow  measurements 
community, 

(3)  mechanical  development  and  proof  testing. 

(4)  scaling  tests. 

(5)  transient  tests. 

b.  Facility  requirements 

(1)  conventional  calibration  (metering  factors)  - 
single  phase  loops,  full  flow, 

(2)  prototype  calibration 

(a)  gas-water  loop, 

(b)  geometry  effect, 

(c)  scaling. 

(3)  proof  testing  -  simulation  of  full  range  of  conditions 
to  test  Instrumentation  survivability.  May  be  skipped 
If  structurally  strong. 
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(4)  sealing  test  -  Urge  scaled  gas-water  systea 
properly  Instruaentcd,  with  visual  observation. 
Density  ratio  Mist  be  slMilated* 

(5)  transient  test  facilities  -  need  a  known  standard 
fast  response  neasuring  device. 

c.  Facility  survey 

Salient  finding  Is  the  need  of  a  high-flow  high  aonwitun 
flux  air-water  loop. 

3.  Discussion  on  standard  calibration  techniques  and  reference 
Instruments. 

The  discussion  on  this  topic  Is  lively  and  soir«t1ines 
controversial,  many  times,  opposing  views  were  voiced. 

The  following  Is  a  sunmary  drawn  up  by  G.  D.  Lassahn  of 
INEL.  with  additions/corrections  made  by  Y.  Y.  Hsu. 


Discussion  on  Standard  Calibration  Techniques 
and  Reference  Instruments 

Sumnary  Notes 

The  discussion  addressed  the  general  problem  of  calibrating  Instruments 
for  use  In  measuring  steam-water  flow  parameters,  with  partlr**1ar 
attention  to  establishing  a  standardized  system  of  callbratlr'  for 
use  throughout  the  Industry. 

General  Considerations 

Instruments  must  be  calibrated  Integrated  systems.  It  Is  usually 
not  acceptable  to  calibrate  an  Instnmient  separated  from  Its  companion 
Instruments.  Effects  of  environment  should  be  checked.  Upstream 
flow  perturbations  and  pipe  size  effects  must  be  realistically 
represented. 

Reference  Instruments  should  perturb  the  fluid  conditions  as  little 
as  possible. 

There  should  be  redundancy  In  the  reference  measurements. 

The  calibration  facility  should  be  readily  accessible  (physically 
close). 
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Concenius;  Host  calibration  should  bo  done  In  gas-water,  with 
occasional  checks  In  steasHwater. 

Dissenting  Opinion:  For  some  Instrunents  (drag  bodies  and  some 
transit  time  flowneters  for  example),  there  are  qualitative  or 
quantitative  performance  differences  In  steam-water  and  In  gas- 
water. 

Cownents; 

The  main  problem  with  steam-water  loops  Is  In  knowing  exactly  what 
the  flow  conditions  are.  Freon-water  loops  avoid  condensation 
effects,  but  turbine  bearings  do  not  survive  long  In  freon. 

Air-water  loops  may  be  more  expensive  than  steam-water  loops  with 
similar  flow  rates,  depending  on  the  particular  Installation. 

An  operating  pressure  of  about  2.8  MPa  (400  psi)  In  an  air-water 
loop  may  be  sufficient  to  simulate  a  15  MPa  (2200  psi)  operating 
pressure  In  a  steam-water  loop. 

Gas-water  loops  have  the  advantage  that  the  gas  and  liquid  flow  rates 
can  be  measured  separately. 

There  Is  no  mass  transfer  between  the  phases  In  a  gas-water  loop; 
tnis  Is  an  advantage  In  that  It  slBq>11f1es  a  detailed  study  of  the 
Instrument  performance,  but  It  Is  a  disadvantage  in  that  It  may  give 
•Jlfferent  calibration  constants  for  some  instruments  (different  from 
steam-water  calibrations,  due  to  difference  In  density  ratios). 

Steady  State  Calibration  of  Full  Flow  or  Global  Measurements 

Consensus;  Reference  Instruments  could  be  a  pair  of  full  flow  drag 
screens,  for  pV^  and  transit  time  measurements,  a  scanning  ganna 
densitometer,  for  density  average  and  perhaps  density  profile,  and 
tracers  as  a  backup,  to  measure  either  total  mass  flux  or  separate 
phase  mass  fluxes. 


Pissentlnq  Opinion:  A  drag  screen  upstreaa  froM  the  test  Instrument 
tends  to  homogenize  the  flow  so  that  the  Instrument  does  not  see 
the  full  range  of  flow  regimes;  a  drag  screen  downstream  does  not 
measure  the  momentum  flux  Incident  on  the  test  instrument. 

We  cannot  measure  density  profiles  reliably  with  presently  available 
scanning  gamma  densitometer. 

Tracers  (either  Injected  or  neutron  activated)  are  reliable  and 
meaningful  global  reference  measuret<ient.  and  should  be  used  as  the 
primary  standard. 

Comnents; 

Photographic  techniques  are  a  powerful  tool  In  at  least  some 
applications. 

Neutron  activation  and  radioactive  tracers  are  expensive  and  not 
universally  available. 

Steady  State  Calibration  of  Local  Measurements 


Concensus;  The  reference  Instruments  could  be  retractable  local 
probes  (of  unspecified  type),  checked  by  photographic  techniques. 
Scanning  gaima  densitometers  might  also  be  useful. 

Dissenting  Opinion:  Most  local  probes  Interfere  with  th_-  flow. 
Paragraph  photographic  techniques  are  not  useful  In  some  flow 
regimes. 

Local  measurements  should  be  used  only  as  a  diagnostic  tool; 
accurate  calibration  Is  not  necessary. 

Ccmments: 


Isokinetic  probes  do  not  work  well  for  the  gaseous  phase,  and  they 
perturb  the  flow  too  much. 

Laser>Doppler  shift  anemometry  Is  useful  In  certain  restricted  flow 
regimes. 
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TiNHcelor  lllui. nation  phetograplilc  taehnlquos  holp  ovorcot  tht 
probicas  usually  ancountart4  with  standard  photographic  taehnlquos. 

Flash  radiography  and  gaaaa  radiation  Frasnal  holography  Might  hava 
some  potential. 

Hot  wires  alght  ba  usaful  In  dataralnlng  tha  liquid  phase  local 
velocity. 

Transit  tiae  flowaatars  should  ba  usaful  In  local  velocity  aeasurcaents* 
although  they  alght  causa  local  velocity  perturbations. 

Steady  State  Calibration  of  Rakes 

Concensus;  The  separata  probes  of  a  rake  can  first  be  calibrated 
separately  but  the  whole  rake  asscably  Must  be  callurated  with 
piping  and  flow  conditions  which  realistically  represent  the  Intended 
applluatlon  of  the  rake. 

Dissenting  Opinion;  In  soae  opinions,  separate  free  field  calibration 
of  the  Individual  probes  of  a  rake  Is  sufficient.  But  soae  others 
believe  that  the  full  rake  calibration  is  the  only  Meaningful  Method. 

CoMMents; 


Full  scale*  realistic  geoMetry  tests  are  necessary  to  determine  optInuH 
rake  geoaetry. 

All  flow  reglaes  Must  be  represented  In  calibrating  a  rake  assembly. 

Large  rakes  are  being  considered  cnly  for  LOFT.  In  LOFT*  flow  regime 
effects  are  not  Important  until  the  latter  part  of  the  blowdown,  when 
the  pressure  Is  relatively  low.  It  would  be  much  easier  to  calibrate 
the  whole  rake  asseably  only  at  these  lower  pressures. 

Steady  State  Calibration  Checks  In  Steara»Water 

Cownents; 


The  total  aass  flux  through  the  pipe  Is  not  affected  by  phase  changes. 

The  steaa  and  water  should  be  near  thermal  equilibrium  when  they  are 
Injected  Into  the  mixer. 


Scaling  (pipe  size)  effects  are  probably  different  for  stean-«Mter 
and  a1r>Mater. 


Transient  Calibration 


Concensus;  The  major  purpose  of  studying  the  transient  response  of 
an  Instrument  Is  to  know  the  magnitude  of  the  error  of  the  Instrument 
readings  (as  calibrated  In  steady  state)  In  transient  flow. 

It  Is  very  difficult  to  create  a  known  fast  transient  In  steam-water 
flow.  Because  of  this  difficulty,  transient  performance  experiments 
should  be  done  In  single-phase  fluids  or  perhaps  In  air-water,  even 
though  the  Instrument  performance  may  be  different  In  air-water  than 
In  steam-water.  Good  mathematical  modeling  Is  very  Important  In  studying 
the  transient  response  of  an  Instrument. 

Comments ; 


If  the  transient  response  of  a  device  Is  linear  and  Independent  of 
flow  conditions,  then  the  transient  response  expression  can  be 
Inverted  to  obtain  the  real  flow  parameter  from  the  Instrument  readings. 
Turbine  flowmeters,  for  example,  do  not  meet  these  requirements. 


The  only  reliable  calibration  check  we  have  in  transient  flow  Is  the 
total  mass  transport  during  the  transient. 


Scaling  or  Pipe  Size  Effects 

Concensus;  Some  full-flow  Instruments  (a  full  flow  drag  screen,  for 
example)  must  be  calibrated  in  a  full-sized  pipe.  For  free  flow  rake, 
full-size  pipe  Is  not  necessary,  provided  flow  pattern  and  flow  profiles 
are  properly  modeled  and  measured. 

To  assess  pipe  size  effect  v«hich  has  been  observed  In  several  experiments 
calibration  snould  be  done  with  increasing  pipe  size  until  pipe  size 
effect  levels  off.  Full-sizes  pipe  and  realistic  simulation  of  upstream 
geometry  are  necessary  in  at  least  some  calibration  work.  It  is 
important  to  stress  that,  whether  for  calibration  or  for  actual  measure¬ 
ment  application,  the  upstream  geometry  effect  should  be  properly 
modeled  and  accounted  for.  To  achieve  this  objective,  the  selection 
of  instruments;  their  location;  deployment;  and  alignment;  should  be 
carefully  done  with  the  help  of  knowledge  In  flow  pattern  and  flow 
distribution. 


G«naril  Conclusions 


Thfs  problea  of  standardizing  calibration  techniques  Is  linportant 
and  needs  further  study.  Or.  Y.  Y.  Hsu  will  establish  a  task  force 
to  work  on  this  problem. 
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A6EN0A 

January  13.  1977,  Thursday 
AM  >  Progress  Reports  by  NP.C  Contractors 
8:30  •*  Uelcome  and  Introduction,  Y.  Y.  Hsu 

8:35  -  Progress  Report  by  ORNL.  J.  Sheppard 

9:15  *  Progress  Report  by  RPI  (IR,  RF,  X-ray),  R.  Lahey 

9:d5  '  Progress  Report  by  BNL  (local  probes),  0.  Jones 

10:05  -  Progress  Report  by  SUKY  (LDV),  R.  Lee 

1C:25  -  Break 

10:35  •  Progress  Report  by  INEL,  Staff  of  EGAG 

12:40  -  LUNCH 

PM  -  State-of-art  Survey  by  Invited  Rapoorteurs  -  (Significant 
developront  In  design,  operation,  modeling  and  Interpretation 
of  tMO-phase  Instrumentations:  with  emphasis  In  progress  In 
1976  •  30  minute  presentation,  10  minutes  discussion,  each) 

1:30  -  Local  probes,  (hot-wire  anemometer,  Impedence  probes, 

optical  probes,  mlcrothermocouple,  etc.),  J.  K.  Delhaye, 
Grenoble 

2:10  -  Radioactive  beams  (attenuation  or  scattering), 

S.  Benerjee,  McMaster  University 

2:50  -  Break 

3:00  -  Drag  body,  turbine  and  other  Impact-activated  devices, 

L.  Davis,  ^'easurement,  Inc. 

3:40  -  Tracer  technique,  P.  Rehler,  ANL 

4:20  -  "ovel  techniques  (acoustic,  PMR,  IP.,  RF,  X-ray, 

Auburn,  fluidic,  Storz  lense,  LW,  etc.),  P.  Griffith, 
Kassachusett  Institute  of  Technology 

5:00  -  Adjourn 
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J«iHiT»  14.  1977.  Friday 
M  •  Discussion  and  Presentations  - 
8:30  -  Discussion  on  papers  presented  January  13,  1977 

10:00  >  Break 

10:10  >  State-of*art  report  on  calibration  tecbr.ttpies  and 

procedures.  M.  Stanley,  EG&G 

11:10  -  MPR  report  on  requirements  and  procedures  for  calibration, 

(MPR  report  will  be  mailed  to  the  addresses  prior  to  the 
meeting  to  serve  as  frMeowrk  for  discussion  In  the 


•  •  wci  lawwaa/ t  w*  «*w«a  e 


12:00 


LUNCH 


PK  -  Discussion  on  Calibration  Procedures  and  Techniques  - 
1:00  -  Discussion  -  Following  questions  will  be  discussed 

1)  How  to  calibrate  full-field,  air-water,  steady 
state  instruments  In  small  channel? 

2)  Now  to  calibrate  free-fleld  Instrumentation? 

3)  How  to  scale  up? 

4)  Is  It  necessary  to  calibrate  In  transient  conditions? 
If  so,  any  difficulties? 

5)  Is  It  necessary  to  calibrate  In  steam-water  mixutre? 
If  so,  any  difficulties? 

6}  Should  (Item:,  3,  4)  (Items  4,  5)  (Items  3,  5)  be 
combined? 


7)  Should  Items  3,  4,  S  be  all  combined,  I.e.,  Is  It 
necessary  to  calibrate  Instruments  In  transient, 
steam-water,  large  loop? 

4;30  -  Summary 

S:00  -  Adjourn 
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2:10  -  Radioactive  beams  (attenuation  or  scattering), 

S.  Benerjee,  HcMaster  University 
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nOGBXSS  REPORT  ON  AOVAMCEO  lUO-PHASE  IRSHUNBIIMTUM 


J.  D.  Shspparl 

W.  H.  LmvcII  F.  Shahrokhi 
M.  C.  Uynu 


ABSTRACT 

A  Venturi  meter  was  teated  aa  a  two-phase  momentum  moni¬ 
toring  device  in  the  ORNL  air-wat'r,  two-phase  flow  facility. 
Void  fraction  was  monitored  at  the  Venturi  inlet  with  a  ro¬ 
tating  electric  field,  conductance  meter,  and  Venturi  pressure 
drop  was  correlated  with  mixture  momentum  determined  by  a 
separated-flow  model. 

A  computer  algorithm  was  developed  to  apply  cross-corre-  / 
latlon  techniques  to  the  measurement  of  nonstationary  two-  / 
phase  velocities.  The  algorithm  was  tested  with  computer- 
generated,  two-phase  velocity  transients.  Preliminary  reaults 
show  good  agreement  between  algorithm  estiautes  and  the  known 
velocities. 


1.  ADVANCED  SPOOL  PIECE  STUDIES 
J.  D.  Sheppard  M.  C.  Wynn 

The  determinatiwi.  .?  transient,  two-ohase  mass  flow  rrtes  during 
blowdown  heat  transft  (poriments  generally  requires  the  combination 
of  two  or  more  characteristic  fluid  parameters  (e.g.,  velocity  and 
density)  with  appropriate  two-phase  flow  models.  For  measurements  in 
the  primary  piping  of  experimental  loops,  flow-monitoring  instruments 
have  been  concentrated  in  piping  spool  pieces  located  at  selected  sites 
in  the  loop.  Spool  pieces  at  the  (MNL  Thermal -Hydraulic  Test  Facility 
(THTF)^  and  the  INEL  Semiscale  Facility^  utilize  full  pipe  turbine 
meters  for  mixture  velocity  measurement,  gamata  densitometers  for 
mixtuie  density,  drag  disks  for  mixture  momentum,  and  te^>erature  and 
pressure  transducers  for  determinations  of  liquid  and  vapor  phase  states. 


1. 1-1 


Th«  scop*  of  th*  Advoncod  tWo«Phos*  Instniaontoclon  Protra*  includes  con* 
sideration  of  new  or  alternative  instruaents,  or  conbinations  thereof, 
for  what  is  referred  to  as  an  advanced  spool  piece  in  order  to  laprove 
the  accuracy  and  precision  of  transient  aass  flow  rate  deterainations. 

1.1  Two-fhase  Moaentua 

Many  nuclear  reactor  safety  studies  in  this  country  utilize  drag 
disks  (target  flowmeters]  for  two-phase  aoaentua  aeasurtaents.  Drag 
disks  are  capable  of  indicating  direction  as  well  as  aagnitude  of 
aoaentua  (i.e.,  bidirectional];  however,  they  have  the  undesirable 
characteristic  of  saapling  a  relatively  saall  fraction  of  the  flow 
cross-sectional  area.  One  solution  to  this  problea  is  to  eliainate  the 
disk  target  and  utilize  a  screen  target^  that  saaples  more  of  the  flow 
area.  However,  other  aoaentun  aeasureaent  techniques  are  applicable  and 
should  be  considered.  For  example,  tt  phase  aoaentua  has  been  related 
to  the  pressure  drop  of  resistance  C!  .?ents  such  as  woven  wire  and  per¬ 
forated  plate  screens.**  This  report  presents  soae  preiiainary  results 
of  steady-state,  air-water  studies  with  a  Venturi  meter  and  the  ap.  li- 
cation  of  a  separated-flow  model  to  data  correlation. 

Steady-state,  two-phase  flow  tests  were  conducted  in  the  air-water 
loop  with  a  Venturi  .seter  installed  for  vertical  upflow.  The  Venturi 
had  a  major  inside  diameter  of  7.62  cm  (3.0  in.)  and  a  3.81-ca  (l.S-in.) 
throat  diameter;  the  total  angle  of  convergence  upstream  of  the  throat 
was  21*,  and  the  total  angle  of  divergence  downstream  of  the  throat  was 
7*.  The  pressure  drop  was  measured  with  differential-pressure  trans¬ 
ducers.  Due  to  the  fluctuating  character  of  the  transducer  signal  during 


t«o»pii«M  tMting*  M  intttrcting  digital  iraltaatar  used  to  mrarag*,^ 
th«  signal  for  poriods  of  approxinatoly  1  nin.^ 

VoluMt-avorago  void  fraction  was  nonitorod  iModiatoly  t^stroaa  of 
tho  Vonturi  by  a  rotating  oloctric  field  conductance  notor^  (Auburn 
International.  Inc.,  nodal  1080).  Further,  a  two-phase  flow  disperser 
consisting  of  a  stack  of  four.  20-aesh  woven  wire  screens  was  located 
at  the  upstroaa  flange  of  the  void  fraction  neter  to  nininise  radial 
variations  in  the  two-phase  velocity  distribution. 

Venturi  pressure  dr«9  and  void  fracti<m  wore  nonitorod  as  a  function 
of  air  and  water  flow  rates.  Axial  separation  of  the  vertically  oriented 
pressure  taps  was  10.2  cn  (4  in.),  so  the  neasured  pressure  drop  included 
gravitational  effects  (ogh)  which  were  considered  in  the  data  reduction. 

F' the  tests  with  wnter  only,  the  gravitaticnal  tern  was  1  kN/n^  (0.14 
f»si). 

Figure  1.1  shows  the  pressure  dr^.  corrected  for  the  gravitational 
tern,  as  a  function  of  sir  flows  fron  0  to  0.24  nVsec  (0  to  S12  scfn) 
and  water  flows  fron  6.3  «  10***  to  2.S  »  lO'^  nVsec  (10  to  400  gpn). 

As  shown  in  the  figure,  the  gravitational  effects  were  of  the  saaw  Magni¬ 
tude  as  the  Bernoulli  effects  only  at  the  lower  water  and  air  flows.  At 
the  highest  air  flows  (0.24  n^/sec).  the  water  flow  rate  was  United  to 
‘vS  e  10'^  n’/sec;  at  higher  water  flows  the  total  systen  pressure  drop 
was  too  high  to  naintain  critical  flow  in  the  air  neters.  Venturi 
pressure  drqp  ranged  fron  'vl  to  200  kN/n^. 

A  separated-flow  nodal  was  used  by  Collins  and  Gacesa*  and  nodified 
by  Fouda  and  Rhodes^  to  detemine  two-phase  quality  fron  Venturi  &P,  nass 
flow,  and  Venturi  geonetry.  While  their  nodel  assumes  no  slip  between 


I  pkuM,  it  coMidm  Mck  •ffwts  at  thadiffarMca  ia  iacfprtittblr  aai 

ctapfattiblt  flow  through  tho  Vouturi.  Tho  hya*  t«parato4>flou  aodol 
coatldort  tho  offoctt  of  slip  by  atsitaag  that  oaeh  of  tho  two  phasos 
caa  bo  r^rosaatod  with  a  charactoristic  volocity,  or  aad  a  givan 
iastruMat  raadiag  is  tho  rosponsa  of  tha  iastnaant  to  aach  phasa  of  tka 
fluid  aiatura.  This  aadal  has  rosultad  ia  distinct  iiyrewsaant  ovar  aora 
caaaon  houogaaoous  aodals  in  datarainiag  aass  flow  frou  IHrP  spool  pioca 
iastnsMBt  rosdiags.* 

According  to  tha  Ays  aodal.  two-phasa  noaantua  is  givan  by 

Cad  -  a)OtVj^  .  d-1) 


wharo  and  ara  aaan  valocitias  of  tha  vapor  and  licpiid  phasas,  rospac 
tively,  0|  and  are  phase  densities,  and  a  is  void  fraction. 

For  tha  data  shown  ia  Fig.  1.1,  void  fraction  was  naasured  near  the 
Venturi,  aad  individual  phases  wore  naterad  into  the  systos;  therefore, 
phasa-cantinuity  relationships  pemit  the  calculation  of  the  phase  veloc¬ 
ities,  aad  V^,  ontaring  tha  Venturi.  That  is. 


(1.2) 


•  (I  -  ®)pjV^  .  (1.3) 

where  a^  and  are  the  swaa  flows  of  air  and  water,  respectively.  After 
solving  for  the  phase  vslocltlas,  the  ulxture  nouentun  can  be  calculated 
fron  tq.  (1.1). 


I.X-5 


Picurt  1.2  skOMS  tk*  VMturt  pr««suM  drop  m  •  fuactlen  of  tko  fluid 
mtuB  dotoxuiaod  kjr  Eqs.  (1.1)  through  (l.S),  oloag  with  void  fraction 
phaso  aast  flow  rata  data.  Tko  faaily  of  curvos  shown  in  Fig.  1.1 


roducod  to  a  singlo  band  of  data,  with  (AF  —  pgh)  varying  linaarly  with 
pV^|,  for  noarly  throo  orders  of  aagnltudo.  Tho  width  of  the  scatter  band, 
tSO%,  can  probably  be  significantly  reduced  by  inproved  veld  fraction 
■easursaents  since  the  void  fraction  aeter  is  sensitive  to  flow  regiae. 

Figure  1.2  and  Eq.  (1.1)  give  Che  relationship  between  pressure  drop 
and  aiature  aoaentua  for  this  Venturi  as 

AP  -  e^gh  -  Clop^V*  +  (1  -  a)p^V*J  ,  (1.4) 

where  C  is  a  constant.  Rewriting  Eq.  (1.4)  gives 

AP  -  p^gh  -  Cavjlp^s'  ♦  pjd  -  a)/al  ,  (1.5) 

where  S  i4  the  slip  ratio,  V^/V^.  If  the  Venturi  reading  is  coabined  with 
turbine  aeter  and  gaaaa  densitoaeter  data,  aass  flow  and  quality  can  be 
calculated  froa  the  Aya  aodel.  However,  if  only  Venturi  pressure  drop 
and  void  fraction  data  are  avail- ble,  an  assuaption  of  the  value  of  the 
slip  ratio  is  required  to  give  aass  flow  rates.  Due  to  the  fact  that  p^ 
is  auch  less  than  for  air-water  and  steaa-water  aixtures,  the  actual 
value  of  the  slip  ratio  is  not  iaportsnt  in  Eq.  (1.5)  except  at  very  high 
values  of  void  fraction.  For  exaaple,  for  air  and  water  at  aabient  condi¬ 
tions,  p  ■  0.07  lb  /ft^  and  p.  ■  62  lb  /ft^.  If  the  void  fraction  were 

g  II  t  m 

as  high  as  95%,  the  slip  ratio  would  have  to  be  nearly  50  for  the  g^s-phase 
effects  to  be  of  the  saae  aagnitude  as  the  liquid-phase  effects.  As  dis¬ 
cussed  above,  Fouda  and  Rhodes^  assumd  unit  slip  and  achieved  good  corre¬ 
lation  for  3  wide  range  of  stean-water  flows. 
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TWe-phas*  Vanturi  studies  will  centiiuM  with  eaphssis  on  optinus 
throat  design  for  two-phase  flow  and  the  use  of  void  fraction  data  deter- 
nined  from  a  triple-bean  gaana  densitoeMter.  Further,  many  blowdown 
experinents  include  flow  reversals,  and  Venturi  designs  that  will 
accowodate  flow  raversal  will  be  considered. 

1.2  Fluidics 

A  subcontract  has  been  initiated  with  Harry  Diaaond  Laboratories 
(HDL),  Departnent  of  the  Amy,  to  investigate  the  application  of  fluidics 
technology  to  two-phase  flow  neasurenent  problens.  The  scope  of  the  work 
to  be  perforaed  by  HDL  includes  (1)  a  study  of  the  feasibility  of  using 
fluidics  technology  to  aonitor  mss  flow  rate,  nixture  quality,  and  other 
key  transient  two-phase  flow  parameters;  and  (2)  the  fabrication  of 
devices  denonstrating  the  application  of  fluidics  principles  to  the 
neasurenent  of  (1)  two-phase  noaentun,  (2)  taaperature,  rnd  (S)  pressure 
drop.  These  devices  will  be  tested  in  ORNL  two-phase  flow  facilities. 
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.  nmstmanuei  stem.  comtEUTioN 

N.  H.  Lwv«ll*  F.  Shahrokhi* 

2.1  Thaery 

On*  tachnii^  of  fluid  velocity  Moturoaeot  is  to  directly  Masure 
the  traasit  tiae  of  a  quantity  betueea  two  points  A  and  B,  fixed  a  dis¬ 
tance  0  apart  la  space,  and  thea  divide  Che  distance  D  by  the  transit 
tins  to  calculate  Che  velocity.  If  the  asasured  process  is  randea,  the 
statistical  technique  of  cross-correlation  is  often  used  to  neasure 
transit  tiae.  Usinc  this  technique,  sifnals  A(t)  and  B(t)  are  asasured 
for  T  seconds  at  points  A  and  B,  respectively,  and  their  cross-correlation 
function  is  approxiaated  using  the  following  general  expression: 

f  jT  A(t)  B(t  ♦  T)  dt  .  (2.1) 

where  t  is  the  value  of  tiae  delay  between  points  A  and  B.  The  transit 
tiae  between  points  A  and  ■  of  the  asasured  process  will  correspond  to 
a  unique  value,  t^,  of  t  that  aaxiaites  the  cross-correlation  function 
given  in  Eq.  (2.1).*  If  the  asasured  signals  A(t)  and  B(t)  are  contaai- 
nated  by  additive  noise,  tiae-averaging  of  h^(T)  will  iqirove  the 
signal-to-noise  ratio  if  the  noise  is  uncorrelated  because  the  cross 
products  between  the  noise  and  the  signals  in  Eq.  (2.1)  will  average 
tero.  An  inportant  assuagttion  of  this  averaeina  process  is  that  the 
fluid  velocity  being  asasured  is  constant  throughout  the  tiae-averaging 
interval.  If  the  fluid  velocity  of  the  process  is  not  constant,  as  in 
the  case  of  a  blowdown  test,  conventional  cross-correlation  analysis 
• 
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is  no  loafor  spplicsbls  to  ths  problcn  snd  a  nedifis^  procedurs  is 
rstpiirsd. 

If  ths  ntssurenent  procedures  were  not  altered  for  a  nonstationary 
fluid  velocity,  what  would  be  the  eapected  results  if  the  cross-correla¬ 
tion  function  (Eq.  (2.1)]  were  calculated?  First,  assune  for  sinplicity 
that  the  fluid  velocity  decreases  with  tine  at  a  constant  rate  (in 
practice,  the  algorithm  does  not  require  that  this  assumption  be  true) 
and  that  the  measurement  at  points  A  and  B  of  the  random  process  proceeds 
for  T  seconds.  Daring  those  T  seconds,  the  fluid  velocity  has  changed 
from  a  maximum  value  at  t  •  0  to  a  minimum  value  at  t  *  T.  If  the 
cross-correlation  function  (Eq.  (2.1))  were  not  coaq>uted  over  the 

entire  T  seconds  of  data  and  examined,  the  peak  value  of  Rj^(t)  corres¬ 
ponding  to  the  unique  transit  tine  in  the  constant-velocity  case 
would  appear  broadened  to  include  the  range  of  transit  tines  that  were 
momentarily  present  during  the  period  T;  the  minimum  transit  time  would 
correspond  to  the  maximum  velocity  at  t  «  0,  and  the  maximum  transit 
tine  would  correspond  to  the  minimum  velocity  at  t  >  T.  However,  if 
the  measurement  tine  T  could  be  nade  sufficiently  snail  so  that  the  flow 
velocity  could  be  considered  constant  (to  first  order)  over  the  interval 
0  to  T,  then  the  velocity  could  be  tracked  as  it  decreases  on  a  longer 
tine  scale  by  continually  readjusting  T  so  as  to  resolve  each  new  mini¬ 
mum  transit  tine  corresponding  to  the  changing  velocity.  Such  step¬ 


wise  alteration  of  the  neasurenent  tine  T  forms  the  basis  of  a  modified 
cross-correlation  algorithm,  which  is  presented  next. 


2.2  HpiuMt lowary  Cro<s-»Corf  litioii  Alioritlw 

TIm  croos'CorrolotiOR  altorltiM  that  has  Imoo  davalopod  to  track  a 
aowstationary  flow  velocity  caw  bo  dividod  into  four  wain  parts  (soo  Fig. 
2^):  initialising,  partitioning,  voraging.  and  data  shifting.  Tbo 
a5^ua|^iows  aado  in  tho  algoritbn  aro  (1)  tho  aajtinuw  nunbor  of  data 
points  oporated  on  by  tho  'algorltha  during  any  one  iteration  is  1024 
(this  is  only  a  linitation  of  tho  particular  conputer  core  site  used  and 
not  a  basic  linitation  of  the  algorithn);  (2)  tho  distance  D  between  the 
two  neasuring  points  is  fixed;  and  (3)  the  initial  fluid  velocity  V(0) 
is  known. 

The  initializatiou  portion  of  the  algorithn  (see  Fig.  2.1)  consists 
of  three  tasks:  (1)  choosing  an  initial  sanpling  rate  f^;  (2)  choosing 
an  initial  partition  P  of  the  neasureswnt  tine  T;  and  (3)  sanpling  the 
signal  A(t)  at  point  A  and  the  signal  B(t)  at  point  B  in  the  fluid  nodiun. 
In  choosing  the  initial  settling  rate  f^.  the  neasurenent  tine  T  is  fixed 
since  1/f^  ■  At.  the  tine  interval  between  sai^>led  data  points,  and  the 
nunber  of  data  points  is  fixed  at  1024;  thus. 

T  ■  10244t  -  1024/f^  .  (2.2) 

Further,  the  resolution  of  transit  tines  between  points  A  and  B  is 
United  to  lAt  seconds.  The  saiqiling  rate  is  chosen  based  on  two  criteria 
which  are  functions  of  the  initial  velocity  V(0).  The  tisw  interval  At 
nust  be  snaller  than  the  initial  transit  tine  between  points  A  and  B. 
T^(0)  »  0/V(0);  thus,  the  first  criterion  is  that  f^  nust  be  greater 
than  ^/'^^(O).  The  second  criterion  states  that  the  nininun  total  neas- 
urenent  tine  T  nust  be  sufficiently  large  that  t^(O)  can  be  adequately 
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resolved  using  the  cross>correl«tion  elgorltlsB.  This  is  ecco^;»lished 
by  defining  T  ss  soae  nultiple  of  ^^^(0)  can  be  resolved  using  the 
first  2S6  date  points  of  the  total  1024  points  saapled.  Therefore  the 
■ininuai  T  has  to  be  equal  to  or  greater  than  seconds,  and  accord¬ 

ing  to  Eq.  (2.1), 

f^  i  1024/4t^(0)  .  (2.3) 

Thus  the  upper  and  lower  bounds  of  f^  are 

1/Ta,(0)  «  f,  S.  256/t^,(0)  .  (2.4) 

Although  the  selection  of  f^  is  based  on  the  assunption  that  .'lO)  is  the 
Mxinua  velocity  neasured  during  the  T  seconds,  this  does  not  prevent 
velocities  greater  than  V(0)  to  be  resolved  if  they  occur  d’^ring  the  T 
seconds,  since  T  seconds  is  nore  tine  than  is  necessary  to  resolve 
higher  velocities  than  V(0).  As  an  initial  selection,  f^  is  norsMlly 
set  at  one-half  the  naxisMoi  saapling  rate,  or  128/t^(0).  Since  the 
sasqiling  nave  is  fixed  for  the  entire  set  of  1024  data  points,  variation 
of  the  neasurenent  tine  T  to  resolve  the  nininun  transit  tine  is  accon- 
plished  by  taking  fractions  of  the  total  tine  T  and  using  these  to  cal¬ 
culate  each  transit  tine  ^^3(1^)  fmn  the  set  of  transit  tines  covered 
by  the  total  tine  T. 

The  second  task  acconplished  during  initialization  is  the  selection 
of  an  initial  partition  P,  which  is  typically  set  equal  to  the  first  512 
data  points,  since  the  sanpling  rate  is  one-half  the  aaxinun  saapling 
rate.  Having  selected  f^  and  P,  the  algoritha  sanples  sensors  A  an^  B 
for  T  seconds,  yielding 

A  *  {ai,  aj,  .  .  .,  aio2<«)  (2.5) 
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B  ■  {b|,  b2,  .  .  ••  b|02i*)  * 


(2.6) 


ifhere  and  b^  are  sa^>led  values  of  A(t)  and  B(t).  the  sensor  signals. 


(2.7) 


The  partitions  of  A  and  B  have  the  fora 
P(A)  *  (aj,  a2(  •  •  "i  A5J2I 
and 

P(B)  »  (bj,  b2»  •  •  ..  b5|2)  .  (2.6) 

The  second  ujor  part  of  the  algorithm  is  the  partitioning  routine 
(see  Fig.  2.1).  Its  task  is  tc  deteraine  the  ainiaua  size  of  partition 


P  for  the  total  tine  T  such  that  only  the  ainiaua  transit  tine  between 


sensors  A  and  B  will  be  resolved  and  all  other  transit  tines  greater  than 
will  be  rejected.  Since  the  data  for  A  and  B  have  already  been  parti¬ 
tioned  by  P,  the  algorithm  will  only  consider  the  first  512  data  points  in 


A  and  B.  The  naxiaun  transit  time  t  that  can  be  resolved  by  the  »lgo- 

BMX 


ritha  is  initially  equal  to  (S12/4)At>  or  128&t  seconds.  By  the  pirti- 
tioning  of  A  and  B,  the  algorithm  has  already  rejected  any  transit  tine 
greater  than  128At. 

The  next  step  is  to  cross-correlate  the  partitioned  A  and  B  data 
using  the  following  equation: 


R 


i-j 


AB^^J^  "  (N  -  j)At  ^  *n‘’n  ^ 


(2.9) 


where  «  jAt  with  j  *  0,  128;  and  N  >  512,  the  initial  partition  size 


of  P.  The  reverse  cross-correlation  function  Ro»(t4)  is  now  searched 

DA  j 


for  a  naxiaun  value  such  that 

*^AB^'aB^  ”  MAX[R|^g(Tj)] ,  j  *  1,  t^a,/At 


aax' 


(2.10) 


This  determines  a  ainiaua  transit  time  for  the  initial  partition  size 


1.1-15 


th«  auii 


P  ■  S12  points.  The  valuo  is  now  eoaparod  to 
rasolvabla  transit  tiM  for  partition  P,  usint  tha  followinf  expression: 

'■ax"’AB-*T»  -S 

If  the  above  coaparison  is  true,  then  the  algoritha  has  found  the  ainiaua 
partition  site  P  equal  to  S12  points.  It  has  also  detensined  the  next 
■ininuai  transit  tiae  following  the  initial  ainiaua  transit 

tine  corresponding  to  V(0).  The  algorithn  would  not  proceed  to  the 
averaging  routine. 

If  the  difference  [Eq.  (2.11))  between  and  is  greater  than 
seconds,  then  there  exists  an  even  saaller  site  for  the  partition  P 
of  the  data  set,  and  the  calculated  for  P  ■  S12  points  is  not  the  true 
ainiaua  transit  tine  for  the  data  set.  The  algoritha  now  coaputes  a  new, 
saaller  sire  for  the  partition  P  bx  setting  it  equal  to  ^t^At  data 
points.  It  also  nakes  ■  t^,  the  previously  calculated  ainiaua 
transit  tiae.  The  algorithn  again  calculates  [Eq.  (2.10)),  but 

with  j  ■  0,  t^/At,  and  N  ■  ^t^/dt.  A  search  is  again  aade  of 
according  to  Eq.  (2.11)  and  a  new  ainiaua  transit  tine  is  found  and 
coapared  to  the  new  using  Eq.  (2.11).  Nith  each  interaction  of  the 
parcicionlng  rouclne,  the  larger  transit  tines  that  were  aoaentarlly  pres¬ 
ent  during  the  total  T  seconds  are  rejected  until  the  true  ainiaua  transit 
tine  for  the  T  seconds  of  data  is  deteiained  by  the  criterion  in  Eq.  (2.11) 
and  set  equal  to  At  this  point  in  the  routine,  the  partitionr.  of 

both  the  A  and  the  B  data  set  appear  as 

P(A)  ■  *2*  *3»  •  •  •»  (2.12) 
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r(i)  -  {b,.  ba.  b, . b^»  . 


whcrt  M  m  4t  /at. 


(2.11) 


Since  the  partitioni'ig  pert  of  the  elgoritha  hes  yielded  only  one 


estiaete  of  the  Binieun  delay  tine  the  statistical  precision  is  very 
poor  and  soae  for*  of  averaging  aust  be  performed  in  order  to  increase 
the  confidence  in  this  estiaete.  Tiae  averaging  of  cross>correlation 
functions  calculated  for  a  nonstationary  process  ordinarily  introduces 
large  amounts  of  error  in  the  estiaated  transit  tiae,  but  in  this  case 
the  algoritha  can  orily  estiaate  the  values  of  the  transit  tiae  T^,(t)  in 
fixed  increaents  of  At  seconds,  the  tiae  betueen  data  saaples.  Thus,  the 
nuaber  of  tiae  averages  k  that  the  algoritha  can  perfora  without  intro¬ 
ducing  any  aore  error  in  the  transit  tiae  than  At  is  a  function  of  the 
resolution  At  and  the  rate  of  change  of  the  transit  tiae 
respect  to  tiae.  The  averaging  routine  of  the  algoritha  (see  Fig.  2.1) 
cakes  advantage  of  the  At  resolution  of  fgg(t)  !>/  u»ing  data  shifting  as 
a  substitute  for  tine  averaging  in  the  following  way.  Starting  with  the 
ainiaua  size  of  the  partition  P  resolved  by  the  partitioning  routine, 
the  algoritha  shifts  this  partition  by  one  At.  The  shifted  partitions 


<'f  the  A  and  B  data  sets  now  appear  as 
P(A)  ■  (*1*  *3»  *  *  *•  ^  ♦ 


(2.14) 


P(B)  -  (b,,  b2.  bj . 

The  algoritha  now  calculates  a  cross -correlation  function  for  the 

an  } 

M  data  points  using  Eq.  (2.9).  The  new  R.n(T.)  is  now  searched  for  a 

A®  j 

aaxiaua  value  using  Eq.  (2.10),  which  has  a  corresponding  unique  transit 
tine  T...  This  is  coapared  with  the  t  that  was  calculated  by  the 
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purtttioninc  routia*  using  ths  followint  squstioii: 

|t^  -  T^l  <  e^.  At  <  <  2At  .  (2.16) 

Th«  coapsrison  in  Eq.  (2.16)  chscks  to  see  if,  by  shifting  P  by  on*  At, 
tho  now  cnlculntod  is  still  within  tho  partition's  ability  to  resolva 
it.  If  the  result  of  the  eonpariscm  shows  that  the  naw  is  outsida  tho 
partition's  resolution,  than  a  new  partition  will  be  necessary  and  the 
algoritha  stops  averaging  and  proceeds  to  the  next  step,  data  shifting 
(see  Fig.  2.1).  If  Che  comparison  given  in  Eq.  (2.«6)  shows  that  the  new 
is  within  tc^  seconds  of  the  algorltha  adds  the  calculated 
K,b(t)  to  a  running  sun  and  proceeds  to  shift  P(A)  and  P(B)  by  one  wore 

AB  j 

At.  This  procedure  continues  until  at  the  X  ♦  1  shift  of  P(A)  and  P(B), 

Che  criterion  of  Eq.  (2.16)  la  not  net.  At  this  point  the  algorlthn 

calculates  an  average  cross-correlation  function  for  the  K  shifts 

of  the  partition.  This  function  is  then  searched  using  Eq.  (2.10)  for  a 

naxinun  value  with  its  corresponding  T^,  the  average  nininua  transit 

tine  for  the  tine  interval  0  ^  t  ^  (H  ♦  K)At,  where  M  eqiuls  the  site  of 

the  aininun  partition  P  and  K  equals  the  nuaber  of  At  shifts  nade  by  the 

partition.  The  algoritha  calls  T^,  the  new  aean  estiaate  of  the  ainiaua 

delay  tiae  ^^(ti),  and  places  this  estiaate  in  tine  at  the  aean  of  the 

M  ♦  K 

tine  interval  covered,  t)  >  — ^ — At.  The  algoritha  usw  calculates  V(ti) 
using 

V{ti)  -  0/7^g(ti)  ,  (2.17) 

where  0  is  the  distance  between  sensors  A  and  B;  V(ti)  is  also  place'i  in 
tiae  at  t|  seconds. 

Now  that  a  aean  velocity  V(ti)  has  been  estiaated  for  the  tine 
interval  0  <  t  <  (N  ♦  K)At,  the  algoritha  is  ready  to  aake  new  velocity 


I  MttMtM.  To  Oo  tkis,  Mch  iMtor't  block  1024  data  poiatt  la  ahlftod 

(N  ♦  K)/2  poiata*  wbich  corraapoada  to  a  tiaa  ahift  of  t**  ^  aacoada. 
Tha  data  bleeka  for  A  aad  •  aaaaora  aow  appaar  aa 

A  •  ♦  1*  *  *  **  *a^  C2»10) 

and 

•  •  {b^i  >»,  ♦  i*  •  •  •• 

Hhara  r  -  (N  «  K)/2  aad  S  •  1024  -  (M  «  K)/2.  Nota  that  thia  ahift  of 
tha  data  blocka  rataiaa  tha  laat  half  of  tha  old  data  uacd  to  aatinaca 
V(ti).  Thia  ia  dona  ao  that  no  poaaibla  now  croaa-corralation  batwaan 
tha  old  and  tha  naw  data  will  ba  loat  aa  a  raault  of  tha  shift.  Tha  algo- 
rithtt  now  chocks  to  saa  if  anough  data  points  (L)  ranain  in  tha  data  blocks 
to  continue  procassing .  It  doas  this  by  conparing  tha  maabar  of  data  points 
needed  to  resolve  V(t|)  with  the  nuaber  ranaining  after  the  shifts  using  tha 

following  expression: 

I 

L  •  s  —  r  >  47^(tj)/At  .  (2.20) 

If  tht  comparison  in  Eq.  (2.20)  shows  that  enough  data  reaain  in  A 
and  i  data  blocks,  the  algorithn  returns  to  tha  partitioning  routine  and 
starts  a  new  iteration  to  resolve  V(t2).  Tha  algorithn  will  continue  as 
above  until  the  criterion  in  Eq.  (2.20)  is  not  net,  at  which  tine  the 
algorithn  will  choose  a  naw  sampling  rate  based  on  the  last  velocity  esti- 
nate  y(t^)  and  proceed  to  sanple  the  A  and  0  sensors  at  the  new  sanpling 
race  foming  a  new  set  of  1024  data  points.  The  algorithn  will  then  con¬ 
tinue  and  repeat  Che  sane  Inner  Iterations  outlined  previously  to  resolve 
new  velocity  estlnates. 
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2.S  CoBClmti 


Aa  mttmft  Iws  b««i  madt  to  Aovolep  a  caafutor  alforitka  batod  oa 
CTOss-cemlatioa  uialytis  tachal^t  that  aill  track  a  aoastatioaaiy 
fluid  flow  volocity  aad  glvo  alaiaally  hiasod,  aaaiaally  prociao  otti« 
aatoa  of  tho  volocity  as  a  fuactioa  of  tiao.  Tho  algeritha  is  basod  oa 
tho  idoa  that  tho  total  aoasuroaMt  tiao  T  caa  bo  so  adjusted  as  to 
rojoct  all  but  tho  aiaiwa  flow  traasit  tiao  betwooa  tho  two  asasuring 
dovicos.  This  proeoduro  is  soaswhat  aaalogous  to  tho  uso  of  spoctral 
filtoriag  of  uawoatod  signals  ia  tho  ftoquoacy  doaaia. 

The  algoritha  is  curroatly  being  tasted  using  eoaputor-geaeratod 
volocity-transieat  data  sets.  Preliaiaary  results  show  good  agroeaoat 
between  the  known  velocities  aad  the  estiaotes  by  the  algoritha. 


Keference 


1.  J.  S.  Bcndat  and  A.  C.  Plersol,  Random  Cota:  Amlytia  and  Proetdureat 
p.  30,  Incerscicnce,  1971. 
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PHOTO  0091  77 


Air-water  downflow  in  the  Coring  bubble  regime 
intersects  a  6.4-mm-diam  (1/4-in.)  rod  transverse 
to  flow  in  a  7.6-cm-IO  pipe 
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PHOTO  0092  -77 


ir*J  •  !•■**  ilOr*!  »l«r'2  Si 

fcOIMUTV 

IWtt-phAM  flttl4  memmtvm  aMSurafrcncs  uslnt  •  drag  screen 
la  tlie  air-uacer  loop.  Curves  tLreush  daca  were  calculated  froa  • 
valacley  aodela  for  separated  flew  using  results  fron  a  turbine  aeter 
dewDstreaa  of  the  drag  screen  (solid  line«  Aye;  dasbsd  line.  Kouhani; 

4-20  flow  disperser  upstreaa  of  dreg  screen). 


Drag  screen  tested  in  BDHT  air-water  facility  shown 
in  coaparison  to  pipe  cross  section  (S-aesh.  0.64  m  dian  wire  screen). 


1 

! 


1.1-24 


MMOtC'Jin  ASSUHFTIM 


K  *  Vt 

•  pV*A|. 

SEPARATED  aW  HOOEL 

*  " 

AYA  HODEL 
Turbiif  ItetT 

CiaPj(Vj  -  V^)*  -  Cid  -  a)p^(V^  -  V^)* 
PfQ  Disk 

— s  -  .2 

^  ♦  C%(l  — 

Ontltoweter 


Ph«M  Continuity 

\ 

"•  "  ^9^9 


1.1-29 


USNKC  Sponsored  Instruacntatinn  Roscarrh 

at 

Rensselaer  Poly technic  Institute  (RPIJ 


R.  T.  Lahey*  Jr. 
Principal  Investigator 


The  current  sponsored  instruaeiitation  research  at  RPI  can  be 
conveniently  divided  into  two  categories. 

(1)  Instrunentation  to  be  used  in  pressure  suppression 
containnent  systea  experiments  to  determine  the  steam  content  in 
a  steam/air  flow. 

(2)  Two-phase  flow  instnimcntation  which  can  he  used  to 
determine  the  volumetric  vapor  (void)  fraction. 

This  paper  is  a  status  report  of  the  progress  during  the  first 
two  quarters  of  a  twelve  quarter  program.  First  the  facilities 
which  have  been  constructed  in  which  to  develop  the  instrumentation 
will  be  described; chen  theory  and/or  instrumentation  will  be  dis¬ 
cussed.  Exploratory  "bench  top"  test  results  for  various  novel 
two- phase  flow  measurement  devices  will  also  he  presented. 

The  first  figure  (Figure  1)  is  a  schematic  of  the  stcam/air 
test  loop.  On  the  left-hand  side  of  the  figure*  one  can  note  that 
the  compressed  air  is  heated  in  a  feedback  controlled  (for  outlet 
temperature)  air  preheater.  This  preheating  is  done  to  avoid  ex¬ 
cessive  condensation  of  the  steam*  which  mixes  with  the  air  in  a 
mixing  tee  at  the  bottom  of  the  test  section.  The  steam  is  pro¬ 
duced  in  a  pressure  vessel  and  is  piped  into  the  mixing  tee.  The 
steam/air  mixture  then  flows  up  the  3/4- inch  square  aluminum  test 
section.  Quartz  windows  are  shown  in  the  top  of  the  test  section. 
The  infrared  (IR)  beams  ore  passed  through  these  windows  to  deter¬ 
mine  the  steam  content.  Figure  2  shows  the  actual  apparatus  (with 
the  insulation  removeci).  The  air  is  piped  in  on  the  right  th;oiip.h  .i 
pressure  regulator,  oiificc  flow  meter  and  t!ie  air  preheater.  The 
3  kw  steam  generator  is  shown  and  the  vertical  insulated  tost  sec- 

Pmidiiit  mi  lilak 


I.M 


tioa  is  sho^  in  pisce.  Not*  ths  quarts  wiiutows  near  tha  top 
of  the  tost  section  aa4  the  plat  fora  on  which  to  uount  the 
l.K.  instruaentation. 

The  principal  behind  the  operation  of  an  l.R.  device  is 
quite  straight-forward.  The  intensity  varies  as 

where, 

>  Monochrooatic  intensity  of  the  beau  (nhotons/m^-ser j 
■  Monochroaatic  intensity  if  no  absorbing  neJiun 
(i.e.  steau)  is  present 
•c^  ■  Monochroust ic  uacroscopic  cross  section 
X  •  Beau  length  through  the  absorbing  uediuB 
(i.e.  stcan). 

It  is  .well  known  that  steau  absorbs  soae  wave  lengths  (x)  much 
aore  than  others.  The  uethnJ  used  is  to  pass  two  wave  lenntlis 
(X|  and  X2)  through  the  stean/air  mixture.  Ne  chose  X2  such  that 
it  is  strongly  absorbed  by  steau  and  Xj^  such  that  it  is  not. 

Figures  3a  and  3a  present  the  theory  behind  the  technique. 

The  result  is  that  we  can  directly  dereruinc  the  steau  flow  frac* 
tian  (y)  using  a  two-beau  l.R.  device.  The  work  during  the  next 
two  quarters  will  be  directed  toward  constructing  and  calibrating 
the  l.R.  device. 

l-et  us  now  turn  our  attention  to  two-phase  flow  ueasureuent 
devices.  As  in  the  previous  tTiscussion,  let  us  first  consider 
the  facilities  which  have  been  constructed  in  which  to  develop 
such  devices.  Figure  I  shows  a  photo  of  the  gPl  squjII  air/watcr 
loop.  This  loop  is  used  forcomepru.il  ’lOp”  testinj*.  of 

various  instriincnts.  Figure  S  is  ;i  schcnatic  of  the  RPT  sicdium 
size  air/watcr  loop.  It  can  he  seen  that  in  this  loop 
Wi!  onl)'  have  air  flow  (i.e.  no  axial  water  flow).  'Ihe  ai>' 
is  intioduccd  into  the  hotton  of  the  Jowci  ptennni  thrniE;;li  h 
sparger.  \  dispersion  plate  iS  installed  which  distrihutes  the 
vapor  and  establ  i.shes  the  bubble  sir.c. 
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A  Dresser  coupling  is  installed  so  that  various  glass  or  pleai- 
glass  test  sections  can  be  easily  installed.  Figure  6  shows 
this  loop  in  operation. 

Various  two-phase  flow  regiaes  can  be  estrblished  in  the 
one-inch  I.O.  test  section.  Figures  7a-7d  show  exawples  of:  low 
quality  bubbly  flow,  high  quality  bid>bly  flow,  slug  flow,  and 
annular  flow,  respectively. 

The  aediun  size  air/water  loop  will  be  the  "work  horse"  in 
which  the  instruaentation  is  developed.  In  order  to  carry  out 
this  developaent  work  one  needs  to  accurately  know  what  the  void 
fraction  actually  is  so  that  the  instruaentation  under  develop- 
aent  can  be  properly  calibrated  and  interpreted.  The  standard 
we  have  chosen  is  a  dual  beaa  X-ray  device  shown  scheaatically  in 
i^igure  8.  The  basic  theory  behind  this  technique  has  been 
successfully  used  by  others  (Jones,  Saiih,  Schrock,  Staub,  etc.) 
however,  the  techniques  previously  used  to  process  the  data  were 
cumbersome.  The  theoretical  basis  for  the  experiment:!  1  technique 
we  plan  to  use  is  presented  in  Figures  9a- 9b. 

In  theory,  if  we  chose  the  0.  C.  level  (1>)  of  the  reference 
beam  correctly  we  can  get  rid  of  the  tmdesirable  60  Hi  ripple  in 
both  beans.  It  can  be  seen  that  the  resulting  expression  for 
instantaneous  chordal  average  void  fraction,  a(l),  is  given  by 


where  the  R  are  log  ratios  of  the  test  section  and  reference  beans. 
This  simple  linear  fora  suggests  an  easily  calibrated  circuit 
design. 

Figure  10  is  a  schematic  of  the  electronics  systems.  It  can 
been  seen  that  it  follows  the  logic  flow  presented  in  Figures 
9a 'Ob.  It  can  also  be  seen  that  not  only  is  the  instantaneous 
chordu I  average  void,  a(t),  obtnine«l,  bur ,  by  passing  this 
stochastic  signal  through  a  time-integr  iting  digital  volt  netcr 
(OWt)  the  tine  averaged  chordal  average  void  fraction,  a,  is  also 
obtained. 

The  merit  of  such  a  system  is  expected  to  he  the  ease  with 
which  it  can  be  calibrated.  With  a  test  section  full  of  water. 
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the  R|^  potMtiuaeter  can  be  adjusted  until  a  ■  0.  Then  with  only 
air  in  the  test  section,  the  potcntioeeter  can  be  adjusted  to 
yield  a  •  1.0.  After  that  one  can  obtain  a  real^tlsM  readins  of 
'  on  the  DVM. 

This  calibration  oust  be  done  at  each  transverse  position, 
but  it  is  expected  that  this  electronic  design  should  greatly 
speed-up  data  acquisition. 

In  addition  to  a  we  arc  also  interested  in  looking  at  the 
statistical  inforaation  in  the  void  signal  to  see  if  we  can 
develop  a  criterion  for  interpreting  the  signature  of  various 
flow  regis«s.  In  order  to  accoaplish  this  the  stochastic  a(t) 
signal  is  passed  through  a  low  pass  filter  having  a  100  Hz  break 
frequency.  The  purpose  of  this  filter  is  to  reject  unwanted 
high  frequency  noise  and  to  avoid  alaising  errors  when  the  signal 
is  digitized.  ^ 

It  is  planned  that  the  probability  density  function  (PDF) 
will  be  acciiaulatod  in  a  wultichannel  analyzer  (HCA) .  The  PDF 
will  be  taken,  via  paper  tape  output,  to  a  PlJP-15  computer  for 
further  analysis  (e.g.  determination  of  the  moments^  In  addition 
to  the  PDF,  the  power  spectral  density  function  (Pfil>)  will  also 
be  constructed.  This  will  he  accomplished  by  digitizing  the 
signal  with  an  (ADC)  analog-to-digital  convertor,  and  Fourier 
transforming  the  signal  with  a  FFT  algorithm  on  the  PPP-lS  data 
acquisition  computer. 

N-'>w  that  we  have  discussed  the  facilities  in  which  to  ^ 
develop  multiphase  instrumentation,  and  the  Xray  standard  against 
which  the  measurements  will  be  calibrated,  wo  are  ready  to  dis¬ 
cuss  various  two-phase  instrumentation  under  development  at  RPl. 

Figure  11  is  a  schematic  of  a  visible  red  fV'.ll.)  test  appara¬ 
tus.  This  device  was  pulsed  at  S3  Kiiz  and  passed  through  the 
transparent  test  section  in  the  small  air/water  loop.  The  top 
view  of  the  oscilloscope  trace  shown  in  Figure  12  shows  the  basic 
signal  (on  an  expanded  time  scale)  when  only  water  is  in  the 
test  section.  The  middle  and  bottom  figures  show  how  the  basic 
signal  was  modulated  by  low  quality  bubbly  flow  and  high  quality 
bubbles  flow  rc  i^ctively.  It  can  bo  seen  tiiat  the  signal  is 
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'  atteaiMtvd  by  bubbles  and  is  very  sensitive  to  voids. 

Figure  13  is  a  schenatic  of  an  experiaent  to  detemine  wby 
bubbles  cause  l.R.  bean  attendation  (i.e.  how  much  intenaitv  loss 
is  due  to  refraction  and  hou  nuch  by  reflection  at  the  interface. 
Four  l.R.  transaitter'* sensor  units  are  shoun.  One  test  was  run 
in  which  unit  11  was  illuninating  unit  #2.  In  this  test  we  were 
looking  for  l.R.  bean  attenuation  and  scattering  (i.e.  refraction 
and  reflection)  on  the  various  units  shown. 

Figure  14  shows  the  actual  unit  in  place.  Near  the  top  of 
the  test  section  can  also  be  seen  a  local  probe  which  we  will 
discuss  shortly. 

Figure  IS  shows  the  oscilloscope  trace  of  the  various  l.R. 
units.  For  low  quality  bubbly  flow  we  had  fine  (‘«>l/32  inch  in 
diaaeter)  air  bubbles*  and  for  higher  quality  bubbly  flow*  we  had 
■ediuB  size  (‘«>3/16  inch  in  diaaeter)  bubbles.  Note  in  the  two 
top  figures  in  figure  IS  that  the  presence  of  bubbles  attenuate 
the  l.R.  bean*  thus*  since  the  exponential  attenuation  gets  less 
(i.e.  air  has  a  auch  smaller  aacroscopic  cross  section  than  water)* 
we  are  seeing  the  effect  of  scattering.  Obviously  the  sensor  in 
I  unit  t4  has  l.R.  scattered  in  but  this  does  not  appear  to  be 

correlated  with  the  scattering  which  occurs  on  bubbles  in  the 
aain  (1*2)  beaa. 

In  the  lower  left  picture  of  figure  15*  the  1*1  back  scatter* 
ing  signal  is  shown.  Surnrisingiy  sensor  fl  is  noraally  illu- 
ainated  and  any  scattering  apparently  causes  a  loss  in  illuaination. 
The  reason  for  this  is  not  clear  at  the  present  time.  The  final 
picture  on  Figure  14  shows  that  there  is  very  little  scattering 
into  sensor  #3. 

It  is  clear  that  wc  do  not  have  isotropic  jicattering  by  the 
bubble.  It  is  also  clear  that  further*  single  bubble  experiaents 
arc  needed  to  quantify  the  scattering  acclianisns.  One  likely  test 
apparatus  for  further  testing  is  show.,  in  Figure  16. 

Several  local  probes  arc  under  development.  Figure  17  shows 
a  new  frequency  nodulated  inpedance  (Z)  probe.  The  probe  itself 
is  constructed  from  gold  alloy  wires  inserted  in  a  ceramic  tube. 

Gold  is  used  to  avoid  possible  corrosion  probleos.  The  probe 
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forM  ui  activ*  alaMnt  !■  the  eecillator  circuit.  Vheu  the 
probe  is  ie  water  a  lower  frequency  pulse  is  ohsenred  than  when 
it  is  in  the  vapor  phase  (in  the  case  shom  625  Hz  for  water  and 
769  Hz  for  air.  althoufh  these  frequencies  can  be  set  as  high  as 
you  like).  These  pulses  are  picked  up  by  a  phase* locked-loop  (Pit.) 
which  sends  out  a  signal  of  S  volts  for  the  air  frequency  and  0 
volts  for  water.  This  technique  is  known  as  frequency*shi ft*keying 
(FSK)  and  is  known  to  be  rather  insensitive  to  ninor  pulse  fre¬ 
quency  changes.  This  inplies  that  the  probe  nay  not  be  sen¬ 
sitive  to  water  purity  and  tenperature  changes,  a  highly  desirable 
feature  for  a  local  Z  probe.  Finally,  as  shown  in  Figure  17.  the 
count  statistics  can  be  quite  accurate  since  the  lobes  of  the 
square  wave  coning  out  of  the  PLL  can  be  filled  with  high  frequency 
pulses  using  a  pulse  generator.  The  signal  is  then  passed 
through  a  single  channel  analyzer  which  dutermines  the  number  of 
coimts  in  air 

The  final  probe  to  be  reported  on  is  shown  in  Figure  18. 
i'he  probe  itself  is  the  sane  as  that  used  in  the  Z  probe  tests. 

The  difference  is  how  the  probe  is  excited,  in  this  case  the 
probe  is  driven  by  a  radio  frequency  (RF)  generator.  This  genera¬ 
tor  is  coupled  to  the  probe  through  a  tuned  passive  circuit.  The 
probe  is  a  capacitive  element  in  this  circuit.  In  the  case  shown, 
the  circuit  has  been  tuned  to  17  MHz.  although  this  is  quite 
a.-bitrary.  If  the  RF  generator  is  set  at  17  tSIz  then  we  have  no 
indication  of  phase.  Instead,  if  we  set  the  RF  generator 
below  17  Miz  (e.g.  16  MIz) ,  we  get  pulses  when  we  have  bubbles 
hitting  the  probe.  In  contrast,  if  we  set  the  RF  generator  above 
17  MHz  (e.g.  18  MIz)  the  vapor  phase  causes  a  reduction  in  carrier 
amplitude.  Actually  scope  readings  are  shown  in  Figure  19. 

The  RF  probe  is  also  expected  to  be  quite  insensitive  to 
water  purity  and  temperature  and.  due  to  the  high  frequency, 
should  have  quite  accurate  count  statistics. 
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FIQME  3« 

ULTIlkMINlNC  STn/VM/AIR  RATIO  BY  I.R.  TI-CHWIQUI’S 


(For  »,)  I,  -  e'-l* 

II or  A,)  1,-1  e'‘'2* 

•  A  Oy 


(1) 

(2) 


TakiiiK  log  of  Equations  (I)  and  (2),  and  subtracting 
tndj/ip  -  ‘n(Io/*o,J  *  <"2  '  “l^** 

where: 


C3) 


Ai  i 


i  -  i  (when  no  stcaa) 

A-. 


(4) 


(51 


rquatiens  (3)  and  (4),  with  the  assumption  Oj*0  (i.c.,  Hj-0), 
yield: 


^1?  "2^  "2  •  E 


or. 


**2  “  fT^JU  ‘«nj/i2) 


O  2 

Thus, 

•’steaitt  "  4hM,/I,J 


(h) 


(7| 
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Recalling  Dalton's  Law: 

^  "  Pair  *  Pstoaa 

Pstaao  *  ’'stean  *^steaa  ^sat 

Conbining  (R)  aiwl  (9): 

Pair  ■  P  ■  *itean  *^stean  ^aat 

But, 

Pair  *  ‘‘air  *air  ’’sat 
Thus,  conbining  (10)  and  (11): 

‘‘air  *  Rjijy  L  P^^sat  ’  ^stean  *^stean  ^ 
The  steaa/air  fraction  (y)  is: 

y  R  _ *^steaa  stean 

^■^air^air  *  **$rean^*steani 

assuning  •  U,.^. 

y  ■  **stcan 

I'^air  *  ^steanl 


(•) 

(9) 

(10) 

(It) 

(12) 

(IS) 

(Hj 


Hquations  (7),  (12)  and  (14)  yield  the  required  result. 
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FIGURE  7b 

HIGH  QUALITY  BUBBLY  FLOW 
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SCHEMATIC  OF  X-RAY  EQUIPMENT 
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DERIVATION  OF 
VOID  FRACTION  EQUATION 


IN  GENERAL; 
I(t)  • 


FOR  THE  TEST  SECTION  BEAM; 

(1)  A*(t)  •  A(t)  L  I  ♦  •  sin«*  1 

For  the  Reference  Beaa  (with  D.  C.  Bias); 

B*(t)  ■  B  (_  I  ♦  b  sinut  )  ♦  D 

or, 

(2)  B  (t)  "fB^D]  !!♦  sinut  | 

taking  the  log  of  (1)  and  (2),  and  subtracting; 


(3)  (n 

B  (t) 


tn 

fB^D] 


f  1  a  sinut  I 

^  *  *  tIW  3 


If  we  set  the  D.  C.  bias  at, 

D  -  B(b/a-\ 

the  last  term  in  uisappears,  hence  (3)  can  be  written 


in 


1  (t) 


■  in 


[  J 

I  B*D  j 


in 


1  *  ux(t) 
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Th«i», 


FIQIK  9b 


14)  x(t)  -  i 


tn 


B  (t) 


The  chordal  average  void  fraction  is  given  by 

(5)  o(t)  -  1  - 

where  C  is  chord  length  through  fluid. 

•  For  an  capty  pipe,  a-1,  and  x*0, 
thus  (4)  yields: 


A  1 

91 

(6)  tn 

■  tn 

.8  (t) 

•  For  a  pipe  full  of  water,  a*>0 
Equation  (5)  yields: 


(7) 

(8) 

(9) 


x-C 

and  Equations  (4)  and  (6)  yield: 


where. 


Combining  Equations  (4),  (5),  (6)  and  (8); 


(10) 
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ABSTRACT 


A  program  was  begun  on  July  1,  1977,  the  first  part  of  which  haJ  as  Its  main 
objective  the  development  of  analytical  expressions  for  the  non-equilibrium  rates 
of  vapor  generation  under  various  circumstances  of  Interest  In  accident  analysis. 

In  support  of  this  development  effort,  an  experimental  program  Is  being  undertaken 
to  measure  the  actual  vaporization  rates  In  flashing  flows.  Instrumentation  Is 
being  constructed  and  developed  In  support  of  this  effort  to  obtain  both  global  and 
local  measurements  necessary  to  provide  the  desired  Information.  This  Instrumentation 
Includes  a  global  densitometer  using  radiation  attenuation  techniques  and  miniature 
optical  probes  for  obtaining  local  Information  regarding  vapor  residence  time  fractions 
and  Interfaclal  passage  velocities.  This  presentation  sunmarlzes  the  work  to  date 
and  the  current  status  of  this  portion  of  the  overall  program. 


OVERVIEW  OF  PROGRAM 


The  overall  program  consists  of  two  parts  as  shown  in  Figure  1: 

Development  of  Constitutive  Relations;  Development  of  Instrumentation  Techniques. 

Task  1.1  Involves  the  development  of  calculatlonal  procedures  applicable  to 
advanced  confirmatory  codes  which  will  allow  accurate  predictions  of  the  nor.'X|u111br1ifln 
rates  of  vapor  generation  to  be  made  under  such  conditions  as  post-dryout  iiiit 
transfer,  sub-cooled  boiling,  downcomer  condensation,  and  decompressive  flaiMng. 
Experiments  are  being  planned  and  constructed  to  provide  a  broad  base  of  direct 
measurements  of  vapor  generation  In  flashing  flows.  This  effort  Is  identified  In 
Tasks  1.2  and  1.3.  In  support  of  this  experimental  program,  certain  two-phase 
Instrumenration  methods  are  being  developed.  Development  of  these  techniques 
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constitutes  Pert  2  of  the  overall  program  and  It  Includes  both  local  and  global 
Instrumentation.  The  global  Instrumentation  of  direct  application  to  the  current 
experiments  Includes  the  global  densitometer  which  Is  based  on  gamma-attenuatlon 
techniques,  and  Involves  no  major  extension  of  previously  existing  technology. 

The  local  Instrumentation  required  for  measurement  of  local  void  time  fractions 
and  Interfaclal  passage  velocities.  Involves  extension  of  existing  methodology 
Into  conditions  beyond  those  of  current  application.  These  conditions  Include 
both  the  application  to  higher  te<.>perature  flows  and  higher  speed  flows  than 
have  been  encountered  In  previous  Investigations. 

In  addition  to  the  Instrumentation  development  work  of  direct  support  to  the 
current  experiments,  loi.ger  term  goals  of  analytical  and  experimental  evaluation 
of  global  Instrument  techniques  applicable  to  LOFT/seml-scale  type  Instruments  are 
Included  In  this  program.  Work  In  this  vein,  however.  Is  not  expected  to  begin 
until  after  fiscal  1977. 


IDEOLOGY 


The  experimental  measurements  are  designed  ro  provide  Information  for 
direct  determination  of  the  volumetric  rate  o*  vapor  generation,  r  ,  In 
flashing  flows  in  ducts.  In  order  to  do  this  u  i  Is  made  of  the  continunity 
equation  for  the  vapor  phase  given  by 


?t  3Z 


(la) 


or  alternately,  the  Ci'n  'uitv  equation  for  the  liquid  phase  given  by 


3<(1  -  a)p  >  3 

- —  +  — <M  -  a)p,>/ >  ’  r,  =>  -r 


(lb) 


3t 


3Z 


Under  steay-state  conditions,  it  Is  seen  that  the  determination  of  the  actual 
volumetric  rate  uf  vapcr  producti»«.  Involves  both  the  measurement  of  the  local 
phase  volume  fraction  as  well  as  the  local  phase  velocity,  which  is  then  integrated 
over  the  cross-sectional  area  to  obtain  cross-sectional  averaged  values.  These 
va’u’s  must  then  be  obtained  at  different  axial  locations  so  that  effective 
detfcri.dnation  of  the  axial  ^-idient  of  the  phase  mass  flux  can  be  determined.  In 
actual  practice,  especially  at  low  reduced  pressures,  interfacial  velocities  will 
most  nearly  follow  the  liquid  velocity  as  a  result  of  the  interfacial  mass 
balances.  Tha*-  Is,  the  vapcr  produced  must  move  away  from  an  interface  with  a 
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for  the  case  of  a  constant  area  duct.  The  errors  Involved  In  determination  of 
r  are  thus  associated  with  errors  In  determination  of  the  mass  velocity,  quality, 
and  ar'<al  positions  at  which  the  quality  Is  determined.  We  must  In  effect  measure 
the  axial  quality  gradient.  The  relative  standard  deviation  In  the  vapor  generation 
rate  Is  thus  given  by 


where  If  both  the  mass  velocity  and  the  actual  positions  are  accurately  known, 
the  largest  source  of  error  occurs  In  the  measurement  of  the  quality  itself.  In 
terms  of  the  volume  fractions  and  velocities,  the  quality  may  be  given  as 
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X  ■ 


1 


(5) 


1 


♦ 


such  that,  neglecting  property  errors,  the  relative  standard  deviation  In  the 
quality  Is  given  by 


V2 


(6) 


The  results  of  Equation  (6)  are  plotted  In  Figure  2.  It  Is  seen  that  In  the 
range  of  atmospheric  pressure  to  450  psi  that  even  for  quite  large  values  of 
relative  velocity,  v  ,  the  error  In  quality  determination  can  be  maintained  at 
a  relatively  low  level  providing  the  relative  standard  deviations  In  the 
measured  void  fraction  and  liquid  velocities  are  kept  low. 


METHODOLOGY 


The  methods  to  be  used  In  measuring  the  non-equilibrium  vapor  generation 
rates  In  flashing  flows  are  based  directly  on  the  previously  described  Ideology. 

That  is.  we  intend  to  provide  a  method  for  establishing  a  steady  circulation  of  water 
at  a  well  controlled  level  of  temperature  and  flow,  and  utilize  a  test  apparatus 
Inserted  In  this  circulation  loop  which  allow  determination  of  the  flashing 
conditions  under  the  following  requirements: 

a)  Minimum  power  requirements; 

b)  local  Information  may  be  obtained; 

c)  global  reference  data  may  be  obtained; 

d)  minimum  effects  due  to  upstream  Instrumentation  disturbances. 

In  addition  to  the  loop  and  the  test  section,  suitable  sensing  devices  must  be 
provided  to  yield  the  appropriate  Information  regarding  phase  volume  fraction  and 
velocities,  and  suitable  methods  must  be  developed  for  handling  and  analysis 
of  the  signals  derived  from  the  Instrumentation.  Each  of  these  Items  shall  be 
discussed  In  turn  below. 


Test  Facility 


A  scheratic  diagram  of  the  test  facility  to  be  used  In  these  experiments 
Is  shown  In  Figure  3.  This  test  facility  Is  patterned  after  the  HOBY  DICK  loop 
currently  *-  use  at  CENG.  The  facility  Is  being  designed  and  constructed  for 
maintalnli  .le  appropriate  circulatory  conditions  with  static  temperature  and 
pressure  ratings  of  150  psi  at  350  F. 

Status  of  the  test  facility  to  date  Is  as  follows 

a)  The  major  test  facility  modules  are  on  hand  Including  pumping  modules, 
heat  exchangers,  water  purification  system,  water  filtration  system, 
steam  drum,  control  panels,  and  miscellaneous  valves  and  fittings; 

|>)  The  thermodynamic  design  Is  approximately  80  to  90X  complete; 

c)  the  component  Interconnection  design  Is  approximately  50t  complete. 


Signal  Analysis  System 


The  signal  analysis  system  Is  shown  schematically  In  Figure  4.  This  system 
I  Is  a  mini-computer  based  system  which  will  acquire  slow  speed  and  Intermediate 
I  speed  Information  directly,  and  will  acquire  high  speed  Information  through  play- 

{  back  of  previously  recorded  analog  tape  data  at  a  reduced  speed.  In  the  design  of 

'  this  system,  emphasis  has  been  placed  on  capability  for  future  expansion, 

rapid  turn-around  and  quick  visual  display  of  acquired  Information,  and  Interactive 
!  flexablllty  with  users.  The  system  will  operate  In  a  real  time,  multi¬ 
user,  multi-program  environment,  which  will  allow  for  simultaneous  development 
of  different  capabilities  by  more  than  one  user.  It  Is  Intended  that  the  calculatlonal 
and  operational  capabilities  of  the  HP-54518  Fourier  analyzer  will  be  emulated. 

The  status  of  this  system  Is  as  follows: 

a)  Honeywell  9600,  4MHZ  tape  recorder  ha*  been  received,  checked  out,  and  Is 
operational ; 

b)  The  digital  sys.em  control  and  display  modules  have  been  received,  Intergrated, 
and  are  operational; 

c)  Driver  design  and  development  Is  In  progress; 

d)  The  multi-user,  real  time  operating  system  Is  operational; 

e)  The  specifications  for  potential  vendors  of  the  multiplexed  analog  digital 
conversion  system  are  under  study. 
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Scheutic  diagram  of  the  test  section  are  shoM  In  Figures  5  and  6.  It 
Is  seen  that  this  test  section  Involves  a  converg1ng*d1 verging  region  through 
which  a  probe  rake  having  a  number  of  local  probes  can  traverse  the  test  section  In 
both  radial  and  axial  directions.  While  relatively  simp'e  In  concept,  the  test 
section  Is  quite  complex  In  execution  and  major  development  efforts  are  associated 
with  the  probe  rake  Itself,  the  pressurized  linear  flow  seal,  and  the  positioning 
slide  mechanism,  as  well  as  methods  for  mechanical  support  yielding  minimum 
obstruction  to  the  gamma*dens1 tome ter. 


Global  Densitometer 


A  bidirectional,  scanning,  gamma -densitometer  Is  being  designed  to  provide 
reference  Information  for  the  local  void  fraction  measurements  which  are  known 
to  be  Inaccurate  by  as  much  as  10  to  15X.  The  local  measurements  will  be  corrected 
by  comparison  with  tiie  global  measurements.  The  designs  specifications  for  the 
gamma  densitometer  are  as  follows: 

a)  multi -beam  capability  with  radial  coverage  for  pipe  sizes  up  to 
2  Inches 

b)  transient  accuracy  shall  be  wltnln  St  for  any  1.0  m  sec  time  slice 

c)  output  shall  be  linear  with  fluid  density  within  *  21 

d)  The  maximum  steady-state  error  shall  be  2t 

e)  axial  scanning  capability  will  be  24  Inches 

f)  transverse  scanning  capability  will  be  over  6  Inches 

g)  a  calibration  mode  shall  be  Included  in  the  device 

h)  the  maximum  source— detector  distance  shall  be  4  Inches 

Design  studies  of  the  system  have  been  completed  to  date.  For  the  test 
section  geometries  shown  In  Figures  5  and  6,  the  energies  have  been  optimized 
to  the  30-35  Kev  range.  Since  no  suitable  nuclear  sources  exist  In  this  exact 
range,  an  Iodine  or  thulium  source  shall  probably  be  used.  The  detection 
system  shall  be  based  on  cadmium-tellurlde  detectors  used  in  the  pulsed  mode  with 
pulse-to-volt  conversion  and  logarithmic  linearization  of  the  signal. 

Local  Probes 


The  local  detection  of  phase  volume  fraction  and  velocity  shall  be  measured 
using  dual  optical  probes  in  a  set  up  similar  to  that  shown  in  Figure  7.  Because 
It  is  expected  that  under  intermediate  and  high  speed,  two-phase  flow  conditions 
the  phase  structure  will  be  quite  amorphous,  a  high  interfacial  passage  frequency 
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Is  llktiy.  If*  for  Instance,  a  pair  of  liquid  vapor  Interfaces  separated  by 
a  distance  of  0.5  m  passes  by  the  sensor  at  a  velocity  n/s*  tne  pulse- 
width  Is  10  vns.  By  examination  of  the  Fourier  represer*  n  of  a  square  wave 
havlngtdomlnant  frequencies  of  100  khz*  It  Is  seen  that  unoistorted  recording  of 
the  20  harmonic  Is  necessary  to  obtain  accurate  zero  and  unity  crossing  times 
within  approximately  3X.  Transient  measurement  and  recording  capabilities  flat  to 
2.0  Hhz  are  thus  seen  to  be  required.  Indicating  measurement  bandwidths  of  4.0  Hhz. 

As  a  result*  attention  must  be  paid  to  the  fine  details  of  Interaction  of  the 
probes  with  the  flow  structure  Itself.  The  schematic  design  shown  In  Figure  7 
attempts  to  minimize  these  Interactions. 

The  optical  sensor  design  currently  being  pursued  Is  based  on  the  Joining  of 
two  rectangular  optical  wave  guides  having  been  previously  beveled  at  the  end  to 
an  angle  of  45°  with  respect  to  the  axis  of  the  wave  guide.  The  resultant  90° 
angle  at  the  tip  of  the  probe  will  thus  have  the  capability  of  reflecting  light 
back  through  a  second*  seperate  wave  guide  from  the  source  wave  guide,  and  having 
a  very  narrow  zone  of  fluid  Interaction  with  a  high  degree  of  rejection  of  light 
sources  external  to  the  sensor  tip  Itself.  The  Ideal  signals  are  shown  In  the 
lower  half  of  Figure  7.  The  phase  volume  fraction  Is  obtained  by  a  simple  time 
average  of  the  Ideal  signal,  while  the  interfaclal  passage  velocity  Is  obtained  by 
determining  the  time  lag  for  the  Interface  to  pass  between  the  two  sensors  seperated 
by  a  known  distance.  In  practice,  since  the  leading  probe  may  also  act  as  a 
nucleatlon  source  of  suffflclent  Intensity  to  mask  the  Inoividual  transient  time 
Information,  It  may  prove  necessary  to  utilize  cross  correlation  between  the  two 
signals  to  determine  this  velocity  on  the  average.  This  would,  of  course.  Imply 
the  assumption  of  Independence  between  passage  time  and  Interfaclal  velocity,  or. 

In  other  words.  Independence  between  Interfaclal  velocity  and  void  structure.  For 
high  speed  flows,  this  should  not  be  a  bad  approximation.  At  this  time  vendors  have 
been  contacted  who  Indicate  they  can  supply  the  material.  Necessary  procedures  to 
effect  the  shaping  and  joining  of  the  Individual  fiber  should  not  pose  difficulties. 
It  Is  expected  that  quotations  for  the  specific  design  shall  be  obtained  by  the  end 
of  March. 

With  regard  to  the  local  Instrumentation,  there  are  several  areas  that  require 
development  to  one  extent  or  another.  These  are  shown  In  Table  1  along  with  possible 
solutions  being  persued  In  each  area.  As  far  as  the  mechanical  design  Is  concerned, 
alternate  designs  are  elso  being  examined  as  shown  ir  Figure  8.  The  emph*sis  in 
these  alternate  designs  Is  mainly  on  producing  a  single  fiber  detector  so  as  to 
minimize  the  internal  volume  requirements  for  the  probe  rake,  and  reduce  interaction 
with  the  flow  structure  to  a  minimum.  The  problem,  of  course.  In  sensors  of  such 
small  size.  Is  In  obtaining  a  suitable  bifurcation  method  that  will  allow  seperatlon 
of  source  and  detector  to  an  acceptable  degree.  There  Is  no  doubt  that,  should 
all  previous  methods  fall,  the  U-bend  sensor  developed  by  Delhaye  can  be  used.  - 

With  respect  to  the  other  development  areas,  the  work  effort  Is  In  various 
stages  of  progress.  Methods  of  accounting  for  electronic  signal  distortion  have 
been  examined,  and  signal  reconstruction  techniques  based  on  Fourier  inversion 
processes  have  been  tried  and  found  to  be  acceptable.  The  methodology,  shown 
schematically  in  Figure  9,  simply  Involves  the  reamplification  and  back  phase 
shifting  of  each  frequency  component  in  the  Fourier  spectrum  of  a  particular  signal 
according  to  the  previously  determined  transfer  function  of  the  measurement  system. 
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This  Is  tquillvaltnt  to  division  of  tht  signal  hy  tho  coaplax  transfer  function 
In  the  frequency  doMln.  Difficulties  Involved  In  this  Methodology  arrise  due 
to  aliasing,  frequency  broadening  due  to  the  sanpllng  freouency  not  being  an 
Intergral  Multiple  of  a  characteristic  frequency,  errors  due  to  truncation  of  the 
tines  series,  and  decreasing  signal  to  noise  ratio  as  the  frequency 
conponent  Increases. 

In  addition  to  sinple  electronic  signal  distortion  on  the  signal  neasurenent 
channel,  problems  associated  with  time  Incoherrancles  when  examining  correlations 
between  two  parallel  channels  have  posed  some  difficulties  In  the  past.  When  the 
auto  correlation  of  a  particular  signal  Is  examined.  It  Is  seen  to  have  a  maximum 
value  equal  to  the  mean  square  of  the  signal  at  zero  time  delay.  If  an  Identical 
signal  Is  recorded  on  two  parallel  channels,  then  the  cross  correlation  between 
two  channels  will  have  a  maximum  equal  to  the  mean  square  value  of  the  signal, 
(assuming  Identical  amplification  In  each  channel),  but  having  this  maximum 
occuring  at  a  value  for  time  delay  equillvalent  to  the  recording  Incoherrancy  between 
the  two  channels.  It  Is  expected,  then,  that  this  recording  Incoherrancy  can  be 
determined  within  the  required  accuracy  by  a  simple  cross  correlation  of  Identical 
signals  recorded  on  two  or  more  channels. 

With  regard  to  the  latter  two  development  areas,  hydrodynamic  signal  distortion, 
and  the  relation  between  Interface  and  phase  velocity.  It  is  anticipated  that  these 
problems  can  be  attacked  by  the  use  of  Independent  global  references.  That  Is, 
hydrodynamic  disortlons  lead  to  Inaccuracies  In  both  the  phase  volume  fractions, 
and  the  measurement  of  the  actual  velocities.  As  previously  Indicated,  the  ganna- 
densltometer  Is  expected  to  be  used  as  an  Independent  global  reference  for  the 
local  void  fraction  measurements.  In  addition  to  the  void  fraction  reference, 
there  are  other  references  which  have  accuracies  at  least  as  good  as  that  expected 
to  obtained  from  the  gamma  densitometer.  These  include  the  actual  mass  flow  in 
the  test  section,  and  the  total  energy  flow  through  the  test  section.  Since  the 
mass  velocity  and  energy  flow  are  both  Integral  values  of  the  localized  phase  values, 
it  is  seen  that  there  are  two  additional  Independent  references  that  can  be  used  for 
calibration  and  correction  procedures  should  the  need  be  determined. 


SUMMARY 


This  presentation  has  summarized  In  general  the  overall  purpose  of  the 
current  BNL  Light  Water  Reactor  ThermohydrauHc  Development  Program,  and  emphasied 
those  particular  aspects  of  this  program  Involving  instrumentation  development. 

It  has  been  shown  that  the  development  procedures  are  not  so  much  ones  of  new 
methods  or  new  techniques,  but  more  of  extension  of  existing  methods  and  techniques 
to  conditions  beyond  those  previously  Investigated. 

With  respect  to  the  program  as  a  whole,  the  major  test  facility  modules  have 
been  received  and  the  design  of  the  test  apparatus  is  in  progress.  Development  of 
the  methodology  to  be  used  for  obtaining  the  desired  measurement  of  non-equilibrium 
vapor  generation  rates  in  flashing  flows  has  determined  the  direction  that  this 
work  has  taken  in  both  the  hardware  and  the  analytical  aspects  of  the  program.  The 
high  velocities  and  amorphous  phase  structures  expected  to  be  encountered  have 
Indicated  to  us  that  considerable  attention  to  the  details  of  the  sensor  design 


1.3-8 


tnd  slgiMl  handling  and  analysis  methodology  had  to  be  utilized.  Jetermi nation  of 
the  system  requirements  and  specifications  for  our  long  lead  development  items  has 
occupied  a  major  portion  of  the  experimental  development  effort.  During  the  next 
six  months  It  Is  expected  that  the  actual  hardware  will  be  obtained  and  the  local 
probes  required  will  be  tested. 


DEVELOPfOT  AREAS 

SOLUTIONS 

•  (kCHANICAL  DESISN 

Alternate  Designs 

a  Electronic  Sicnal  Distortion 

Signal  Reconstruction 

a  Recordino  Incoherencies 

Cross  Temporal  Elimination 

a  Hydrodvnahic  Signal  Distortion 

Independent  Global  Reference 

a  Relate  Interface  to  Phase  Velocities 

Independent  Global  Reference 

Table  1:  Development  Area  for  Local  Probe  Design 
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mr  I:  OEVUOMCItT  Of  COlSTITUTlyt  hclatims 

Talk  I.  1  AMl/t<ca1  ilMaltiig  •  kM-CquOibrliiB  Vkwr  GtntratlM 

•  CffkCU  «f  HMt  F1i»  •  Pott  ChF 

•  CFFoct*  of  Prttturo  ChMQt  -  FlkthtfiQ 

Ttfk  1.2  Dttlfn  and  CoMtructlen  of  Fsporlotoul  Cowipaont 
0  Tait  Facility 
a  Tott  Sactlon 
o  ftlobal  DcMltoaotar 
Task  1.3  Eiparlaantal  Cvaluatloo 

Voluactric  Vapor  (ianorattoo  In  Flatninp  Stcao>tfatcr  FIoms 
o  Steady  State 
0  Trantlent 


PMT  !1:  OCVELOFMENT  OF  INSTRUKtKTATIOH  TtCHNigUFS 
Task  II. 1  Analytical  Hoaellng 
0  Local  Initruaentt 
0  Global  Instrunents 

Task  II. 2  Deslon  and  Construrtlon  of  Cxperloental  Eqjipwnt 

0  Local  Void  Fraction 

a  Local  Phase  Velocity 

0  Global,  LOFT/Seniscale-Type  Instnments 

Task  Ii.3  Experimental  Evaluation 

a  Local  Instrunents  (Part  of  Task  1.3) 

e  Global  Instrunents  -  Provide  Calibration  Techniques 
(Conblned  with  Task  11.1) 

a  Steady  State 

a  Transient 


Figure  1:  Program  and  Task  Overvievf 
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:?«einat1c  of  the  Test 


Section  in  Axial  Cross-Section 


Figure  2 


o  0.1  az  as  a*  as  ojt  ar  a*  o*  u> 

VOID  rNACriOM-O 


Relative  Error  In  Quality  Determination  Due  to  Errors  In 
Determining  Parameters 


—  MAIN  lines 

- INTERWCDIATE  lines 

- SMALL  Lines 


Figure  3:  Schematic  Diagram  of  the  Test  Facility 
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ABSTRACT 

The  ultiaate  aia  of  this  experiaental  investigation  is  to  study 
droplet  hydrodynaaics  under  conditions  siaulating  steaa  water  droplet 
flow  predicted  to  occur  during  the  reflood  portion  of  LOCA.  This 
research  will  support  the  developaent  of  advanced  systeas  code  for 
which  constitutive  law*  aust  be  developed  ultiaately  for  predicting 
heat  transfer  during  core  rcflooding  in  the  vicinity  of  the  rising 
water  level  where  droplets  are  foraed  by  violent  boiling  and  broken 
up  by  the  Leidenfrost  phenoamon.  The  technique  involved  the  utiliia- 
tion  of  SUKY-SB  developed  Laser-Doppler  aneaoaeter  and  experiaental 
hardware  for  flow  aeasurcaents .  In  order  to  do  this,  an  entirely  new 
aethodology  and  the  accoapanying  special  electronic  instruaentation 
would  have  to  be  developed.  To  test  and  verify  this  new  scheae,  a 
solid  particle-water  flow  systea  has  been  used.  Soae  preliainary 
results  are  presented. 
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nmooucTiQN 


Although  tho  use  of  leser-Doppler  sneaoaetry  in  single  phsse 
flows  has  been  established  for  soae  tlae.  Its  application  in  two-phase 
dispersed  flow  still  presents  considerable  probleas  in  obtaining 
statistics  on  such  flow  properties  as  particle  size  nuaber  density 
and  velocity  distributions.  One  of  the  probleas  in  particle  aeasure- 
■ent  is  the  a^iguity  in  the  Doppler  signal  due  to  the  presence  in 
the  aeasuring  voIum  of  aore  than  one  scattering  particle  at  the  sane 
tiae.  If  the  dispersed  phase  is  reasonably  dilute  in  a  usually  extreaely 
saall  optical  aeasuring  'olusie,  such  an  aabiguity  can  be  readily  reaoved. 

A  far  aore  serious  problea  is  the  aadiiguity  originating  froa  the 
inherent  nonuniforaity  of  illuaination  in  the  optical  measuring  volxiae. 

In  order  to  extract  inforaation  on  sizing,  in  addition  to  inforaation 
on  velocity  of  particles  froa  Doppler  signals,  it  is  desirable  to  devise 
a  scheae  by  which  a  central  core  region  of  the  optical  aeasuring  voluae, 
where  the  illuaination  is  aore  unifora,  can  be  isolated  out  to  serve  as 
a  new  probing  voluae  in  the  flow. 

Laser-Doppler  Aneaoaetry 

Laser  aneaoaetry  is  a  technique  which  utilizes  scattered  light 
froa  particulates  in  a  fluid  to  aeasure  the  velocity  of  that  fluid. 

In  theory  the  laser  aneaoaeter  is  linear,  needs  no  calibration,  and 
aeasures  velocity  independent  of  fluid  properties.  A  relatively  small 
measuring  voltaae  and  inherently  fast  response  give  its  ability  to 
follow  rapidly  changing  velocities  in  the  fluid.  Only  light  needs  tc 
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enter  the  fluid  to  the  Mesurins  point.  General ly  tuo  basic  aodes  of 
operation,  e.f.  dual  beta  and  reference  beaa  are  used. 

Dual  beaa  is  the  crossing  of  two  laser  beaas  of  equal  intensity 
at  a  point  in  the  fluid  to  be  aeasured.  Mhere  these  beans  cross,  they 
interface  with  each  other  to  fora  'fringes'.  A  particle  noving  through 
the  crossing  point  in  the  plane  of  the  two  beans  then  goes  through  a 
region  of  very  low  light  intensity  (light  cancelling)  to  a  region  of 
high  intensity  and  back  again.  A  photodetector  is  used  to  pick  up  the 
scattered  light  fron  the  particle.  The  frequency  of  the  signal,  so 
obtained,  is  proportional  to  the  velocity  with  which  the  particle  is 
noving  across  fringes. 

Reference  beaa  operation  is  generally  described  in  tens  of 
Doppler  shift  of  the  scattered  light  and  it  requires  only  one  light 
beaa  at  the  neasuring  point.  The  scattered  light  fron  a  particle  is 
aixed  with  a  reference  beaa,  to  detect  the  shift  in  frequency.  The 
photo-detector  responds  to  the  difference  in  frequency,  which  is 
proportional  to  the  velocity  of  the  particle. 

In  both  operations,  the  shape  of  the  neasuring  voluae  is 
approxinately  ellipsoidal  and  is  as  shown  in  figure  1.  The  laser 
beau  generally  have  a  Gaussian  intensity  distribution.  This  causes 
the  aiqplitude  Variation  shown,  with  naxiaua  anplitude  occurring  in 
the  center  of  the  neasuring  voluae.  The  three  aeasurable  characteristic 
paraneters  of  a  Doppler  signal  are  as  shown  in  figure  2. 
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Iht  frlngt  spacing*  and  other  Inharaot  nature  of  dual  bean 
operations  Halt  the  application  of  tills  systaa  only  to  very  snail 
particle  sixes. 

The  amplitude  of  Doppler  slpial  Is  nalnly  a  fwietlMi  of  particle 
site,  nuaber  of  particles  In  nsasurlng  volua*  and  particle  length  in 
the  nsasurlng  voliaM.  Most  neasureaonts  In  two-phase  suspensions  are 
restricted  to  soae  rather  over-siaplified  cases^'^.  Considerable 
effort  has  been  devoted  to  this  difficult  tash  of  precise  siting  and 
velocity  noasureaent  of  individual  particles  by  concentrating  on 
particles  of  large  enough  site  for  the  technique  of  geoMtrlcal  optics 
to  be  rendered  useful^  and  by  asking  use  of  an  additional  property  of 
the  signal  '  ,  for  exaaple.  One  additional  serious  problen  in  signal 
aablguity  originates  froa  the  inherent  non-imifoxalty  of  illuaination 
in  the  optical  aeasuring  voluae*.  Soae  effort  has  been  atteopted  to 
by-pass  this  difficulty,  for  Instance,  by  incorporating  a  second  optical 
systea  designed  for  particle  sizing  which  is  focussed  around  the  saae 
preoing  point^^.  It  still  reaains  desirable  to  devise  a  scheae  by  which 
a  central  core  region  of  the  aeasuring  voluee  where  the  illuaination 
is  aore  imifora.  can  be  isolated  out  to  serve  as  a  new  probing  voluae 
in  the  flow. 

Theory  of  the  proposed  scheae 

Ve  selected  the  reference  beaa  node  of  operation  as  this  can 
be  extended  for  larger  size  particles.  We  decided  to  operate  in  a 


411iie«  dicpcTMl  flam  tfttm  mU  hmc»,  tkm  slpul  Mi^lltuia  will  mm 
b«  Mialx  a  Aaictioa  of  portlclo  also,  4,  and  particlo  path  laagth  U 
tha  aaasttiiag  woIom.  Ihls  rtatrictioa  of  diluto  flow  to  g&t  mm 
particlo  la  tha  aaaauriag  voltaH  ia  aot  aorlewa  aa  oar  ■aaawrlng 
voluat  ia  ao  aaall.  tha  aatl—  allowahlo  atahor  doaaitioa  poaaiblo 
ara  quit#  largo  (particloa  caa  havo  a  awabor  doaaity  of  op  to  about 
10*/cc). 

If  wo  focus  our  attoatioa  thea  to  oaly  coatral  coro  of  tha 
aaaauriag  voliow,  a  vary  difficult  foot  ao  far  aot  obtalaad  by  othars, 
tha  aigaal  aaplituda  will  dapoad  oaly  oa  tha  sixa  of  tha  porticla,  d. 
This  caa  ba  aehiavad  by  alactroaicailly  iaolatiag  a  caatral  oora  of 
aaaauriag  woluoe,  Figura  S. 

Tha  path  laogth  of  a  scattariag  particle  through  this  aaaaur- 
log  wolvow,  as  showa  ia  Figura  3,  caa  ba  detesoiaad  froa  tha 
product  of  tha  tloo  of  duratloa  of  tha  sigaal,  t,  oad  the  velocity 
of  the  particle,  v^.  which  is  dataroinsd  by  the  Doppler  froqueacy  of 
tha  signal  itself,  ty*Tv.  By  then  ignoring  all  signals  whose  corres¬ 
ponding  particle's  path  length  is  soaller  than  a  chosen  lower  Uniting 
value  2^,  we  con  essentially  eliaiaate  the  dependency  of  signal 
a^ilitude  on  particle  position  in  the  aeasuring  voluae.  The  aasiaiai 
path  length,  (the  one  going  through  the  geonetrical  center  of  tha 
aeasuring  voluae)  is  obtained  eopirically,  by  detemining  the  path 
lengths  of  aoay  particles  through  calibraticn.  The  lower  liaiting 
value  of  the  path  length  for  tha  central  core,  I^,  is  deteraiaed  by 
the  choice  of  a  window  sixe  for  tha  aaplituda  discriaiaation. 


After  tipwl  aaplltu4«  disariaiaatioB,  •  discriaiaetioa 
•isiast  tk#  path  laagth,  with  the  llaitiaf  path  laagth,  of 
tha  eontral  cora  of  tha  ■aasuriag  voliaa  aa  tha  lowar  bound,  only 
signals  within  tha  salactad  anplituda  window  from  particles  passing 
through  tha  central  core  will  be  allowed  to  go  throu^.  Then  an 
additional  discrinination  on  velocity  will  produce  required 
statistics  on  the  particle  site  and  velocity  distributions  in  tha 
flow  of  a  suspension. 

By  setting  a  lower  Unit  on  the  tine  duration  T,  say  T,  using 
a  tine  discrinination  sche*,  and  a  related  lower  Unit  on  the  velocity 
Vj,l^/T,  using  a  velocity  discrinination.  we  will  retain  signals  within 
the  velocity  range  which  can  only  be  produced  by  particles 

going  through  the  central  core  region  of  the  neasuring  volune, 
as  shown  by  the  sketch  of  Figure  4. 

However,  this  is  not  to  say  that  these  signals  represent  the 

signals  produced  by  all  the  particles  of  the  specified  size  d|^ 

(corresponding  to  signals  of  aaplitude  A^)  which  pass  through  the 

central  core  region,  because  .he  si|^als  fron  those  particles  which 

have  a  velocity  above  the  upper  linit  of  the  local  adnissable  velocity 
£ 

range,  that  pass  through  the  central  core  region,  will  not  be 

T 

registered.  In  order  to  obtain  the  correct  nunber  count  appropriate 
for  the  detersdnation  of  the  velocity  probability  density  distributions 
and  particle  size  niaber  densities  fron  the  neasured  nunber  coimt  so 
obtained,  a  correction  factor  K,  the  ratio  of  actual  nuaber  count  to 
correct  number  count,  has  to  be  detemined.  This  correction  factor 
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is  strlctljr  s  tactioa  of  fsoastry  only  and  can  bn  doterninnd  fion 

cnllbrntlon  with  a  flow  of  known  mabor  donaity  (N.)  and  a  laiifpxn 

'  o 

pnrticlo  siso  d.  with  tha  aaaa  voloeity  (v.)  . 

*  Jo 

bat  N  .  ■  no.  of  paxticlas  of  sito  d.  per  tmit  velvmt,  per 
^  •  J  ^ 

(Bit  velocity  raafo  and  laie  site  ranse  with  the 
canter  velocity  Vj. 

Ihen 


(1) 


and 

>*i  •  I  >*i  *  •  Z  ^  (2) 

1  ^  l.J  j 

where 

■  mni»er  density  of  particles  of  size  d^. 

But  for  the  over>siaplified  calibration  nai,  the  value  of 

(Bj  j)  aiaply  stands  fbr  the  mniber  count  of  particles  of  size  d^ 

having  velocity  (v.)  per  laiit  cross-sectional  area  of  flow.  By 
J  o 

comparing  this  to  tho  total  nuhber  count  for  the  whole  central  core 

regioB,  obtained  by  a  schene  of  gradually  raising  the  lower  linit  on 

tho  dumtion  tine  discrinination,  t,  the  value  of  the  cross- 

sectional  area,  can  be  obtained,  since  i  is  siaply  the  ratio  of 

o 

this  total  ma*er  count  to  the  maibor  count  per  unit  cross-sectional 
area  of  the  central  core,  (n.  .}  .  For  a  aeasuring  volune  approainately 

o  .2 

elliptical  in  shape,  the  plot  of  the  nuaber  count  a  against  (t)  should 
be  Bade  up  of  two  approainately  straight  segnents  as  shown  in  Figure  S. 
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It  is  of  iatorost  to  not#  that  sinco  the  nuiber  count  n  end 
the  path  length  1^,  for  this  over-sinplified  case,  are  proportional 
respectively  to  and  t,  the  correction  factor  K,  needed  for 
converting  the  actual  niari>er  count  to  corrected  nunber  count,  can 
be  readily  obtained 


-  (2t  -t  f  -f*J 
£  ■  n  V  >r 

^  (T  -T*) 

'  n  v' 


(3) 


Instrunentation  and  Experinental  Apparatus 

Instead  of  neasuring  the  aoplitude  (A^),  tine  duration  (t) 

and  velocity  (v^]  of  each  signal,  the  output  of  the  LOA  is  filtered 

as  shown  in  Figure  6.  After  an  anplitude  discrinination  set  on  signal 

anplitude  A^  and  a  tine  discrinination  set  on  the  lower  Unit  on  the 

signal  path  duration  tine  f,  only  the  signals  having  an  anplitude 

within  the  pre-selected  narrow  anplitude  window  in  the  neighborhood 

of  A^  and  a  duration  tine  t  i  t  car  reach  the  tracker.  The  tracker 

then  converts  the  Doppler  frequency  of  the  signals  into  voltage  signals 

proportional  to  the  ~'<rticle  velocities.  The  output  of  the  tracker  is 

passed  to  a  velocity  discrininator  cl.  .uit  which  allows  througti  only 

those  signals  with  a  voltage  (correspcndiiig  to  the  particle  velocities) 

above  a  certain  value  (equal  to  This  results  in  signals  for  which 

f 

the  product  Vj  •  t  >  is  a  value  greater  than  K^,  and  hence  only  for 
those  particles  passing  through  the  central  core  region.  The  nuBi>er 
count  of  such  signals  is  executed  on  an  electronic  counter  to  give  the 
■ensured  nui6er  count  of  particles  of  a  particular  size  and  a  particular 
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vtlocity  as  pravlously  dtscribad.  By  dianginf  tha  tiaa  discriaiastor 
lovar  liai£  t  and  tha  corraspaadlnf  valocity  dlacriniaatar  lowar  Unit 

K 

nunhar  caimts  at  different  velocities  can  be  obtained  for  particles 

T 

of  tha  saae  site  ranfe.  Than,  by  systeaatically  changing  tha  setting  of 
the  aaplitude  discrininator  and  repeating  tl:e  aforenentioned  tine  and 
velocity  discrinination  schenes,  statistics  for  the  conplete  particle 
size  range  can  be  readily  accoaplished.  Furthemore,  by  analyzing 
sipials  froa  particles  of  snail  enough  size,  the  velocity  probability 
distribution  for  tha  fluid  phase  can  also  be  obtained,  using  these 
snail  particles  as  tracers.  An  instruaentation  block  diagran  is  shorn 
in  Figure  7. 

The  operational  arrangenent  of  our  reference-node  Laser-Doppler 
Anenonetar  is  shown  in  the  sketch  of  Figure  8.  The  onconing  laser  bean 
from  a  ISan  He-Ne  Ol  laser  is  split  into  two  beans;  the  scatter  bean 
which  foms  an  angle  of  7.5*  with  the  transverse  axis  in  a  horizontal 
plana,  and  the  reference  bean.  The  neasuring,  scattering  volune  is 
approxinately  ellipsoidal,  250y  in  dianeter  and  800ii  in  length.  The 
receiving  lenses  are  synnetrical  with  the  sending  lenses,  and  the 
scattering  angle  is  fixed  at  noadnally  15  degrees.  The  received 
scattered  bean  fron  a  noving  scattering  body  passing  through  the 
scattering  voliaw  is  brought  to  nix  with  an  mshifted  reference  bean 
to  generate  a  heterodyne  on  the  photonultiplier  tube  surface  to 
produce  the  Doppler  frequency  shift  signal. 

The  flow  apparatus  consists  nainly  of  a  precision  close-loop 
water  channel  which  has  an  effective  length  of  6.10a  and  inside 
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diamsicns  of  45.7ca  wide  and  30. Sea  deep  with  walls  nade  of 
optical  srade  plate  glass.  The  entire  optical  system  is  mounted  on 
top  of  a  heavy  traversing  base,  switching  and  aligned  with  the  water 
channel,  which  is  capable  of  positioning  the  optical  measuring  voliose 
inside  the  water  channel  accurate  to  25y  in  all  three  directions. 

The  water  channel  receives  its  water  supply  from  a  Hiexonite  filter 
with  a  maximum  pore  size  of  3)j.  Precision  neutrally  buoyant  hollow 
glass  spheres  of  sizes  in  a  distributed  range  of  up  to  3Sy  are  added 
to  the  channel  water  loop  to  serve  as  the  suspended  phase.  The  Doppler 
signals  from  the  optical  anemometry  system  are  first  recorded  on 
magnetic  tape  which  is  afterwards  played  back  at  a  slower  speed  for 
data  procession  in  the  instrumentation  system. 

Calibration 

According  to  the  calibration  scheme  as  shown  in  Figure  5, 
particles  of  one  single  precise  size  are  needed,  a  feat  not  easily 
attainable  because  of  the  difficulty  of  obtaining  such  particles.  A 
substitute  way  has  been  found  in  experimenting  with  very  dilute  sus¬ 
pension  of  particles  of  distributed  sizes. 

The  optical  measuring  volume  was  placed  in  that  portion  of 

the  water  channel  where  a  laminar  unif-  rm  free  stream  flow,  (v.)  » 

^  o 

6.43  cm/sec,  was  established  as  shown  in  Figiire  9.  In  a  preliminary 
trial  run,  the  amplitude  discrimination  window  was  deliberately  made 
very  narrow  to  scan  over  the  whole  amplitude  range  to  obtain  a  rough 
picture  of  the  signal  amplitude  distribution  pattern.  Due  to  the 


dilution  of  portieios  nd  tho  vozy  norrow  uplitudo  window  sist  uood, 
tho  tifnnl  uplitudo  diaployod  noor  disczoto  sito  diotribution  tt  • 
vory  fow  iaolotod  oqilitudoa.  Ono  of  thon,  ■  770*110  nv  was 
finally  adoptod  and  its  suitability  for  this  task  is  daaonstratod 
in  tho  wall  bohavad  calibration  curva  of  signal  nunbar  count,  n,  vs. 
square  of  lower  limit  on  duration  tisw,  t*,  as  shown  in  Figure  10. 

The  duration  tine  based  on  tape  speed  of  IS  inches/sec.,  is  twice 
as  large  as  t,  based  on  the  recording  speed  of  SO  inches/sac.,  and 
the  operational  amplitude  window  size  O.OSA^.  The  total  tape 
length  is  60S  ft.  The  two  important  quantities  characterizing  the 
central  core  of  the  measuring  volume  for  this  case  are  found  to  be: 
ly  -  228u 

-  267m. 

And  the  correction  factor  needed  for  converting  the  actual  number  count 
to  corrected  number  count,  according  to  Eqn.  (S),  has  the  value: 

X  «  0.536. 

The  next  step  is  to  esta’ lish  a  calibration  between  the  signal 
amplitude  and  particle  size  from  the  signal  number  count  and  particle 
number  density  from  Independent  measurement  of  the  water  sample.  The 
accumulative  actual  signal  number  count,  [n],  for  signals  of  amplitude 
firom  A^  to  A^^,  the  maximum  amplitude,  for  the  various  signal  ampli¬ 
tudes  is  plotted  in  Figure  11.  The  curve  gives  the  limiting  value: 

Inl  ■  1.  at  A^^^  •  1.08v. 


1.4-11 


V«t«r  ••■pi*  taken  lHM41atcl)r  downatinaa  of  tha  iwaauriaf 
velUM  with  an  iseklnatlc  probe  was  analyted  with  a  Coulter  (Model  B) 
counter.  The  resulting  accunulative  particle  nunber  density  count, 
(N],  for  particles  of  size  fron  d^  to  d^^,  the  nariniM  size,  for  the 
various  particle  sizes  is  plotted  in  Figure  12.  The  curve  gives  the 
limiting  value: 

(N]  •  1.  at  d^  •  SS.4U 

ly  definition,  we  have  for  the  usiifon  flow. 


'  "lOj).  •« 


k  N. 


%ihere  n^  ■  actual  signal  niaber  count  of  sip  '*f  aaplltude 

K  ■  correction  factor  for  signal  ntod" 

■  nuMber  density  of  particles  of  siz 

(v.)  ■  unifom  nixture  velocity 
'  o 


■  base  cross-sectional  area  of  central  c 
o 

vol 


of  neasuring 


■  tape  recording  tine 

k  ■  K(v,)  6^  T_,  a  omstant 

Jo®® 

Therefore,  we  have  the  relationship 
[n]  -  k(NJ 

between  the  two  accunulative  quantities.  Furthemore,  if  we  let 
A*  -  Aj/A^  and  d*  •  dj/d, 


we  can  show  that  by  probability  theory 


*n[Nl 


in  In] 


at  d* 


at  A* 
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vliicb  MtablitlMs  •  wilqut  comlatlM  mlatiauUp  bctwMB  d*  nd  A* 

M  slNMi  hif  cfliviitiaMl  ichi  of  Pifuro  IS.  Hm  tofultMt  colibntloa 
curvo  btwi  tho  BonMliiod  portlclo  tlso  d*  and  tipMl  uplitudo  1» 
ihowB  in  Pifuro  14.  Rasult*  froa  diract  a^plituda  aaasuraMnt  uaiag 
particlas  of  fivo  paaoisa  slzo  raagas  ara  also  plottad  for  coaparison. 
Tho  diaansioBal  calibratioa  cunra  botwoaa  particla  slzo  and  sliaal 
aaplitudo  is  plotted  in  Figuro  IS.  Using  this  calibratioB,  tho  rasul- 
tant  particle  size  and  ntabor  density  distribution  in  this  laninar 
unifom  free  stroaa  flow  is  plotte>'  in  Figuro  16. 

Furthart  for  the  previous  calibration  run  for  tho  dotoxnination 
of  ly,  and  K  as  shown  in  Figure  10,  we  have 

(•),  •  CN,)  (..)  »,  I. 

''  o  o 

whore  (n)^  ■  18,  tho  ^.tabilizod  ntad>or  count 

(N.)  *  20/em^t  tho  nuabor  density  of  particles  of  sizes  between 
*  0 

10. Su  and  12. Su,  corresponding  to  amplitudes  of 
770  anr  and  810  nv  respectively,  throu^  the  cali¬ 
bration  curve.  Figure  IS,  at  IS/cn^  -  |i  froa  the 
nunber  density  distribution  curve.  Figure  16. 

(v.)  ■  6.43  cn/sec,  the  wifom  firee  stroaa  velocity 
^  o 

■  484  sec.,  the  tape  recording  tine, 
and  therefore,  the  base  cross-sectional  area  of  central  core  of 
aeasuring  voluae  is  then: 

■  .  _  .a4.  2 


A  MMSitriaf  point  wu  Mloctod  on  the  edge  of  the  wake  behind 
a  tranivene  soctangiilar  body,  as  showi  tn  Figure  9,  where  significant 
particiilate  accunulations  and  sitsbie  velocity  fluctuations  are 
expected.  The  particle  nuiber  coiaiting  rate  per  unit  cross-sectional 
area  per  wit  velocity  range  and  wit  site  range  was  obtained  fro* 
a  length  of  tap*  of  32S  ft.  recorded  at  SO  inches/sec.  for  two 
representative  particle  site  ranges  (4.60-4.9SW  and  S.6S-6.1Qu).  The 
velocity  distributions  for  these  particle  ranges  are  shown  in  Figure  17. 

Concliision 

A  new  experiaental  scheae  has  been  developed  for  the  neasure- 
iwnt  of  the  local  nuaber  densities  and  velocity  probability  densities 
of  a  dilute  two-phase  suspension  which  has  a  distribution  of  particle 
size  and  a  predoainate  direction  of  flow  orientation  by  the  use  of 
laser  Doppler  aneaoaetry.  An  accompanying  indirect  calibration  aethod 
needed  for  this  scheae  has  also  been  developed  based  on  the  probability 
correlation  between  Doppler  burst  cowts  and  the  independently  obtained 
aixture  saaple  sizing  nudber  cowts. 

Results  of  aeasureaents  using  this  scheae  have  been  presented 
for  a  laainar  free  streaa  flow  and  a  turbulent  shear  flow  of  a  dilute 
two-phase  suspension  of  neutrally  buoyant,  spherical  glass  particles 
of  distributed  sizes  in  water.  Reproducibility  of  these  results  has 
been  coafiraed. 


nu*  mak  kM  hmm  4mm  at  tha  Stata  ttiivanltjr  of  Man  York 
at  Staay  Iroak  mdmr  SabcaatTact  Na.  S60St4.S  tor  tko  Iroakhavan 
Natieaal  Irtaratery  of  Caatraet  No.  AS0t4(0SSS4)  fbr  tha  Uiita4 
Statas  Nuelaar  Ragulatary  CoaBiaaiaa.  Ta  thaaa  afaaciaa,  tha 
aatlwra  aould  lika  to  aiqpvaaa  thair  daa^  appraciatiaa. 

Tha  authera  alao  ackaowladca  Mr.  Philip  Pritchard  fbr  hia 
caatribiitiaa  ia  tha  davalopiag  atagaa  of  thia  aork. 
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F16.  2.  SKETCH  OF  PARTICLE  DOPPLER  SIGNAL 
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FIGURE  A  OPERATIONAL  DLOCR  DIAGRAM  FOR  SIGNAL 
0ISCRIN1NAT^«3M  SCHEME 


FIG.  7.  INSTRUMENTATION  BLOCK  DIAGRAM 


FIGURE  0-  OFTICAL  ARRANGEMENT 


riCURE  •  .  FLOW  OF  A  TWO-PHASE  SUSPEJiSIOH  AROUHO  A  TRANSVERSE 
RECTANGULAR  BOOT 
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FIG.  ii|.  couparison  of  sighal  amplitude  calibration  results  FROr 

NUriBER  COUNT  CORRELATION  AND  DIRECT  MEASUREIIENT  USING 
PARTICLES  OF  PRECISE  SIZE  RANGE 
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FIG.  17.  VELOCITY  DISTRIBUTION  FOR  PARTICLES  OF  SELECTED  SIZE  RANGES 
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1.  IimMOUCTlttl 


TIm  Loss-of>F1u1d  Testing  (LOFT)  nxptrlMnt  has  be«n  dnslgnad  to  do 
safety  testing  for  the  Hater  Reactor  Safety  Agency.  The  objectives  of 
the  experlaantat  prograe  are^^^. 

(1)  To  provide  data  required  to  evaluate  the  adequacy  and  Inprove 
the  analytical  Mthods  currently  used  to  predict  the  1oss>of> 
coolant  accident  (LOCA)  response  to  large  pressurized  uater 
reactors  (LPHRs).  The  perforeance  of  engineered  safety 
features  (ESF)  with  particular  enphasis  on  the  eaergency 
core  cooling  systee  (ECCS)  and  the  quantitative  Margins  of 
safety  Inherent  In  the  perfonMnee  of  the  ESF  ere  of  prlnary 
Interest. 

(2)  To  Identify  and  Investigatz  any  unexpected  event(c)  or 
threshold(s)  In  the  response  of  either  the  plant  or  the  ESF 
and  develop  analytical  techniques  that  adequately  describe 
and  account  for  such  unexpected  behavior. 

LOFT  Is  presently  In  the  prenuclear  checkout  phase  and  has  made 
three  nonnuclear  blowdowns  to  date.  The  flow  Instrumentation  Is  of 
major  Interest  In  measuring  the  blowdown  results  since  the  flow  measure¬ 
ment  parameters  are  of  prime  Importance  In  the  code  predictions. 

One  of  the  major  purposes  of  the  nonnuclear  test  series  Is  to 
determine  that  the  equlpment/systems  function  properly.  This  Includes 
the  experimental  Instrumentation  of  which  the  flow  Instruaentatlon  Is  a 
part  and  Is  the  topic  of  this  paper. 


LOFT  Is  •  SS  (mt)  prtssurlitd  Mttr  mctor  fsclllty  Inttndsd  to 
slouUto  tho  Mjor  bthivlorsi  ssptcts  of  gontrfc  1000  (IWt)  LPyRt  In  « 
carefully  cotMkicted  Loss-of-cooUnt  oxporloont  (LOCE).  Tho  imctoir  core 
Is  opproxlntoly  S-1/2  ft  long  ond  2  ft  In  diaaoter  and  contains  1300 
fuel  pins  and  four  control  assaoblles  of  typical  LPyR  design.  The 
prloary  coolant  systoo  Is  designed  with  a  slollar  prloary  syston  voluoe* 
to«core  power  ratio  as  exists  In  typical  LPWRs. 

The  unbroken  LPUR  reactor  coolant  loops  are  sloulated  by  the  single, 
unbroken  circulating  loop  In  the  LOFT  prloary  coolant  systen  (PCS), 
operating  loop.  Figure  1.  and  the  postulated  broken  LPM  loop  Is  sim- 
latcd  by  the  LOFT  blowdown  loop.  Figure  2.  Qu1ck*open1ng  valves  (with 
opening  tines  adjustable  fron  approxinately  10  to  50  nsec)  sinulate  the 
Initiation  of  prloary  coolant  piping  ruptures.  Prloary  coolant  blowdown 
effluent  Is  collected  In  a  blowdown  suppression  tank  which  can  nodel  the 
significant  portions  of  the  various  LPW  containnent  back>pressure 
transients.  An  energency  core  cooling  systen  (ECCS)  Is  provided  to 
nodel  the  1oss*of*coo1ant  ESFs  In  LPyRs.  The  principal  conponents  of 
the  reactor  plant  are  located  on  the  noblle  test  asscnbly  (NTA). 

Heat  Is  transferred  fron  the  PCS  to  the  secondary  side  of  a  vertical 
U-tube  stean  generator.  Stean  generated  at  approxinately  750  psig  Is 
fed  to  an  air-cooled  condenser  and  Is  condensed  at  approxinately  300 
psia. 

Fluid  pressure,  tenperature.  velocity*  and  density  are  nonitored  at 
key  points  at  the  reactor  systems  by  extensive  Instrunentatlon.  Core 
themocouples  are  provided  to  monitor  fuel  pin  clad  temperatures  and 
support  tube  temperatures  at  In-core  locations.  Fixed  nuclear  detectors 
and  the  traversing  In-core  nuclear  detector  systen  measure  core  nuclear 
response  and  neutron  flux  shapes. 
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Fig.  1  Operating  loop  vltli  1nstn»ent  locations. 


Of  prlaary  Interest  to  us  for  this  discussion  are  the  flow  locations 
There  are  three  Instrunent  penetrations  in  the  PCS  (See  Figure  1)  and 
t«o  penetrations  In  the  blowdown  1o(^  (See  Figure  2).  At  each  of  these 
locations*  there  Is  a  drag-disc  turbine  transducer  (OTT),  a  fluid  therwo- 
couple.  an  absolute  pressure  transducer  and  a  ga— a  densitoneter.  There 
are  also  five  on(s)  located  Inside  the  reactor  vessel  for  nuclear 
testing  to  provide  flow  Infonaatlon. 


Tht  drag>41sc  turbine  trensduccr  (OTT)  essentially  is  a  drag  ■eter, 
a  turbine  flOHeeter«  and  a  therwcouple  coebined  in  a  single  body  and 
■ountad  in  series  as  shown  in  Figures  3  and  4. 

The  drag  disc  is  aechanically  linked  to  a  variable  reluctance 
transducer.  The  associated  electronic  systea  generates  an  output  voltage 
proportional  to  the  displaceaant  of  the  drag  disc»  which  is  proportional 
^  the  force  on  the  drag  disc.  The  polarity  of  the  output  voltage 
indicates  the  direction  of  the  fluid  flow. 

The  turbine  consists  of  the  turbine  and  an  eddy  current  transducer 
which  senses  the  passing  of  each  turbine  blade.  The  electronic  systca 
then  generates  an  alternating  voltage  with  frequency  equal  to  the 
frequency  of  the  passing  of  each  turbine  blade,  and  also  an  analog 
signal  proportional  to  this  frequency.  Friction  in  the  turbine  bearings 
detereines  the  winiiui  fluid  velocity  that  can  be  weasuredi  this  is 
about  1.S  ft/sec  for  water. 

Data  frow  these  weasurewents  will  be  used  to  calculate  the  therao> 
dynawic  properties,  mss  flux,  and  mss  flow  rate  of  the  fluid  in  «diich 
the  OTT  is  installed.  Three  transducer  Models  are  available.  One  Model 
is  used  in  various  vessel  or  plenuM  installations;  the  other  two  are 
used  in  piping  installations.  The  transfer  is  constructed  so  that,  in 
the  event  of  a  transducer  Mechanical  failure,  pieces  of  the  unit  will  be 
contained  so  as  to  not  affect  the  nuclear  systCM. 
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OQNMIISQi  OF  OPBtIt  FROR  BTT 
HITN  OQBlKMETEt 
(TEST  11-1) 

? 

Ttat 

imi 

OamltOMtar  . 
1Ii^ft3(kA^) 

OTT  . 
llM/ft3  (ko/a)) 

2 

44.2  <708) 

47.6  (762) 

4 

41.0  (657) 

43.9  (703) 

48.4  (ns) 

42.4  (679) 

8 

35.8  (573) 

40.9  (655) 

33.2  (532) 

40.2  (644) 

12 

32.0  (513) 

37.3  (597) 

14 

26.0  (416) 

35.1  (562) 

16 

22.1  (354) 

30.1  (482) 

18 

18.0  (288) 

27.0  (432) 

20 

16.8  (269) 

25.3  (405) 

22 

12.0  (192) 

26.0  (416) 

24 

9.5  (152) 

25.4  (407) 

26 

7.0  (112) 

23.7  (380) 

28 

4.8  (77) 

19.6  (314) 

30 

4.0  (64) 

15.3  (245) 

32 

3.0  (48) 

8.1  (130) 

j 

2.2  (35) 

2.6  (42) 

36 

1.2  (19) 

1.2  (19) 

38 

1.0  (16) 

1.0  (16) 

40 

0.8  (13) 

1.2  (19) 
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TIm  9mm  tfMHitowttr  used  fn  MiMuelMr  ttttlng  to  Mte  tho  LOFT 
dtMSlty  aMswramit  Is  as  sImm  In  Figun  5. 


(lytMCHl 


OMMWr  Calk*  (lypicail 


LOTTPl 

gnMTA 


Fig.  5  GiMi  d* 
The  Model  FN-UgaM  de 
electrical  and  Mechanical  el 
In  a  pipe. 


densItOMeter  systoM  IsoMetrIc. 

densItOMeter  systen  contains  all  the  necessary 
elenents  for  Monitoring  the  density  of  a  fluid 


The  basic  eleoents  of  the  systcM  are  (1)  Cs>137  gaaMS  source  hinised 
In  a  shield  cask,  (2)  three  each.  Integral  scintillation  detectors,  (3) 
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ttHTM  Mdit  tliltidad  casks  fbr  tht  daUctors.  (4)  three  each,  high 
voltage  poMor  supplies*  (S)  three  each,  preaapllflers,  (C)  three  each, 
aepllfler/control  units  and.  (7)  an  actuator  control  panel. 


The  Cs>137  source  consists  of  30  Curies  of  cesim  chloride.  Cs-i37 
Cl.  The  source  Mterlal  Is  pelletized  and  a  shio  Is  used  to  prevent 
■oveaent  of  the  pellet  In  the  source  housing. 

The  gams  source  lead  casks  Unit  radiation  to  less  than  50  nr/hr 
at  contact  with  the  source  In  the  stored  position  and  neets  the  U.S. 
Oepartnent  of  Transportation  (DOT)  regulations  as  an  approved  shipping 
container.  The  cask  Is  designed  to  be  rigidly  nounted  on  a  pipe  fixture. 
The  unit  Is  8  In.  In  dianeter  and  18  In.  long  exclusive  of  the  aounting 
device.  An  air  operated  source  positioner  with  a  fall-safe  (stored) 
nechanisa  provides  reaote  control  of  the  source.  An  electrically  oper¬ 
ated  solenoid  driven  latch  asseably  Is  also  provided  to  lock  the  source 
In  the  exposed  position. 

Each  scintillation  detector  consists  of  a  2  x  2  In.  Nal  (thalllua 
activated)  crystal  and  a  photoaultipller  tube  aounted  as  an  Integral 
unit. 

The  scintillation  detectors  are  Installed  In  shielded  detector 
casks*  see  Figure  5.  The  cask  Is  16  In.  long  by  6  In.  In  diaaeter  and 
contains  a  cooling  coll  for  water  cooling.  The  unit  Is  constructed  of 
lead-filled  stainless  steel  and  tungsten  to  provide  shielding  and 
colllaatlon.  Two  reaotely  controlled  calibration  shlas.  which  are 
attached  to  each  of  the  detector  casks  are  provided  to  slaulate  a  known 
change  In  density.  The  SHIN  LO  slaulates  a  change  of  approxiaately  10 
1bs/ft^  while  SHIN  HI  slaulates  a  change  of  approxiaately  45  1bs/ft^ 

In  the  14  In.  schedule  160  pipe  for  Bcaa  *8"  (diaaetrical  beaa).  For 
Beaa  "A*  and  *C”.  the  SHIN  LO  and  SHIN  HI  values  are  11.5,  52.3,  and 
16.2,  73.3  1b/ft^,  respectively,  beaa  locations  are  shown  In  Figure  6. 

A  NENA  box  Is  supplied  to  house  the  air  actuator  valves  that  control  the 
air  to  the  gaaaa  source  positioner  and  calibration  shia  air  cylinders. 
Also  provided  In  the  unit  Is  an  air  regulator,  filter,  and  relief  valve 
to  prevent  buildup  of  pressure  Inside  the  box. 
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Fig.  6  Thr«®  bew  ga»«  daisitoaeter. 


The  high  voltage  poi*er  supply  is  •  calibrated,  high  voluge  power 
supply  designed  for  use  with  photoewUi pH er  system  requiring  a  high 
stability,  low  noise,  well  regulated  power  source.  It  uounts  in  a 
standard  NIN  Bin  as  in  AtC  (ERDA)  Specification  TID-20893  (Revision  3). 
"Standard  Nuclear  Instnaaent  Modules."  and  occupies  two  nodule  widths. 
It  is  supplied  with  an  integral  three-wire  power  cord  and  NDiR  plug  and 
does  not  use  the  HIM  Bin  Power  Source  for  operation. 

The  high  voltage  power  supply  incorporates  over  voltage  and  arc 
protection  circuits  which  prevent  dawage  to  the  supply  and  its  load  if 


ttm  supply  ^plls  or  ooltoft  troosloots  ort  cmmd  bf  sottslooi  orcfuf  at 
tlM  loop  aoP  Its  asaociatop  wlrlop. 


Tlw  preaoplincr  sappIloP  oitii  tte  tpstao  piwiPB  curraat  to 
voltaga  conworslon  at  tM  Potortor  asiMiMy.  THs  caooortar  is  a  «ofy 
stable,  low-naiso  lorit  that  pr— iPes  am  acavate  aoP  liaoer  caooersion 
as  nail  as  pale  to  boost  aoP  Prfee  the  stpal  thramfit  Has  cable  to  the 
ai^llfler. 


The  amplifier  has  tmo  stages  (A^  and  of  —11*11.  at  too  to 
provide  flexible  ootput  neeei  Th—  — 1  i  rb  i  i  so—  stable  aaP 
10M-no1se  operatiao. 


A  potentiimtei  ZEBO  JUBUST  f-iiiPri  lo  ■tjiirtlii  0111111x1  vol- 
tape  to  the  differeetial  — 1  rfirr,  «inc»  tvsiUtts  ma  aro  shift 


(bla:*)  of  the  reLioag'  notpf  so  that  the  aaer  am  — aatch  the 
three*po1e  dcxiblethRaM  tepph.  gartch  pooviPes  *■  —  —  aoilerities  for 


the  ZfRO  s^^JilST  amitage  tprfle  m  siogie-pole  doable  thrao  aotch  grounds 
the  aapllfler  pasttive  iapec  (reference  voltage  is  swttehad  out).  A 
front  panel  neter  anartors  theoatpat  of  aatplifier  A^  —  pvavides  the 
operator  with  a  risaal  iadicatian  of  iqataa  operation. 


The  filter  seiecsar  settches  the  aatput  of  aapllfler  A2  directly  to 
the  recorder  output  or  SB  a  10  Uz  plug-in,  s1x>po1e  Bessel  filter  The 
OUTPUT  ADJUST  potcncaa— er  controls  the  full-scale  output  of  the  re¬ 
corder.  Internal  iiljB''amiHii  are  provided  for  the  panel  neter  full-scale 
indication  and  to  zero  aSl  anplifiars. 


6-  GAHWA  OEWSITOMETER  UhCERTAIUn 


The  gaana  densitoneter  system  was  thoroughly  tested  with  simulated 
flow  regimes  of  various  materials.  These  included  annular,  bubbly,  and 
stratified.  Calculations  were  made  from  the  known  densities  of  the 
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Mtcrfals  and  then  Inserted  in  a  section  of  LOFT  piping  and  aMSured 
with  the  ganM  densitoaeter  systea.  one  such  configuration  for  annular 
flow  is  shown  in  Figure  7.  The  systea  was  also  tested  with  various 
depths  of  water  (stratified  flow  regiae)  which  was  accurately  aeasured 
and  the  density  calculated  for  each  beta.  This  was  coapared  with  the 
actual  aeasure  density.  A  suaaary  of  the  uncertainty  of  each  beaa  is 
shown  in  Table  III. 


TABLE  III 

SIMWRY  OF  GAItM 
OEflSlTOHETEft  UNCERTAINTIES 

Beam 

(Percent  of  Full  Scale! 


A 

B 

C 

Statistical  Error 

0.55 

0.36 

0.88 

Nonmonoenerget i c 

0.5 

0.5 

0.5 

Stability 

0.1 

0.1 

0.1 

Calibration 

0.6 

0.6 

0.6 

Linearity 

0.2 

0.2 

0.2 

Noise 

0.4 

0.4 

0.4 

OAVOS  (OOAPS) 

0.13 

0.13 

0.13 

B— wOawaily 
ima'N'l 

'A' 

'S' 

C 

Maaiutad 

Oanatly 

2104 

14  00 

44  31 

Calcuiatatf 

Dwtsitv 

2100 

14  47 

44  20 

Fig.  7  Three  been  densItcMter  static  test  data  for  annulai'  flOM 
configuration. 


Although  uncertainty  Is  of  najor  laportance  ««e  can  also  deteralne 
flOH  reglae  by  utilizing  the  three  beans  In  a  simple  truth  diagram. 
Figure  8. 


A-B««n  112IB 
9661 


Bmh  Locations 


Fig.  8  Ganma  densitometer  flow  regime  logic. 
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7.  mA 

Tilt  OTTt  and  9«M  dcnsltOMttrs  tovt  provtn  to  bt  good  transductrs 
for  tlw  ■oasuroMfit  of  tha  flow  paraaetara  to  tha  first  thraa  oomuclaar 
tasts. 

Typical  of  tha  data  froa  tha  gaaaa  daosltoactar  Is  that  shown  In 
FIguras  9  and  10.  This  shows  tha  thraa  baaas  displayed  on  tha  saae 
graph  of  brokan  loop  reactor  Inlet  location.  Test  L1>1  as  In  Figure  9. 

It  can  readily  be  seen  that  tha  flow  Is  stratified  during  the  saturated 
blowdown  froa  the  respective  readings  of  each  of  the  three  beaas  and 
applying  the  readings  to  the  truth  table  shown  In  Figure  8.  Later  during 
CCC  Injection,  aora  stratification  Is  In  evidence.  At  the  Intact  loop 
locations  slug  flow  Is  very  prevalent  as  shown  In  Figure  10. 

MIth  this  Inforaatlon  on  flow  reglae.  It  Is  auch  easier  to  Inter* 
pret  the  OTT  data.  Since  the  OTT  Is  a  point  aeasureaent  located  In  the 
center  of  the  pipe.  It  Is  laportant  to  know  the  flow  regiae  In  the  pipe. 
Typical  data  froa  Test  L1>3A  as  recorded  froa  the  OTT  are  shown  In 
Figures  11  and  12  which  depict  respectively  the  aoaentw  flus  and  the 
velocity.  Coaplete  data  sets  are  avallable^^*^**^.  Since  the  OTT 
Is  a  point  measuring  system,  this  1$  the  primary  factor  In  the  measured 
differences  In  density  from  the  drag  disc  turbine  versus  density  from 
the  gaana  densitometer  previously  shown  In  Table  II. 
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Fifl.  10  LOFT  exptrlwnt  LI-1:  chort«l  amag*  «« 
tfansitonetar  at  InUct  loop  rtactor  Inlat  cold  leg. 


Mtl  A 


TtfntlMCI 


Fig.  11  LosS'Of-coolant  experinent  L1*3A  monentua  flux  >•  broken  loop 
cold  leg. 


•to  0  10  30  30  40  so  60  70 

TMIWIMCI  Mtt  •  <V7 


Fig.  12  LosS'Of'COolent  experlnent  L1-3A  velocity  —  broken  loop  cold 
leg. 


8.  samun 


Mt  iMv*  b«M  successful  In  atkinf  density,  velocity,  and  wentMS 
fliM  ■easurewnts  on  LOFT.  This  Is  partly  due  to  the  large  aanunt  of 
perforaance  and  calibration  testing  of  the  transducers.  Phich  additional 
separate  effects  testing  uas  also  done  to  characterize  the  flow  trans* 
ducers.  Our  present  transducer  perforaance  Is  well  docuaented  and  un¬ 
certainties  have  been  derived  for  both  single-  and  two-phase  flow.  The 
present  coablnatlon  of  the  drag  disc  turbine  and  the  aultlbeaa  gin 
densltoaeter  have  provided  very  good  correlations  of  the  data.  Re¬ 
liability  of  the  transducers  has  been  good  although  we  fee'i  It  can  be 
Improved  and  are  actively  pursuing  an  laprovcaent  progrca.  This  Is 
especially  true  In  the  case  of  the  gin  densltoaeter  which  oust  be 
changed  for  nuclear  operations. 
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TIm  MiclMr  cMMunlty  f«CM  •  particularly  difficult  problaa  ra- 
lating  to  tlw  calllratlofi  of  Instruoantatloii  In  a  tMO>phasa  flow  stoaa/ 
Mtar  anvIrcMwnt.  TIm  rationale  of  the  approach  to  water  reactor  safety 
questions  In  the  United  States  deaands  that  accurate  ■easureaants  of 
■ass  flaws  In  a  decoepresslng  two-phase  flow  be  aade.  An  accurate 
■easuroHent  dictates  an  accurate  calibration.  This  paper  addresses 
three  questions  relating  to  the  state  of  the  art  In  two-phase  calibration: 

(1)  Uhat  do  we  wean  by  calibration? 

(2)  Mhat  Is  done  now? 

(3)  Hhat  should  be  done? 

It  Bay  not  be  possible  to  answer  all  three  of  these  questions,  partic¬ 
ularly  the  last  one.  Perhaps  the  stage  can  be  set  for  thinking  and 
discussion  to  follow  this  afternoon. 


2.  MHAT  IS  CALIBRATION 


According  to  the  dictionary,  calibration  neans  to  check,  adjust, 
or  systcMtIcally  standardize  the  graduations  of  a  quantitative 
■easuring  Instruaent.  There  are  soae  key  words  to  this  definition; 

They  are:  standardize  and  quantitative. 

To  calibrate  anything  It  Inplles  a  standard  of  comparison  exists, 
and  further,  that  a  known  set  of  conditions  can  be  recreated  at  will. 

Before  we  can  have  a  two-phase  flow  calibration,  a  standard  nist  exist. 

The  title  of  this  paper  Infers  this,  but  It  has  been  submitted  that 
perhaps  this  Is  ludicrous.  Let  us  examine  the  other  key  word,  quantitative. 


Aflalv,  tl»  41ct1oMry  tells  vs,  partetefnf  to  or  sosctetlblo  to  ■Mt«ro> 
■Mt.  Vt  iMvo  ottsitelod  litro  for  tlio  purposo  of  dlscusslof  ooiv 
croitivo  and  losplMtlvo  atteopts  to  test  tlia  suscopteblllty  to  Mosora- 
■aot  of  SOM  ratlwr  alusiva  paraaatars.  A  tMo-phosa  flOH  callbratloo 
oust  aabody  som,  parhaps  all.  of  tkasa  Boasuraamt  taclMl^iias. 


3.  WHAT  IS  DOME  WOW 


Tachniquas  aoployed  currently  Involving  a  steady  state  tuo-phasa 
flow  calibration  Involve  weasurewent  of  flow  rata  of  each  coaponent  In 
single-phase  Mixing,  taking  a  heat  balance,  and  calculating  the  quality 
assuMing  equal  phase  velocities.  What  Is  known  about  these  kinds  of 
systcMS? 

The  mass  flow  In  steady  state  Is  known,  the  pressure  and  the 
teMperature  are  also  known,  provided  we  arrange  for  a  saturated  condition 
to  exist.  What  Is  not  known  Is  phase  velocity  and  distribution,  flow 
reglaes,  and  temperature.  Where  there  are  nixed  streams  of  different 
temperatures,  this  Is  of  particular  Interest  In  view  of  the  ECC  problem. 

Or.  Estrada  Indicates  that  a  flow  Instrument  which  Is  to  be  tested 
nay  actually  consist  of  several  measuring  devices  from  which  Individual 
phase  mass  flows  and  regipw  can  be  Inferred.  This  would  seem  to  Indicate 
that  any  two-phase  flows  standard  would  not  be  a  single  reference 
Instrument  but  rather  a  group  of  Instruments  from  which  a  two-phase 
standard  would  be  Inferred.  This  thinking  Is  more  or  less  universal  In 
the  Instrumentation  coemunity  charged  with  the  two-phase  flow 
measurement  problem.  The  types  of  measurements  that  are  made  are 
numerous.  A  number  of  good  reports  have  suaaarlzed  single-phase  and 
two-phase  flow  measurement  capabilities.  Among  them  are  two  EPRI  reports. 
(EPRI  NP-195  by  Brockett  and  Johnson  and  EPRI  RP-446  by  Hewett.) 

Tables  I  through  VI  present  a  broad  brush  listing  of  the  few  ^ 

two-phase  flow  loops  known  and  briefly  discuss  the  kinds  of  measurements 
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St«Mi/ttot«r  or  Alr/Votor 
Froosurot  to  22S0  pti 
Tooporoturet  to  3S0*C 
Nut  flou  to  S  kg/toc 
Horixontol  Tut  Soctloo 

Nutyrwwti  (Sloolo-Plioto) 

Inlot  StOM  Floif  (Orffico) 

Inlot  Uator  Flow  (Orlfico) 

Pruturo 

Tooporotiono 

Nutoruuts  (Two*  Photo) 

Drop  Body 

Tnio  Nott  Flow  Netor 

Radlonuclldt  Tracer 

Nuclur  Hagnotlc  Rotononco 

Drag  Body.  TurtolM.  y  DomItOMtor  CoablMtlon 

FIom  Pattom 

High  Spead  Photography 
Nultibeaa  y  Oansitoaater 
lapodancc  Prote 


boing  aodo  In  thue  laboratorlu.  Tlilt  Infonutlon  It  by  no  Mont 
Intandod  to  be  all  Inclutive.  f'^rrectlont  or  addUlont  are  Melcone. 

The  tMO-phate  flow  loop  facility  located  In  Karltruhe.  Genaapy.  It 
a  tteady  state  flow  loop.  Prusure  capabilities  to  2250  psi,  teoperatures 
to  350*C,  and  mss  flows  to  5  kg/sec  utilize  a  horizontal  tut  section. 

The  unique  future  of  this  loop  Is  that  It  Is  capable  of  also  operating 
In  air/water.  Nusuremnts  mde  are  s1ng1e>phase  Inlet  steaw  flow. 

Inlet  water  flow,  pressure,  and  teaperature.  Husuraaents  mde  on  the 
two-phase  test  section  are  drag  body  type  wusurments,  the  true  mss 
flouaieter  which  they  have  developed,  the  radio-nuclei  tracer  for  phase 
velocity  nusurewent,  nuclur  mgnetic  resounce  techniquu.  and  a 


TMLE  !I 

im-msi  FLQM  TEST  LOOP  -  CISF.  ITALY 


StMi4y  SUf  Two>Pfc»f  Flow  (So—  Trtnslgwt  toptblllty) 

SteM/Uater 
Pressure  to  3689  psi 
Tenpereture  to  300*C 
Ness  Flows  >  2S  kg/sec 
Vertlcle  Test  Section 

NeesuraMfit  (Slwole-Phese) 

Inlet  Steea  Flow 
Inlet  Ueter  Flow 
Pressure 
Teapereture 

Neesureaents  (Two-Phes^) 

Quick  Closing  Values 
lepedence  Probes 
Pressure 
Flow  Pattern 

No  special  technique  is  at  present  available. 
Heasureaent  of  local  and  Instantaneous  flow 
rates  In  two-ptese  conditions. 


coablnatlon  drag  body  turbine  gaaea  densltoaater.  For  flow  pattern 
recognition,  they  use  high  speed  photography  and  aultlbeaa  ymm  densltoa* 
eters,  plus  soae  work  has  been  done  with  «n  Inpedance  probe. 

The  two-phase  flow  test  facility  in  Italy  has  soaa  transient 
capabilities.  The  aaxlaua  pressure  i»  3689  psi.  tevveratures  to  300*C. 
and  aass  flows  In  excess  of  cl  hg/sec  yith  a  vertical  test  section. 
Neasureaents  of  the  Input  to  the  tao- phase  test  section  Include  Inlet 
steaa  fltw  and  inlet  water  flow  pressure  and  teapereture.  Unique 
aeasurrwents  within  their  two-phase  test  section  are  aade  with  quick 
closing  valves.  lapedance  probes  pressure  neasureaents  ere  also  aade. 

At  the  present  tiae.  no  special  technique  exists  for  flow  pattern 
Identification. 
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TMLE  II 


TVO-PtUSC  FUM  TEST  LOOP.  tCL.  CMMIM 


St— <|y  State  Two-Phat«  Flow 

St— M/Witer 

Pr— sure  to  900  psio 

T«^ior4tiire  to  350*C 

Hiss  i^lOM  22  kg/s«c 

Horiz— U1  or  Vortical  Test  Section 

Weasure— fits  (S1nqle«Phase) 

Inlet  Ste—  Flm 
Inlet  Water  Flow 
Pressure 
Tenperature 

M— sure— nts  {Two-Phase) 

Quick  Closing  Valves 
Hultibe—  Y  Densito— ter 
Turbl— 

Pitot  Tubes 
Flow  Pattern 
Not  Known 


At  the  WCL  Loop  In  Ca—da.  — asure—nts  are  — de  at  the  Inlet  ste— 
flow,  the  Inlet  —ter  flow,  pressure,  and  tenperature.  In  the  two-phase 
test  section,  quick  closing  valves,  —Itlbe—  ga— ■  d— site— ters, 
turbines,  and  pitot  tubes  are  utilized.  Types  of  flow  pattern  1nstru> 

— ntatlon  are  not  known. 

The  French  have  nu— rous  facilities.  A  sun— ry  of  their  total 
capability  of  the  various  facilities  Is  giv—  In  Table  IV.  Two^phase 
flow  —asure—nts  Include  ■Icrother— coupl— ,  optical  probes,  hot  wire 
and  hot  flla  ane— ters,  a  unique  liquid  b— rli^g  turbl—  -ter,  n— tr— 
be—  attenuation,  quick  closing  valves,  and  —Itlbe—  ga— a  d— sltoM- 
eters. 

During  a  recent  visit  to  Japan,  several  large  test  facilities  —re 
visited.  These  test  facilities  gave  the  Iwpresslon  that  Instru— nt 
develop— nt  Is  very  United. 
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i  .]  Winy  Loops  Aw<1«b1* 

FrciiMic.  Gru1«11a  (Fnnmi  12) 

Patricia.  OMQa,  Eraac.  Nobjr  Dick,  Canon  (Water) 
Dadlf  (Water/Ar9on) 

Wtesuranants  (Two-Phasa) 

fHcrotheraocoupla 
Optical  Probas 

Hot  Wire  or  Hot  Flla  Anonoaotry 
Electrical  Probos 

Turbina  FlouMtars  (Fluid  Bearings) 

Venturi  ->  Orifice  Neters 
Heutron  Bean  Attenuation 
NuUlbeaa  Gin  Denslteaeters 
Quick  Clo'jing  Valves 


The  tuo-phase  flou  test  loop  In  Idaho  1?  an  a1r/«ater  loop. 

Again.  Inlet  air  flou.  Inlet  uater  flou.  tanperature.  and  pressure 
awasureaents  are  aade.  Neasureaents  aide  In  the  tuo>phase  test  ;^ect1on 
of  the  loop  are  turbine  aater.  drag  body,  both  discs  and  screens,  aultl 
beaa  gaaaa  denslteaeters.  high  speed  photographs,  and  differential 
pressure. 

Also  located  at  IHEL  Is  a  transient  tuo-phase  flou  loop,  featuring 
a  quick  opening  valve,  and  a  partial  bloudoun  capability  for  rapid 
testing  turnaround  tiae.  Neasureaents  aade  In  the  transient  tuo-phase 
loop  are  aultlbeaa  gaaaa  denslteaeters.  single-beaa  gaaaa  denslteaeters 
quick  closing  valves,  turbines,  drag  body,  pressure,  and  teaperature. 

The  calibration  of  all  mass  flou  Instruaentatlon  for  Sealscale 
occurs  In  ambient  uater.  All  tee  Instnawnts.  turbine  neters,  and  drag 
discs  are  calibrated  In  place.  A  USNBS  traceable  turbine  aeter  serves 
as  a  standard.  The  alr/uoter  loop  previously  aentloned  Is  utilized  to 
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HOLE  ¥ 

TW-PHASE  FLOM  TEST  LOOP.  HCL.  10(00 


Ste«<|y  State  Two-Phase  Flow 

Alr/Uater 
Pressure  60"F 
Teapereture  130  psi 
Mess  Flow  22  kg/sec 

Measurecnewts  (Slngle-Phese) 

Inlet  Air  Flow 
Inlet  Hater  Flow 
Tenperature 
Pressure 

Measurewents  (Two-Phase) 

Turbine 
Drag  Body 
Oise 
Screens 

Nultlbeam  y  Densitometers 
High  Speed  Photographs 
Differential  Pressure 


study  the  relative  accuracy  of  water  calibrations  when  applied  to  two> 
phase  steam/water  environments  such  as  occurs  during  the  Sealscale 
blowdowns. 


Steam/water  two-phase  flow  testing  of  the  Semiscale  spool  pieces  In 
the  Karlsruhe  two-phase  flow  water/steam  loop  Is  being  planned.  Presently, 
much  the  same  approach  Is  taken  for  the  calibration  of  LOFT  drag-disc 
turbine.  Calibration  Is  made  against  the  standard  flowmeter  In  ambient 
water. 


The  Semiscale  calibration  In  aii4)1ent  water  has  proven  to  give 
reasonable  results.  There  Is,  however,  some  criticism  of  the  Semiscale 
experiment  In  that  the  piping  sizes  are  not  typical  -  that  the  size 
almost  guarantees  a  homogeneous  two-phase  flow.  Perhaps  we  have  “lucked 
out"  on  Semiscale  because  of  this.  If  one  does  a  mass  balance  using 


TMIE  VI 

TMHSICNT  TMO-PfUSC  FLOM.  INCL,  IQAMO 


Trawtiwt  Two«Ft>i»«  Flow 


Quick  Opening  Vtivt 
P«rtl«1  BloHdOHR  Upcblllty 

Weesuru— ntt  (Transient  Two-Ph>t>> 

HuUlbtM  Y  OcnsltoMUr 
Single  Bmh  v  OcncltoMter 
Quick  Closing  Vtivts 
Turbin* 

Drag  Body 

Pressure 

Teopereture 


Sealsceic  date,  the  results  are  fnund  vlthln  reason.  Two^phase  mss 
flOM  aeasurcaents  appear  to  have  been  Mde. 

This  Is  good  for  such  a  conpUcated  Integrate  type  experleent.  Ue 
are  now  obtaining  mss  flow  data  froa  LOFT.  The  flow  Is  definitely  not 
hoaogeneous.  Ue  know  from  our  nuttlbean  densftOMters  that  It  Is 
stratified.  Ue  currently  have  no  Mthod  of  Maturing  phase  velocities, 
slip  velocities,  or  density  distribution. 


4.  UH4T  SHOULD  BE  DOHE 


These  MasurcMnts  are  difficult  even  under  Ideal  laboratory 
conditions.  How  can  such  MasureMnts  be  perforMd  In  a  LOFT?  Those 
InstnMnts  which  can  be  ruggedized  to  take  the  environment  must  be 
selected  and  then  It  must  be  figured  out  how  to  Interpret  what  they  say. 
Continued  bench  type  laboratory  testing  will  be  required  If  personnel 
are  to  get  a  handle  on  the  subtleties  of  the  two-phase  flow  measurement 
problems.  Such  testing  requires  that  certain  of  these  unique  Mthods  be 
adopted  and  agree,  to  utilize  them  so  that  data  My  be  compared  from  one 
laboratory  to  another. 


It  wst  to  Ittmto  ton  to  Mtofcitto  itoMttolt  cotoUiMt  -  Mt 
Mljf  fra  tM  toy  to  tto  Mat  witiito  •  Utoratory,  tot  fra 
Ittorttory  to  Utortttvy.  A  m1v«rMl1y  KCtotobU  ttontortf  or  «to  of 
otootorto  ast  to  otttolittod.  Thit  itoltn  tto  tolllty  to  oMtin 
ropottoblo  condUiem.  To  to  tlilSt  ptoM  volocity*  flw  rogla.  ond 
dontity  dlstrftotlm  ast  to  kooa. 

CufTOot  toclmolegy  toai  to  fndlcoto  ttot  tto  host  tots  for  ptoso 
voloelty  aituraif  ts  art  trocar  tactolquot  or  local  protet.  Tto  ast 
attractivo  of  ttoto  aro  tto  trocar  tactolouos  tocauM  ttoy  do  mt 
disturb  tto  onviroamt. 

For  fla  ragfa*  altfboa  gaaa  donsftoatars  aro  swggostad  sfnca 
noarly  all  laboratorlos  havo  soa  oxpartonco  with  tha.  but  local  protes 
could  bo  usod.  Again  tto  donsltoatof .  to  mt  porturb  tto  onvlron, 
Instoad  thoy  aasuro. 

Furttoraro*  ttoso  saa  dovfcos  could  prorldo  tto  donslty  distribu¬ 
tion.  It  can  bo  sold  ttot  fim  regia  and  density  distribution  Im^y 
tto  saa  thing  and  those  tuo  can  bo  cotolned. 

To  suairizo*  calibration  has  been  defined*  although  the  specific 
question  of  utot  a  noan  by  a  t«o-ptoso  calibration  ay  mt  have  been 
rnsarod.  A  look  has  boon  taken  at  what  techniques  are  available  to  us 
and  It  has  been  concluded  ttot  a  standard  Is  needed.  Right  now,  om 
does  mt  exist.  The  proi-las  aro  largo  and  complicated.  Universal  two- 
phase  flow  standards  mst  bo  adopted. 
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Tto  In  dtvtloplnf  •  full  flw  druf  fcrunn  Is  to  obtain  a 

total  —ntii  flui  MaMiroMot  ubicb  uban  coibilnad  with  a  sultablo 
Indopondont  voloelty  or  donsity  ■aasuraaont  will  ylold  a  total  nats 
flux.  Tho  oajor  dotign  consldaratlona  art  prodicatad  by  tbo  fact  that  an 
accurato  ■ouantun  flux  uoasurMant  oust  bo  nado  owtr  a  wfdo  rango  of 
flou  cenditlont.  Tho  dovico  should  oxhibit  a  constant  calibration 
rtgardloss  of  Roynolds  nuubor*  void  fraction,  slip  ratio,  or  flou 
roglua.  Tho  dynaulcs  of  drag  davicos  art  util  undarstood  In  singlo* 
phastflous.  This  Is  not  truo  for  tixHphaso  flous.  Tho  prosant  dovalop- 
■Mit  prograu  Is  diroctad  toward  gaining  an  undorstandlng  of  tho  dynaulcs 
of  drag  davicos  which  saopit  tho  total  arta  of  a  plpo  In  tuo>phast  flow 
and  dovtioping  a  aathod  for  daducing  nass  flow  rata  using  such  a  dovico. 

Various  goonotric  arranganants  ara  to  ba  Invastigatad.  Tasting  to 
data  has  shown  excallant  rasults  using  a  round  wira  wash  scraan  In  tha 
Saulscala  alr/watar  loop.  Future  air/water  tasting  will  Include  per¬ 
forated  Platts  and  wire  washes  with  both  rectangular  and  dianond  shaped 
cross  sections. 

Analytical  uodals  of  the  hydrodynanics  of  tha  drag  screen  as  well 
as  the  associated  density  or  velocity  ueasuring  device  are  being  used  to 
select  tha  optlnua  configuration.  Alternate  force  sensing  uethods  are 
also  being  considered.  These  Include  single  and  nultiple  transducer 
arrangeuents.  Multistage  springs  and  pressure  drop  across  the  body  are 
be  evaluated  for  extending  the  dynaulc  range  of  the  drag  body. 
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A  full  flow  drag  scrotn  Is  boing  dtvtiopod  to  obtain  total  wowantii 
flux  ■Msurwonts  In  Samlscalt*  Lo$s-of-F1u1d  Tost  (LOFT},  and  Battls 
Flask.  Naasurowont  roqulrawants  appaar  In  Tabla  1.  Tfia  davica  Itsalf 
consists  of  a  scraan  or  parforatad  Plata  nhlch  sanplas  tha  antira  cross* 
sactlon  of  a  given  flow  conduit  and  a  aachanlcal  arrangoMnt  to  waasure 
tha  force  on  tha  screen.  A  pressure  d^  ■aasurawent  across  tha  scraan 
Is  Mda  at  tha  saaa  tiaa. 

Tha  prototype  Sawlscala  unit  utilizes  a  three  point  nount  and  three 
Independent  force  sensors.  Figure  1.  This  davica  has  bean  used  to 
evaluate  tha  parforaanca  of  a  round  wire  scraan.  and  Is  currently  being 
used  to  test  rigid  perforated  plates  supplied  by  tha  LOFT  davalopawnt 
prograai. 

Tha  first  prototype  unit  built  by  LOFT  will  use  a  single  force 
sensing  device.  Figure  2.  The  entire  sensor  assembly  will  be  contained 
within  the  pressure  boundary.  An  eddy  current  displacement  sensor  will 
be  used  to  obtain  a  measurement  of  spring  element  displacement. 


2.  OBJECTIVE 


The  objective  Is  to  obtain  a  better  mass  flow  measurement  by 

suitable  combination  of  two*phase  momentum  flux  measurement  with  an 

Independent  density  or  velocity  measurement.  For  homogeneous  flow 

without  slip  an  average  density  measurement  p  from  a  gamma  densitometer 

2 

may  be  simply  combined  with  pV  to  yield  a  mass  flow  rate  pV.  For  flows 
with  slip  where  the  momentum  flux  of  each  phase  Is  of  the  same  order  of 
magnitude  the  combination  Is  not  so  straightforward.  This  problem  will 
be  addressed  both  experimentally  and  analytically  In  the  present  develop¬ 
ment  program.  Suitable  Independent  velocity  measurement  could  be  made 

iBifi  Itek 
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Fl9*  2  Proposed  LOFT  drag  screen  design  configuration. 

by  e  turbine  or  pitot  tube.  Present  date  Indicate  that  aorc  detailed 
velocity  field  neasurenents  My  be  needed  to  properly  Interpret  drag 
screen  data  and  deduce  a  mss  flow  MasureMnt. 


3.  PRESENT  hWmm  TESTIHS 


To  date  the  only  data  which  have  been  analyzed  In  detail  are  round 
wire  screen  data  In  the  Sealscale  alr/weter  loop.  The  coablnatlon  nf 
aoaentua  flux  values  with  average  density  froa  a  two>beaa  densItOMter 
Is  seen.  Figure  3  yields  an  adequate  mss  flow  aeasureaent  for  the 
range  tested.  Other  types  of  screen  eleoents  to  be  tested  are  perforated 
plates  with  and  without  beveled  edges. 


4.  FUTURE  AIR/WTER  TESTING 


First  generation  LOFT  prototypes  will  be  fabricated  in  3-in. 
schedule  160  pipe.  Phase  I  development  testing  will  be  performed  in  the 
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Mim  now  flaw  inpwi  wtt-A-HM* 

Fig.  3  Semiscale  drag  screen  air/water  data  (0.02-in.  diameter  round 
mesh  on  0.2-1n.  centers). 

Semlscale  alr/nater  loop.  The  objective  of  these  tests  Is  to  check 
mechanical  function  and  to  develop  screen  elements  which  provide  a 
constant  drag  coefficient  over  a  wide  range  of  flow  conditions.  Figures 
4.  5.  and  6.  Parameters  which  should  be  considered  for  screen  design 
are  Insensitivity  of  drag  coefficient  to  Reynolds  number,  flow  regime, 
slip  ratio,  and  void  fraction.  The  flow  restriction  should  be  minimized 
as  much  as  possible  and  still  make  an  accurate  total  momentum  flux 
measurement. 
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S  iir  JD  40MM10*  a  aMMto*  a  a  aai9  a  aaaia  a  aaaio' 

Wawem  Htimtm  R»  •  ^dir/|| 

Fig.  4  Reynolds  nuabtr  dependence  In  s1ng1e>phese  flow  (drag  coefficient 
for  a  circular  cylinder  and  thin  disc  and  pressure  drop  coefficient  for 
a  nuUlholed  plate). 

5.  FUTURE  $TE/W/HATER  TESTING 


Phase  11  development  testing  will  be  transient  flows  of  steam/water 
performed  In  the  Bettis  Flask  Blowdown  Facility.  The  objective  of  these 
tests  is  to  verify  performance  in  a  transient  steam/water  flow. 


6.  CONCLUSIONS 


Preliminary  data  have  shown  that  a  full  flow  drag  screen  can  make  an 
accurate  measurement  of  total  momentum  flux  over  a  wide  range  of  flow 
conditions.  For  certain  ranges  these  data  can  be  combined  with  average 
density  from  a  gamma  densitometer  to  yield  an  accurate  mass  flow  rate. 
More  development  testing  is  required  to  be  confident  that  a  good  mass 
flow  measurement  can  be  obtained  over  the  wide  range  of  conditions 


present  in  a  LOFT  blowdown. 


Fig.  5  Alr/wter  te»t  section  insert. 


Fig.  6  3-in.  diameter  drag  screens 
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OMfi  SCREEN  DESIGN  PARMCTER  CONSIDERATIONS 


APPtMOlX  A 


DRAG  SCREEN  DESIGN  RARAHHER  CONSIDERATIONS 


Soat  blunt  bodies  ere  known  to  exhibit  constent  dreg  coefficients 
for  Reynolds  nuabers  ebove  e  certein  alnlaua  velue.  for  exenple*  thin 
circuler  or  squere  ptetes  for  Reynolds  nuabers  ere  besed  on  height  or 
dleaeter  ebove  ebout  3000^^^.  For  e  given  Reynolds  nuaber  held  constent. 
end  the  length  to  height  retio  of  e  rectenguler  section  verled.  the  dreg 
coefficient  becoaes  e  constent  for  e  length  to  height  retio  of 
epproxiaetely  1/3^^^.  Thus,  for  e  dreg  body  In  en  Infinite  flow  field 
(e  saell  body  In  e  Urge  pipe)  to  beheve  like  e  thin  p1e*e  or  disc.  Its 
thickness  should  be  less  then  1/3  Its  height. 

For  aultiple  bodies  such  es  screens  or  aultiholed  pletes  this 
criterion  Is  probably  unnecessarily  restrictive  for  the  height  to  length 
ratio.  The  plenomenon  which  causes  the  drag  coefficient  for  a  thin  disc 
or  plete  to  becoae  constent  for  e  Urge  enough  Reynolds  number  Is  the 
behavior  of  the  flow  field  around  the  bod/.  The  flow  Is  unable  to 
negotiate  the  sharp  corner  end  separates  from  the  body.  For  Reynolds 
nunbers  above  3  x  10^.  the  flow  field  changes  very  little  end  the  drag 
coefficient  Is  essentially  constant.  If  the  body  1$  made  thicker,  the 
possibility  of  the  separated  fl^w  reattaching  Itself  to  the  body  must  be 
considered.  If  this  occurs,  the  drag  coefficient  will  change,  possibly 
red;red  by  es  auch  as  a  factor  of  2^^^. 

It  Is  well  known  that  a  pressure  rise  Is  required  for  reettechnent 
of  a  separated  shear  layer  to  a  body^^*^*^^.  This  pressure  rise  Is 
achieved  through  the  transport  of  laomentua  or  mixing  (either  laminar  or 
turbulent)  which  accelerates  the  fluid  In  the  recirculation  zone  and 
results  In  a  net  slowing  or  diffusing  of  the  fluid  in  the  external 
flow^^^.  This  slowing  of  the  external  flow  creates  a  pressure  rise. 

The  acceleration  of  the  slower  moving  fluid  In  the  recirculation  zone 

bwiriini  mil  Hal 
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IncrMstt  Us  klnsttc  energy  whicli  eventuelly  enables  U  to  overceae  tbe 
adverse  pressure  gradient  and  leads  to  reattactaent  as  sboun  In  Figure 
A-1. 


Fig.  A-1  Sketch  of  the  recirculation  zone. 


As  the  Reynolds  number  Increases,  a  turbulent  shear  layer  attains  a 
limiting  form  and  thus  becomes  Independent  of  Reynolds  number.  This 
effect  Is  ezsily  seen  In  case  of  constant  drag  coefficient  for  a  thin 
disc  or  plate  for  Reynolds  numbers  above  3000.  The  same  effect  Is 
apparent  for  many  other  types  of  turbulent  flows,  such  as  the  discharge 
coefficient  of  single  hole  orifice  plates^^^  for  large  Reynolds  nianbers. 
For  a  separated  shear  layer  on  a  semi Infinite  blunt  flat  plate^^^  the 
length  to  height  ratio  required  for  reattachment  Is  on  the  order  of  5. 

In  the  Reynolds  number  range  studied  In  Reference  2.  Re  >  2  -  7  x  10^. 
the  shear  layer  is  turbulent  and  the  length  required  for  reattachment 
Is  constant. 

Consider  a  screen  with  square  mesh  or  muHIholed  orifice  plate  with 
square  corners.  If  the  flow  blockage  is  significant,  greater  than  about 
lOS,  there  will  be  a  pressure  drop  across  the  device  due  to  the  acceler¬ 
ation  of  the  fluid  caused  by  a  reduction  In  area  available  to  the  flow. 

Since  a  negative  pressure  gradient  is  Imposed  on  the  separated 
shear  layer  by  the  accelerating  external  flow  the  pressure  rise 
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•ttocUUd  witli  rMtUdMnt  it  ntvtr  ttuirnd  rtgardlm  of  tlw 
Rtyiio14t  nmb&r.  Tlit  4r«t  cMfflcItnU  of  tko  ceoftsurottoos  dtscrtNd 
•bovt  should  thus  ho  consUot  for  ony  Roynolds  ouobor  obovt  about  ICp. 

For  vory  thick  otvlcts.  tho  flow  choimtl  could  concoivtbly  occoloroto 
tho  flow,  thou  «ct  os  a  diffusar  and  rocovar  prossuro.  A  rtviow  of 
subsonic  conical  diffustr  data^^^  shows  this  to  bo  vary  unllkaly  uniass 
tha  Plata  is  on  tha  ordar  of  savaral  scraan  or  orlfica  plata  hola  diaaetars 
thick.  Scraan  alaaants  aada  froa  blunt  bodies  should  not  suffer  froa 
Reynolds  nuabar  affects. 
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SUMIMV 


CvalMtIon  of  data  froa  Loss«of-F1u1d  Teat  (LOFT)  nonnuclear  tests 
has  shown  a  need  for  profile  neasureaents  at  several  locations  In  the 
LOFT  piping.  A  prototype  rake  consisting  of  three  Drag>01sc  Turbine 
Transducers  (OTT)  has  been  designed  and  fabricated  for  Installation  at 
one  location  (FE-BL-1)  In  the  blowdown  loop.  After  successful  operation 
during  a  LOFT  nonnuclear  test  (L1>4)  scheduled  for  Kay,  1977,  additional 
rakes  will  be  Installed  In  the  primary  and  blowdown  loops.  A  research 
program  to  develop  a  pitot  tube  rake  for  measurement  of  steady  state  and 
transient  two-phase  flows  Is  In  progress  at  HcMaster  University,  Hamilton 
Ontario.  A  rake  of  thermocouples  and  pitot  tubes  will  be  developed  for 
Installation  near  the  emergency  core  coolant  (ECC)  Injection  points. 
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1.  IHTRODUCTION 


Evaluation  of  data  from  Loss-of*F1u1d  Test  (LOFl)  nonnuclear  tests 
has  shown  r  for  profile  neasurenents  at  several  locations  In  the 
LOFT  piping.  Instrument  arrays  (rakes)  have  been  designed  and  fabr1> 
cated  to  both  Improve  the  accuracy  and  provide  redundancy  for  measure* 
ments.  This  presentation  will  summarize  the  status  of  the  LOFT  rake 
Implementation  program. 


2.  MASS  BALANCES  USING  LOFT  INSTROMEKTS 


Figure  1  shows  a  view  of  a  Drag-Disc  Turbine  Transducer  (DTT) 
manufactured  for  Installation  In  the  14-1n.  O.D.  schedule  160  LOFT 
piping.  Note  that  the  Instrument  Is  positioned  such  that  a  single  free- 
fleld  measurement  Is  made  at  the  center  of  the  pipe.  A  three-beam  gamma 
densitometer  to  measure  average  chordal  density  Is  strapped  around  the 
pipe  near  the  DTT  location.  Figures  2  and  3  show  DTT  Instrument 
locations  on  the  LOFT  primary  and  blowdown  loops.  Three  DTT  locations 
are  Identified  In  the  primary  loop  and  two  are  located  In  the  blowdown 
loop.  Figure  4  shows  mass  calculations  made  using  data  taken  during  the 
LOFT  LI -2  nonnuclear  test  with  the  combinations  of  Instruments  as 
Identified.  The  solid  horizontal  line  represents  the  total  mass  of 
fluid  In  the  system.  The  dashed  horizontal  line  represents  the  mass 
actually  measured  by  the  various  Instruments.  The  mass  of  fluid 
Initially  In  the  piping  downstream  of  the  DTT  locations  has  been  sub¬ 
tracted.  Note  that  a  mass  balance  resulting  from  Integration  of  drag- 
disc  turbine  measurements  (curve  tagged  with  small  squares)  gives  a 
result  which  Is  almost  twice  the  correct  value.  During  a  LOFT  loss-of- 
coolant  experiment  (LOCE)  test,  stratified  flow  occurs  and  a  single  free 
field  measurement  Is  Inadequate.  Rakes  of  three  DTTs  will  be  Installed 
at  each  of  the  piping  locations  where  single  units  are  now  Installed  to 
Improve  accuracy  and  Increase  redundancy.  The  curves  also  show  a  mass 
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Fig.  2  LOFT  operating  loop  with  instrument  locations. 
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(nstniMiit  iMttlom 


btlMCt  using  calcutstnd  asst  fiM  bsstd  on  donslty  ■otsurod  vis  tlw 
ganw  dansltoaotor  and  a  dlfforontlsl  prtssurt  (aMSurod  st  •  piping 
location  ontitlod  "contraction  a«"  on  Figure  3).  The  calculated  asss 
balance  using  &o  Is  In  good  agreeaent  with  the  actual  Initial  value. 
Pitot  tube  arrays  will  be  designed,  fabricated,  and  Installed  In  the 
LOFT  piping  to  facilitate  calculation  of  flow  profiles. 


3.  RAKE  REQUIREWPfTS 

The  LOFT  facility  Is  a  scaled  design  of  a  Westlnghouse  type 
pressurized  water  reactor  (PWR).  Each  of  the  rake  designs  will  require 
penetrations  of  the  prlaat  *  coolant  pressure  boundaries.  The  saae  codes 
and  standards  required  for  PIffi  reactors  will  be  followed  during  the 
design,  fabrication  qualification,  and  Installation  of  the  rakes  on  the 
LOFT  systea. 

3.1  DTT  Rakes 

A  fonnal  requirenents  document  wa^  r spared  and  approved  after 

design  reviews.  Table  I  summarizes  ;  .  Uonmental  conditions  and 

required  ranges.  The  maximum  momentum  rlux  value  for  design  purposes 

2 

was  188.000  lbm/ft>sec  as  determined  from  RELAP  computer  models  and 
previous  nonnuclear  tests. 

3.2  ECC  Rakes 

One  of  the  goals  of  the  LOFT  LOCE  testing  program  is  to  define 
emergency  core  cooling  (ECC)  mixing  phenomena.  Temperature  measurement 
arrays  on  either  side  of  the  ECC  injection  points  would  give  indications 
of  ECC  oscillations  about  the  Injection  port  and  the  rate  at  which  the 
water  temperature  increases.  A  similar  array  on  the  broken  loop  will 
Indicate  the  condition  of  ECC  bypass.  Installation  on  both  hot  and  cold 
leg  sides  Is  required  for  potential  use  of  hot  leg  ECC  Injection. 


TAiLt  I 

on  MKi  REQUIREMENTS 


I 


1.  Design  to  Meet  ROT  Stenderds  for  LOFT  Conditions: 
Tampcreture  300>6S0*F  (tSO-34S*C) 

Pressure  22S0  psig  (1S5  BAR) 

KaxInuM  Velocity  100  ft/sec  (30.5  a/sec) 

Redletlon 

Neutrons  > 

Nonnel  Operetlons  1  x  10°  nv 

Total  Exposure  7  x  10^^  nvt 


Games 

NohmI  Operation  1000  R/hr 

a 

Total  Exposure  7  x  10  R 


Ranges: 

(Ib/ft-sec^) 

(Kg/m-sec^y 

11 

Top 

2,000-50,000 

2,976-74,400 

Middle 

2,000-50,000 

2,976-74,400 

Bottom 

2,000-50,000 

2,976-74,400 

#2 

Top 

200-3,500 

298-5.209 

Middle 

800-14,000 

1,191-20,835 

Bottom 

200-3, sr 

1,191-5,209 

In  addition  to  the  temperature  measurement,  it  Is  possible  to  Include 
other  measurements  on  the  probe  to  further  Improve  flow  measurement 
accuracy.  INEL  Is  currently  sponsoring  a  university  Investigation  of 
the  usefulness  of  pitot  tubes  fo^  two>phase  flow  measurements.  Pitot 
tubes  will  also  provide  a  good  measurement  of  the  Initial  flow  profile 
In  the  pipes.  The  Initial  CCC  rake  design  will  thus  Incorporate  pitot 
tubes  and  thermocouples. 

Summarizing  the  ECC  rake  design  requirements  Include  the  following: 

Range: 

Pitot  Tubes  -  20  -  200,000  Ib/ft-sec^  (30-207,600  Kg/m-sec^) 

T/CS  -  100  -  1300”F  {38-722'’C) 

Response  100  msec 
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InsUllatlon:  PC-1  (A  «  B).  PC-2  (A.B)  BL-1.  BL-2 

Location  of  Probts: 

PC-1  (ABB).  PC-2  (ABB) 

10  Pitot  Tubes.  5  T/CS  Each 
BL-1.  BL-2 

S  Pitot  Tubes.  5  T/CS  Each 

Naximun  Forces  During  BIOMcoMn: 

180.000  Ibs/ft-sec^  (279.800  Kg/m-sec^) 

A  total  of  six  ECC  rakes  Mill  be  fabricated  and  qualified  Mith  tm  each 
Installed  at  locations  PC-1  and  PC-2  on  the  primary  coolant  loop  nnd  Mf. 
each  Installed  at  locations  BC-1  and  BC-2  on  the  broken  loop.  McNaster 
University.  Hamilton.  Ontario  (Dr.  S.  Banerjee.  principal  investigator), 
has  received  a  FY  77  subcontract  to  Investigate  the  performance  of  pitot 
tubes  In  tMO-phase  flow.  The  results  of  this  research  will  be  Incorp'*- 
rated  Into  the  ECC  rake  design. 


4.  PITOT  TUBE  RESEARCH  AT  MCMASTER  UNIVERSITY 

Or.  S.  Banerjee  was  In  attendance  at  the  Nuclea*'  Regulatory 
Coamlsslon  (NRC)  conference.  He  mentioned  that  the  progrc.i  would 
Include  evaluation  of  pitot  tubes  In  both  steady  state  and  transient 
two-phase  phenomena.  The  steady  state  data  have  been  taken  uslr*  an 
array  of  five  pitot  tube  probes  In  three  different  pipe  diameters  - 
1-1/2.  2-  end  3-1n.  0.0.  The  data  from  the  tests  are  still  being 
evaluated.  Or.  Banerjee  did  show  Figure  5  and  pointed  uc  that  the 
calibrations  made  with  all  water  (void  fraction  zero)  seemed  to  also 
hold  for  the  higher  void  fraction  test  data.  The  transient  steam/water 
tests  will  be  run  this  spring  and  suinner.  A  summary  report  will  be 
available  by  the  end  of  FY  77. 
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rakt  will  bt  Installed  at  location  Bl-1  for  tht  next  LOFT  nonnuclear 
test,  test  LI-4  scheduled  for  Nay  1977.  The  second  unit  wHi  be  used 
for  additional  tests.  After  evaluation  of  the  LI-4  test  data,  the  rake 
design  will  be  modified  as  required  and  the  remaining  DTT  rakes  fabricated 
and  Installed  at  locations  BL-2,  PC-1.  PC-2,  and  PC-3. 

5.2  ECC  Rake 

The  requirements  for  the  ECC  rake  are  still  being  developed. 

Several  conceptual  designs  have  been  considered.  One  of  the  conceptual 
designs  Is  shown  In  Figure  7.  Fabrication  of  the  ECC  rakes  will  be 
completed  In  time  for  installation  prior  to  the  L1-5  LOFT  nonnuclear 
test. 
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SUWMRV 


Tht  prtscnt  Drag-Disc  Turbin*  Transducers  (DTT)  Installed  In  the 
Loss-of-Fluld  Test  (LDTr)  facility  were  analyzed  and  redesign  goals  were 
established  during  FY  76.  A  prototype  DTT  was  developed  anc:  single¬ 
phase  and  separate  effects  testing  was  completed.  An  alternate  eddy 
current  sensor  for  the  DTT  drag  body  was  developed  and  tested  during 
FY  76.  The  results  of  the  separate  effects  testing  progran  will  be 
Incorporated  In  the  design  and  development  of  a  second  prototype  which 
will  be  tested  under  single-  and  two-phase  conditions  during  FY  77. 
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I.  NaHiig  L1f«  Test  Conditions . 

II.  Turtlno  lotrlngs  Lifd  Tost  RosuUs  .  . 

III.  «W  (Or.  R.  Lohty)  Objoctlvos  of  Progri* 


1.  IliTIMOUCTIOII 


Tftt  d*s1gn  and  perfonnnce  of  the  Oreg>01sc  Turbine  Transducer 
(OTT)  currently  Installed  at  various  locations  In  the  Loss>of-F1u1d*Test 
(LOFT)  facility  Mere  analyzed  and  redesign  goals  Mere  established  during 
TY  76.  The  redesign  goals  Included:  (1)  retain  as  aany  of  the  preseiitly 
used  CQMponents  (sensors,  turbine  blade,  etc.)  as  feasible;  (2)  nake 
separable/ replaceable  turbine  and  drag-disc  transducers;  (3)  extend 
ranges;  (4)  strcanllne  Mhere  possible;  (5)  leprove  turbine  bearing 
llfetlae;  (6)  reduce  friction  In  the  drag-disc  sensor  linkage; 

(7)  define  perfomance  analytically  and  perforaance  of  adequate  calibra¬ 
tions.  A  tMO  year  program  mss  Initiated  to  Improve  the  DTT.  The  design 
and  test  program  stressed  component  Improveamnts  -  particularly  modularity 
of  units.  Improvement  of  turbine  bearings,  and  development  of  an  alter¬ 
nate  drag  body  sensor.  The  results  of  the  separate  effects  testing 
program  Mill  be  summarized.  The  FY  77  goals  for  development  and  testing 
of  a  second  prototype  OTT  Mill  be  presented. 


2.  FIRST  PROTOnPE  OH 


The  OTT  .urrently  qualified  for  Installation  Is  shOMn  In  Figure  1. 
The  single  body  Instrument  Incorporates  a  six  bladed  turbine  mounted  on 
a  17-4  PH  steel  shaft  contained  by  gold  alloy  sleeve  bearings  in  each  of 
tuo  support  posts.  The  rotation  of  the  turbine  blades  Is  sensed  via  an 
eddy  cun'-nt  sensor  located  In  the  body  beloM  the  blades.  The  other  end 
of  the  OTT  Instrument  has  a  circular  disc  drag  body  attached  to  lever 
arms  and  a  pivot  to  position  a  variable  reluctance  transformer  core 
as  the  fluid  momentum  flux  on  the  disc  acts  against  a  torsion  bar 
restoring  spring.  The  prototype  I  design  Is  shomn  In  Figure  2.  The 
sensors  from  the  present  LOFT  TIT.  Mhich  had  previously  been  qualified 
for  LOFT  service,  mere  used  In  the  prototype. 
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Present  LOFT  DTT  with  shroud  removed. 


Fig.  2  Nodular  prototype  DTT. 


2.1  Modular  Design 

As  seen  in  Figure  3,  the  turbine  and  drag  body  are  constructed  as 
Identical  sized  separate  units  with  interlocking  key-ways  in  order  that 
single  or  multiple  units  can  be  joined  together  as  required.  Results  of 
fabrication  and  testing  of  three  copies  each  of  the  turbine  and  drag 
body  units  showed  that  the  manufacturing  costs  could  be  reduced.  Stress 
analyses  shotted  that  the  shroud  and  other  modular  parts  would  withstand 
the  maximum  conditions  encountered  during  LOFT  LOCE. 

2.2  Performance  Tests 

In  the  single  phase  tests  of  the  prototype,  the  LOFT  environmental 
conditions  of  2250  psig  and  SSO^’F  were  simulated  by  the  use  of  a  large 
autoclave  chamber.  Autoclave  tests  on  each  of  the  three  copies  of  the 
prototype  DTT  indicated  that  the  combined  effects  due  to  changes  in 
temperature  and  pressure  from  ambient  to  LOFT  conditions  caused  less 
than  IX  (of  full  scale)  error.  Test  results  showed  that  the  turbine 
unit  had  a  useful  range  of  5  to  40  ft/sec.  The  drag-disc  pivot  spring 


afttr  eyelM  of  full  tctit  Mfluctlon.  Thu  t^lng  will  to 
rtoMlflfito  to  MOt  •  fallurt  lift  of  10*  cyclot. 

3.  TURBINE  BEARIto  ItoBOVDttNT  POOGtoM 


The  service  life  goal  for  turbine  flowmeters  Installed  In  the  LOFT 
system  Is  10*000  hours  of  operation  at  LOFT  conditions.  The  Improvement 
program  Includes  life  testing  on  bearing  materials*  testing  of  alternate 
turbine  bearing  designs*  and  development  of  compound  bearing  concepts. 

3.1  Searing  Life  Tests 

Table  I  shows  the  test  conditions  of  the  test  loop  In  which  various 
combinations  of  bearing  materials  have  been  tested.  A  sumnary  of  the 
results  to  date  Is  shown  In  Table  II. 

TABLE  I 

BEARING  LIFE  TEST  CONDITIOKS 

Flow  Velocity:  13  ft/sec  (4  m/sec) 

Temperature:  600-620*F  (Approximately  320"C) 

Operating  Pressure:  2200  psi  (153  BAR) 

Demineralized  Water  with  Low  Chlorides  and  Oxygen  Content 


3.2  Single  Post/Bearinq 

A  single  post  turbine  has  been  developed  as  shown  In  Figure  4* 
for  evaluation  of  the  design*  bearing  materials  comblnstlons*  and  the 
effects  of  forced  fluid  lubrication.  The  combinations  of  materials  to 
be  tested  Include  sleeve  bearings  of  (1)  H/H  gold  alloy*  (2)  P*658-RCH 
graphite  and  (3)  440C  stainless  steel  each  combined  with  a  one  piece 
440C  stainless  steel  shaft  and  collar.  The  same  combinations  will 
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be  tested  In  a  single  post/bearing  with  water  lubrication.  A  440C 
sleeve  and  shaft  lopregnated  with  diaaond  powder  will  be  tested. 

3.3  Coapound  Bearing 

Or.  Jack  Cole*  University  of  Arkansas*  has  a  research  cuntract  for 
FT  77  to  develop/test  a  coapound  bearing  concept.  Figure  5  shows  a 
coapound  bearing  which  Is  currently  under  developaent  and  testing  by 
the  Mechanical  Engineering  Departaent  of  the  University  of  Arkansas 
under  a  FT  77  subcontract  (Or.  Jack  Cole*  principal  Investigator). 

The  bearing  concept  consists  of  a  Jewelled  bearing  which  acts  against  a 
spring  loaded  thrust  bearing  until  forces  on  the  turbine  asseably 
exceed  the  spring  constant,  then  the  Jewelled  bearing  seats  against  a 
friction  shoulder.  Sleeve  bearings  act  as  rugged  long-life  bearing 
surfaces  as  long  as  the  high  thrust  conditions  exist.  A  patent 
application  has  been  filed  on  the  concept. 
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single  bearing  teats. 


Fig.  5  Coapound  bearing  turbine  floweeter  concept  for  reducing  friction 
at  low  flows. 

IHPROVED  0WA6  r.0DY  SENSOR 


4.1  Design 

The  present  LOFT  OTT  unit  has  a  Variable  Reluctance  Transfomer 
(VRT)  and  core  as  a  sensor.  Ihe  lever  arns.  pivot  bearings*  and  VRT 
clearances  give  rise  to  fabrication,  calibration,  and  service  life 
probleas.  An  alternate  sensor  has  been  designed,  fabricated,  and 
tested  which  utilizes  eddy  current  sensor  colls  to  aeasure  dlsplace- 
aent  of  the  drag>d1sc  am  without  any  contact  or  connection  required. 

The  asscidtled  eddy  current  sensor  Is  shown  In  Figure  6.  One  end  of  the 
sensor  contains  a  sensing  coll  to  sense  a  reference  target  and  the  other 
end  has  an  Identical  coll  to  aeasure  deflection.  The  colls  are  wound 
with  Identical  geometries  and  electrical  properties.  When  connected  into 
an  ac  bridge  circuit,  the  Identical  colls  are  very  insensitive  to  large 
changes  In  temperature. 


Fig.  6  Assenbled  eddy  current  transducer. 

4.2  Test  Results 

Figure  7  shows  a  plot  of  test  results.  The  sensor  has  a  useful 
displaceattnt  neasurenent  range  of  0.0  to  0.150  In.  The  expected  range  over 
which  the  sensor  will  be  used  In  a  drag-disc  Is  0.060  In.  as  Indicated 
by  the  vertical  narks  Intersecting  the  curves.  Test  results  have  shown 
that  the  effects  of  sudden  temperature  changes  from  689  to  73°F  cause 


Fig.  7  The  eddy  current  transducer's  dc  output  versus  target  dlsplace- 
nent. 

less  than  Sf  (of  full-scale)  error.  The  tests  In  an  autoclave  Indicated 
that  effects  of  pressure  changes  from  ambient  to  2250  psig  caused  less 
than  2X  (of  full-scale  e  ror). 


5.  COWCLUSIONS  AND  FY  77  GOALS 


5.1  FY  77  Goals 

The  FY  77  g(^1s  Include  (1)  reduction  In  size;  (2)  new  turbine 
bearings;  (3)  new  drag-disc  sensor;  (4)  removable  sensors,  and  (5)  two- 
phase  calibrations.  A  final  design  review  to  evaluate  the  prototype  I 
will  be  held  In  narch,  1977.  A  second- prototype  will  be  designed  to  the 
FY  77  goals,  fabricated  and  testing/callbratlons  begun  during  the 
remainder  of  FY  77  will  be  completed. 


1.5-88 


S.2  OwlytlMl  HodtU 

RPI  (Or.  R.  Uhcy  Jr.,  principal  investigator)  Has  accepted  a 
subcontract  to  develop  and  evaluate  analytical  perfonnance  Models  for 
the  LOFT  turbine  and  OTT  units.  Table  III  suanarizes  the  objectives 
of  the  prograH  for  FY  77. 

TABLE  III 

RPI  (DR.  R.  LANEY)  OBJECTIVES  OF  PROGRAM 


1.  Develop  a  detailed  analytical  model  for  the  tre-'sient  response  of  a 
turbine  meter  in  sing1e>phase  and  two*phase  flow.  This  isodel  would 
include  such  things  as  rotor  inertia,  bearing  friction,  etc. 

2.  Synthesize  a  method  of  Interpreting  local  readings  in  terms  of 
global  parameters  of  interest. 

3.  Investigate  the  effect  of  virtual  mass  on  the  response  of  a  drag- 
disc  in  transient  two-phase  flow. 

4.  Suggest  experiments  to  investigate  the  interaction  that  occurs  be¬ 
tween  instruments  in  series  (e.g.  How  does  the  turbine  affect  the 
drag-disc?).  Based  on  the  results  of  these  experiments,  to  provide 
analytical  models  and/or  correlations,  to  account  for  these  inter¬ 
action  mechanisms. 

5.  Model  and/or  Instnment  changes  will  be  proposed  based  on  the 
results  of  the  analytical  and  experimental  work  done  in  items  1*4 
above. 


5.3  Conclusions 


Evaluation  of  data  taken  during  three  nonnuclear  tests  of  the  LOFT 
facility.  Indicated  that  the  presently  installed  DTTs  have  provided 
valuable  and  essential  Information;  however,  problems  in  manufacture, 
accuracy,  and  reliability  were  noted.  Test  results  from  single  effects 
testing  and  from  the  design,  development,  and  testing  of  a  prototype  I 
demonstrated  that  most  of  the  problems  could  be  minimized  or  eliminated. 
Defined  design  directions  for  a  second  prototype  are  to  be  developed 
during  FY  77. 
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IteasurcMnt  requlrcaants  and  Initial  constraints  havt  boon  dafinad 
for  tha  nuclaar  fualad  Loss>of-F1u1d  Test  (LOFT)  dansitoaatar.  The 
principal  probleai  Is  to  maka  tha  photon  transmission  ■easurenent  In  tha 
presence  of  an  expected  heavy  background  radiation  field  due  to  nitrogen* 
16  and  other  activated  fluid  and  corrosion  products  In  the  itatar.  Also, 
further  Information  on  the  mass  distribution  over  the  flow  field  Is 
desired.  Different  Instrumental  techniques  are  being  considered  to 
solve  these  problems. 

He  are  conducting  hardware  tests  at  the  Advanced  Test  Reactor  (ATR) 
to  determine  the  feasibility  of  solving  the  background  radiation  problem 
by  straightforward  shielding  and  digital  background  subtraction  (energy 
discrimination)  techniques.  Error  analyses  of  the  chordal  average 
density  measurement  show  the  advantages  to  be  gained  by  selection  of 
detector  position  and  source  energy. 

The  configuration  presently  being  pursued  consists  of  the  Nal  type 
detector,  single  source,  ano  digital  background  subtraction.  Subsequent 
modifications  may  Include  use  of  Germanium  detectors,  multiple  sources, 
and  perhaps  on-line  analog  background  subtraction. 
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!.  IMTHOaUCTIOi 


OtnsItOMters  ustd  during  the  nonnuclear  LOFT  testing  need  to  be 
Modified  In  order  to  perfone  their  Intended  function  idien  subjected  to 
the  high  background  radiation  field  expected  during  nuclear-fueled  LOFT 
wrk.  The  developMent  work  Involves,  first,  deterplnatlon  of  the  design 
requIrcMents  for  the  electronic  and  xiechanical  equIpMent  and  second, 
detenalnatlon  of  data  reduction  techniques  to  account  for  flow  regine 
and  SMall  L/0  flow  disturbance  effects.  The  work  accoig>11shed  to  date 
Is  In  support  of  the  first  objective. 


2.  MEASUREMENT  REQUIROIEWTS  AMO  TECHNIQUES  BEING  EVALUATED 


MeasurcMent  requirements,  revised  to  reflect  data  gathered  from  the 
first  few  nonnuclear  LOFT  tests,  have  been  established  for  the  advanced 
densitometer.  These  requirements,  along  with  the  principal  Measurement 
constraint,  (background  radiation)  are  given  In  Table  I. 

TABLE  ! 

LOFT  DENSITOMETER  MEASUREMENT  REQUIREMENTS 


Required 

Goal 

Kange,  Ib/ft^ 

0.7  -  62.4 

0.06  -  62.4 

Accuracy.  Ib/ft^ 

♦  0.6 

♦  0.3 

MixImum  Sampling  Time,  sec 

0  1 

0.016 

Number  of  Beams 

at  Least  2 

9 

Background  Radiation 

Source  Strength 

140  b  Cl  -  N-16 
cm^ 

Plus  20  u  Cl 

cm 

Mixed 

Activation, 
Fission  and 
Corrosion 
Products 

nWBHH  PV 
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TIm  prlairy  aettiod  of  addressing  the  hackground  radiation  problaa 
Is  to  separately  Measure  and  subtract  It  froM  the  coMbIncd  signal  plus 
background  aeasuraMent.  Photon  energy  analysis  type  1nstn«Mntat1on  Is 
used  to  acconpllsh  this  background  subtraction.  The  standard  Method 
Involves  off-line  digital  background  subtraction,  an  Instrtaent  block 
dlagran  for  which  Is  ..hown  In  Figure  1.  With  this  equIpMent,  the 
Infornatlon  obtained  and  recorded  Is  only  the  portion  jf  the  total 
photon  energy  spectrun  of  Interest  to  the  densitoaeter  aeasurenent. 

This  portion  of  the  spectrun  consists  of  the  full  energy  (photo)  peak 
plus  twenty  to  forty  channels  on  either  side  of  the  peak  used  to 
establish  the  coefficients  o^  an  equation  for  the  background.  This 
equation  Is  used  to  detenalne  background  values  under  the  peak  which  are 
subtracted  In  order  to  obtain  the  net  peak  area  (above  background).  This  - 
net  peak  area  Is  the  nunber  which  Is  properly  related  to  the  chordal 
average  density.  Once  the  off-line  digital  technique  has  been  Instituted, 
and  Its  accuracy  and  capabilities  established.  It  nay  be  possible  to 
sinplify  the  equlpnent  and  data  processing  procedures  by  going  to  an 
on-line  analog  background  subtraction  technique.  The  block  diagran  for 
thif.  Is  shown  In  Figure  2.  Essentially,  the  background  under  the  peak  Is 
detemlned  by  averaging  the  background  on  either  side.  This  can  be  done 
with  analog  techniques  during  the  conduct  of  the  test,  so  that  only  the 
(approxlnate)  peak  area  Is  actually  recorded.  Figure  3  shows  both  the 
existing  single  source/trlple  bean  type  densiteneter  and  a  proposed 
Inprovesent  to  that,  the  triple  source/nine  bean  densitoneter.  Each 
source  enits  photons  of  different  energies  and  radiates  to  all  three 
detectors.  Thus,  each  detector  sees  all  three  sources  and  three  spectral 
peaks  occur  for  each  detector  channel.  These  energy  peaks  are  associated 
with  the  corresponding  fluid  chordal  path,  so  that  a  total  of  nine 
sinultaneous  neasurenents  are  nade.  (Such  a  setup  would  be  more  con¬ 
sistent  with  the  off-line  digital  background  subtraction  technique.) 
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Flo.  3  Instrumental  tKhnIques  being  considered. 
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To  exporinenuily  establish  the  feasibility  of  the  slaple.  direct 
background  subtraction  technique*  tests  have  been  conducted  at  the 
Advanced  Test  Reactor.  Figure  4  sliows  the  overall  test  arrangeaent  of 
the  source/detector  setup  at  one  of  the  reactor's  prinary  coolant  outlet 
lines.  Figure  5  shows  photographs  of  the  actual  Installation.  Figure  6 
Is  a  photon  energy  spectrum  obtained  from  this  test.  The  main  peak  Is 
due  to  the  densitometer  Cs>137  snurce.  The  Nitrogen-16  background 
radiation  leaking  through  the  shield*  as  well  as  radiation  from  other 
activated  fluid  and  corrosion  products  causes  the  observed  background. 
Since  the  N-16  concentration  Is  larger  than  expected  In  LOFT*  the  volume 
of  water  Is  larger  and  the  pipe  wall  Is  thinner*  the  radiation  source 
strength  Is  considered  to  be  substantially  larger  than  will  be  experi¬ 
enced  In  LOFT.  Thus*  this  result  serves  to  establish  the  feasibility  (at 
least  for  steady  conditions)  of  the  use  of  background  subtraction  In  a 
high  background  radiation  field. 
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Fig.  Sb  ATR  background  radiation  test  (view  of  sources,  pipe,  detector  shield). 


F1|.  6  FtotM 


iptctnw  ftroi  advtimdl  ttst  rttctor  and  Ct-137 


im  error  Mwljrtls  of  tlw  cinrdtl  avorofe  domity  MUoroaMit  m% 
ORdortaloR.  Seat  of  tlio  loltlol  roiolts  aro  them  la  Flgoroa  7  aad  8. 
la  both  figiiros,  a  statistical  orror  factor  Is  plotted  idildi,  obaa 
divided  hF  13.67  gives  the  perceatage  error  la  the  chordal  average 
daasItF*  for  the  LOFT  setup  (large  pipe).  As  shoaa  la  Figure  7,  ‘dm 
■latauB  error  occurs  fbr  a  diaaetral  boaa  path  (0  ■  0*)  aad  lacraeses 
rapidly  beyead  about  16*.  la  Figure  8,  the  error  factor  Is  plottad 
agalast  photoa  easrgy  for  too  aagles  aad  four  daasitles.  The  alalM  are 
■ore  Mild  hero  because  the  large  fluid  path  drives  the  optlaaa  photoa 
easrgy  to  high  values  uhere  the  aass  attenuatlea  coefficleats  are  both 
saall  aad  choage  little  with  eacrgy. 


l.S-104 


s.  ^  mL 


S««Mr1i1ii9  a  ft«  of  tht  conclusions  to  data,  tho  datcctcr  shield 
.needed  to  nake  the  neasurenent  feasible  will  be  at  least  8  In.  thick. 
k  hlQh  energy  gaoM  source  (1.5  -  2.0  NeV)  h111  give  a  louer  error 
and  a  better  signal/nolse  ratio  (snaller  background)  than  the  existing 
Cs«137  source.  The  present  design  configuration  being  oursued  Includes 
use  of  a  Hal  type  detector  In  the  pulse  counting  node,  a  single  source 
and  off-line  digital  background  subtraction.  Table  II  details  the 
developnent  work  areas  and  the  near  and  long  tens  directions. 


TABLE  II 
FUTURE  WORK 


Near  Tens 


Long  Tena 


Detector  Development  Improve  Nal 


Investigate 
Germanium  Application 


Signal  Conditioning 


Develop 

Transient  Digital 
Processor 


Investigate 
Online  Analog 
Background  Subtraction 


Densitometer 

Source 


Develop  High  Energy  Consider  Embedded 

Source(s)  Supply  Source 


Software  Development 


Calculate  Flow  Cross  Develop  Local  Denslvy 
Sectional  Average  Density  Calculation 


Background  Radiation  Lawrence  Livermore 
Determination  Laboratory  Study 

t  Evaluation 


Loft  Shielding 
Study 
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ExptrlMnts  and  prellMlnary  theoretical  results  Indicate  that 
transit  tiae  f1o««eters  using  theraocoupje^  as  sensors  are  not  useful  In 
Loss-of-Fluld  Test  (LOFT)  bloudown  experlaents.  ktork  Is  In  progress  to 
detemlne  whether  transit  time  flowmeters  using  other  types  of  sensors 
might  be  useful  In  LOFT. 
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1.  mrnooucTiOi 


Transit  tfm  flonanttrs  am  baing  contidtrad  as  an  altamative  to 
turblno  flOMWtars  for  valocity  ■aasurisnt  In  LOFT  blowdowi  experl- 
■ants. 


A  transit  tlaa  floMoeter  (TTF).  In  this  discussion,  consists  of  tao 
basic  parts:  a  pair  of  sansors  and  a  data  analysis  systea.  The  sensors 
are  located  one  downstreM  froa  the  other.  They  sense  naturally 
occurring  randoa  fluctuations  In  some  fluid  property,  such  as  teapera- 
tura  or  density*  for  uhlch  the  fluctuations  aove  with  the  fluid  rather 
than  propagating  through  the  fluid.  The  data  analysis  systaa  Is  soae 
device,  traditionally  a  correlator,  that  deteralnes  the  relative  tiae 
delay  between  the  signals  froa  the  two  sensors.  A  correlator,  for 
exanple.  geiterates  a  curve  of  correlation  aagnltude  versus  delay  tiae. 
This  curve  hopefully  has  one  doalnant  peak  which  Indicates  soae  sort  of 
average  or  doalnant  delay  tlaa  between  the  two  s1<:^ls.  This  delay  tiae 
Is  Interpreted  as  the  (average  or  doalnant)  transit  tiae  of  the  fluid 
between}  'he  two  sensors.  The  velocity  Is  trivially  obtained  by  dividing 
the  spacing  between  the  two  sensors  by  the  transit  tiae.  it  Is  conve¬ 
nient  to  u:;e  a  device  that  generates  a  curve  of  transit  tiae  versus  real 
tiae.  or  perhaps  velocity  versus  real  tiae.  Such  a  curve  will  always 
have  at  least  a  little  noise,  because  the  peak  In  the  correlation  ct.'ves 
aoves  a  little  In  real  tiae  even  for  perfectly  steady  fluid  flow.  This 
point  will  be  further  discussed  later. 

There  are  several  potential  advantages  of  transit  tiae  flowaeters 
over  turbine  flowaeters: 

1.  TTFs  aay  disturb  the  fluid  flow  nuch  less  than  turbine  flow¬ 
aeters.  depending  on  the  type  of  sensor  used. 

2.  The  hardware  Inside  the  pipe  Is  usually  auch  slapler  and  acre 
reliable  than  a  turbine,  depending  on  the  type  of  sensor  used. 

PtwWkipailbd 
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3.  In  principle,  no  calibration  Is  required  with  a  TTF.  Even  If 
calltratlon  Is  necessary,  the  calibration  prob'iens  are  Mich 
less  severe  than  with  turbines. 

4.  In  principle,  a  TTF  can  Indicate  the  distribution  of  flow 
velocities  present  In  the  measurement  region.  In  two-phase 
flow  with  slip  for  example,  a  good  TTF  could  give  separate 
Indications  of  the  gas  velocity,  the  liquid  velocity,  ar.d  the 
Interface  velocity.  (In  fact,  this  Ideal  Is  far  removed  from 
the  present  state  of  the  art.) 

Transit  time  flowmeters  offer  a  potential  for  either  very  localized  or 
global  velocity  measurements.  Including  the  potential  for  detei^lnlng 
both  the  average  velocity  and  velocity  probability  density  function  for 
an  entire  pipe  with  a  single  TTF. 


2.  PAST  WORK 

We  have  tested  an  experimental  TTF  using  thermocouples  as  sensors^^^. 
The  pipe  section  with  the  TTF  sensors  Is  sketched  In  Figure  1,  and  the 
details  of  the  thermocouples  mounting  are  shown  In  Figure  2.  These 
thermocouples  had  response  times  (0-63S)  of  about  16  msec,  for  a  roll-off 
frequency  of  10  Hz.  It  is  not  known  how  to  make  significantly  faster 
thermocouples  that  will  survive  the  required  3000  hours  In  350”C  water. 

In  an  effort  to  learn  about  the  effects  of  sensor  spacing  In  a  TTF, 
five  thermocouples  were  rsunted  In  the  pipe  with  spacings  between  pairs 
of  adjacent  thermocouples  ranging  from  0.4  to  29.7  cm.  Data  from  all 
five  thermocouples  were  recorded  simultaneously  on  analog  magnetic  tape, 
and  pairs  of  signals  were  analyzed  later. 

The  most  Important  conclusion  of  this  experimental  work  comes  from 
steady  flow,  all-liquid  data.  Figure  3  shows  four  graphs  of  Inverse 
velocity  (transit  time  *  sensor  spacing)  versus  real  time,  for  four 


Fig.  3  Transit  tine  curves  obtained  with  20«sec  averaging  tint. 


ttwr  «pKl«|s.  1l»  wiM  In  Vtm  first  fmt  iicmii  sf  son  sf 
tht  frtpks  It  M  srtlfsct  of  tl»  tystai*  Mtf  timili  ht  Ifsorsi.  TImm 
foyr  grtphs  slioiiM  all  te  tlit  sms.  If  tte  TTF  wrM  Idatlly.  If  • 
fist  vslocitr  proflls  It  sstuMd,  tlit  four  graplit  should  oil  lit  on  tht 
straight  linos  that  art  draun  In  at  2S  muc/m,  to  agrot  with  turhino 
flouMttr  MOsurtBonts.  Thoso  four  graphs  Indicate  a  good  potential  for 
Making  steady  state  velocity  MoasureMcnts  ulth  this  TTF.  Houever,  uhen 
MO  try  to  use  this  systMS  In  transient  flow,  we  cannot  use  a  20*sec 
averaging  time  In  the  data  analysis,  as  we  did  for  the  curves  of 
Figure  3.  The  averaging  tine  for  blowdown  data  nust  be  only  a  fraction 
of  a  second,  because  the  velocity  Is  not  even  approxlnately  constant  for 
tine  Intervals  greater  than  1  sec.  Figure  4  shows  the  results  of  the 
sane  data  and  analysis  procedures  represented  In  Figure  3.  except  that 
the  averaging  tine  was  reduced  fron  20  to  only  1  sec  for  the  Figure  4 
curves.  Obviously,  the  curves  of  Figure  4  are  nuch  wore  noisy,  with 
sone  noise  peaks  approaching  the  full  range  of  the  analysis  systen 
(indicated  by  the  “FULLSCALE"  notation  In  the  figure).  For  steady  flow 
data,  these  curves  can  be  snoothed  to  get  a  reasonable  Indication  of  the 
fluid  velocity.  But  for  transient  flow,  the  curves  cannot  be  snoothed, 
and  It  Is  unknown  whether  the  peaks  In  the  curve  represent  randon 
statistical  noise  or  real  velocity  fluctuations. 

An  approximate  theoretical  treatnent  predicts  that  the  noise 
(standard  deviation  of  the  curve  about  the  mean)  In  a  transit  tine  curve 
obtained  from  a  cross->correlat1on  analysis  Is  about  BOX  of  reading  for 
the  best  data  obtained  In  these  tests.  This  Is  In  reasonable  agreement 
with  the  observed  noise  level. 

The  noise  In  the  curves  of  Figures  3  and  4  Is  Inherent  in  any 
transit  time  flowmeter  that  uses  random  fluctuations  of  the  fluid 
properties.  The  effwt  of  the  noise  may  be  less  severe  In  some  types  of 
data  analysis  systems,  but  It  will  always  be  present  In  at  least  a 
qualitatively  similar  form.  The  noise  magnitude  decreases  as  the 
averaging  time  Increases  and  as  the  bandwidth  of  the  correlated  signals 
Increases,  but  for  a  given  ITT  system  In  a  transient  flow  application, 
the  parameters  cannot  be  Increased  at  will. 
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Fig.  4  Transit  tiM  curvts  obUlnad  with  1-tac  avtraging  tiM. 


It  CM  kt  concImM  tint  tkU  TTf  tyttai  It  mu  tsaftl  It  LOFT 
klwkOMt  mptrlaintt  wtltts  trrort  of  at  loatt  201  of  rootfltg  art 
accoptabla.  Tho  orosant  stato  of  kMMlad9a  dictalot  a  pottfaUtfc 
attitudo  about  possfolo  Itprou— wtt  rotultfnp  fraa  difforoot  data 
analysis  systMS.  Howavtr,  TTFs  usino  tbaraocouplas  as  sansors  should 
bo  vary  useful  In  applications  In  which  longer  averaging  tines  are 
acceptable,  and  In  less  severe  envlrewaents  In  which  faster  response 
themocoui^les  can  be  used. 


3.  COWTItlUIWS  MOWK 


The  negative  results  obtained  with  the  themocouple  type  TTF  do  not 
Inply  that  no  transit  tine  flowneter  systna  can  work  In  LOFT  blowdowns. 
At  this  time  optimism  Is  high  about  the  potential  of  any  TTF  system  In 
which  the  correlated  power  bandwidth  Is  100  Hz  or  greater.  Of  course, 
such  a  TTF  needs  sensors  with  a  response  bandwidth  of  at  least  100  Hz. 

He  are  therefore  preparing  experiments  In  which  four  alternative  sensors 
»  drag  screens.  gaoM  densitometers,  ultrasonic  void  sensors,  and 
electrical  conductivity  void  sensors  —  will  be  evaluated  In  steam/water 
blowdown  tests  In  the  Bettis  Flask  facility  at  INEL.  These  sensors 
will  all  be  designed  with  bandwidths  In  excess  of  100  Hz.  The  most 
Important  Information  to  b'i  gained  from  these  tests  will  be  the  band* 
widths  of  the  signals  obtained  from  the  fluid  under  blowdown  conditions. 
He  also  hope  to  gain  some  Information  about  the  coherence  time  or  the 
coherence  length  of  the  signal  In  the  fluid.  That  Is,  we  hope  to  learn 
how  long  the  signals  last  In  the  fluid,  or  how  far  the  fluid  travels 
before  the  signals  become  uncorrelated  or  Incoherent. 

This  experimental  Information  Is  necessary  but  not  sufficient  to 
predict  the  performance  of  a  TTF  In  transient  flow.  We  also  must 
understand  quantitatively  the  errors  In  TTF  readings  and  how  they  depend 
on  such  parameters  as  the  power  spectral  density  (cross>  and  auto-}  and 
the  record  Ic.ngth  or  averaging  time.  We  are.  therefore,  doing  a 
theoretical  study  concurrently  with  the  experimental  work.  This  study 
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iHII  cMsMtr  «M«1  posclklt  tm»  4tte  Miljrsff  iiicliitflaf 

tkt  trtdUltMl  cms-cprrtUkIw.  tiit  mt*  infiofmtlvt  (potMtIally) 
lipilM  rtipwin  function,  oni  on  o^taolly  uoltiitni  avtmgo  of  tiio 
contributions  of  soporott  fropuoncjr  conpononts. 

Fifuro  S  show  soao  oorly  moults  frm  this  thoemtlcal  stu^. 

This  Is  a  graph  of  vorsus  t^T  for  savomi  values  of  p,  uham 

«2  Is  tho  angular  froguancy  cutoff  of  tho  cormlatad  pouar 

Is'tha  standard  deviation  of  the  position  of  the  peak  In  the 
cross*conre1at1on  curve,  or  ttie  error  In  the  transit  tiM 
reading 

T  Is  the  record  length  or  averaging  tine 

p  Is  the  conrelatlon  coefficient  or  the  nomallzed  cormlatlon 
nagnitude  at  the  peak  of  the  correlation  curve. 


•»T  (angular  fraquancysMaragingiMw) 

Fig.  S  Graph  for  predicting  noise  nagnitude  In  transit  tine  curves. 
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Ill  tlilt  calcuUtloii*  tiM  cross  powr  spoctrsi  iMsI^  ms  ossMti  to  feo 
fist  for  and  zoro  for  Tiw  onto  poMor  spoetrol  doosUlos 

MTO  sssMsd  fist  for  0<M<»|,  and  xtro  far  M|<«,  wItP  m2-**1* 

In  tho  bast  cast  of  tlio  .n'cvlously  discussad  SmIscsIo  data,  m  bad 
p  ■  0.5  aod  m2  •  2«  10  Hz  ■  62.8  sac*^.  If  m  cboosa  T  ■  1  sac.  m  flod 
that  Is  batwaan  11  and  17  as,  tdileh  Is  graatar  than  201  of  roading. 

For  T  ■  20  sac,  Is  batnaan  2.5  and  3.8  as.  Thasa  valuas  are  In 

raasonabla  agreaaant  with  axparlaantally  obsarvad  valuas.  Tha  aMct 
valua  of  In  tha  speciflad  rangas  (l.a.,  tna  amt  position  In  tha 
shadad  band  on  tha  graph)  dapands  on  tha  valua  of  Saallar  values 
give  points  along  tha  lOMr  edge  of  tha  shaoad  band. 


4.  COWCLUSIONS 


It  Is  fait  that  transit  time  floMKtars  using  thenncouples  as 
sensors  cannot  be  used  In  LOFT  bloudown  experiments.  However,  TTFs 
using  other  types  of  sensors  may  be  useful.  Experimental  and 
theoretical  work  to  evaluate  the  potentials  of  TTFs  using  other  sensors 
Is  presently  being  done. 
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Out  of  tho  btsic  qumtltlM  of  lotonst  during  •  1oss*of-coolMt 
ixporloMit  (LOCE)  Is  tho  prlairy  systot  atts  flou  ratt.  Presantly. 
ttoro  art  no  transductrs  aanarclally  avallablo  ulilch  continuously 
■tasurt  this  paraaotar.  Thtrtfort,  c  ^ansductr  uas  dosigncd  at  E6I6 
Idaho*  Inc.  which  coahlnts  a  drag-disc  and  turhlno  Into  a  single  unit. 
The  basis  fbr  the  design  tas  that  the  drag-disc  would  aeasure  acatntua 
flux  (pV^)*  the  turbine  would  aeasure  velocity  and  the  aass  flow  rate 
could  then  be  calculated  froa  the  two  quantities  by  assualng  a  flow 
profile. 

For  two-phase  flow,  the  outputs  are  approxlaately  proportional  to 
the  desired  paraaeter*  but  rather  large  errors  can  be  expected  under 
those  assMptlons.  Prel  Ininary  evaluation  of  the  experlaental  two-  a^ 
single-phase  calibration  data  has  resulted  In  uncertainty  estlaates  of 
jjBt  of  range  for  the  turbine  and  ♦20X  of  ra..je  for  the  drag-disc.  In  an 
effort  to  reduce  the  errors*  further  Investigation*  were  wade  tc  deter- 
alne  what  the  drag-disc  and  turbine  really  aeasure.  In  the  present 
paper*  three  turbine  nodels  for  vertical*  two-phase*  steaa^water  flow 
are  Investloated;  the  ^a  Hodel*  the  Rouhani  Model*  and  a  volunetric 
flow  aodel^^*^*^^*  Theoretical  predictions  are  coapared  with 
experlaental  data  for  vertical*  two-phase  steaa/water  flow.  For  the 
purposes  of  the  aass  flow  calculation,  velocity  profiles  were  assuned  to 
be  flat  for  the  free-fleld  condition.  It  Is  appreciated  that  this  aay 
not  be  true  for  all  cases  Investigated,  but  for  an  Initial  Inspection, 
flat  profiles  were  assuaed. 

It  should  be  noted  that  because  of  the  physical  arrangeaent  of  the 
drag-disc  and  turbine  and  because  of  unique  blade  geometry,  the  results 
should  be  applied  with  caution  to  other  turbines.  See  Figure  la  for 
transducer  geometry. 


[a]  This  paper  Is  condensed  from  Reference  1. 
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Fig.  la  LOFT  drag-disc  turbine  transducer  geometry  (shroud  reiaoved). 


The  experimental  data  employetJ  In  the  Investigation  itere  obtained 
at  Uestinghouse  Canada  Limited  (UCL)  In  two  series  of  tests^^^.  In  one 
series,  all  flow  went  through  the  turbine;  this  Is  denoted  as  Wie  "full 
flow"  series.  In  the  other  series*  a  4. 05- In.  ID  pipe  was  employed; 
this  Is  denoted  as  the  "free-fleld"  series  since  flow  could  bypass  as 
well  as  pass  through  the  1.18-1n.  diameter  turbine.  The  ranges  of  the 
test  variables  are  suauaarized  in  Table  I. 
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2.  DESCRIPTION  OF  AYA‘S  MODEL 


A  model  of  the  behavior  of  a  turbine  meter  was  proposed  by  Aya^^^. 

It  Is  assumed  the  flow  Is  dispersed  with  the  steam  moving  faster  than 
the  water.  The  forces  acting  on  the  blade  are  visualized  as  a  steam 
force  acting  on  the  upstream  portion  of  the  blade  and  a  resisting  water 
force  acting  on  the  downstream  portion.  The  turbine  is  assumed  to  be 
turning  at  a  speed  that  corresponds  to  a  fluid  velocity  which  lies 
between  the  steam  and  water  velocities.  It  is  assumed  that  the  steam  is 
trying  to  drive  the  turbine  Into  the  water  and  the  water  on  the  other 
side  of  the  blade  resists  this  driving  force.  The  derivation  presented 
In  Reference  2  is  followed  below. 

The  blade  is  shown  in  Figure  lb.  The  tangential  velocity  at  a 
point  is  ru.  where  r  is  the  radial  distance  from  the  ax',  to  be  the 
point  and  y  is  the  angular  freguency  of  the  blade. 

T;.e  gas  velocity  is  usually  greater  than  o**  equal  to  che  l'*quid 
velocity  in  two-phase  flow.  For  steady  state  flow,  the  forces  acting 
across  the  blade  and  normal  to  the  turbine  axis  of  rotation  must  be  in 
equilibrium.  Assuming  that  the  forces  acting  normal  to  the  b'ade  are 
drag  forces  and  that  viscous  forces  can  be  neglected,  the  equilibrium 
equation  is  as  follows  (see  Figure  2). 

0/2)  Cjj  Ag  (Vj- (1) 

•  (1/2)  »f  Af  (V,  - 

where 

can  be  obtained  either  from  blade  geometry  or  an  all  water 
calibration. 
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Fig.  2  Velocity  vectors  for  the  Ay*  nodel. 

Equation  (1)  slnpllfles  to 

«,  (»,-  *»)*  c,  ■  (1  - (Vj-V,)* 

For  the  foregoing  equations  the  following  are  defines: 

A^  >  flow  area  occupied  by  the  stean 

A  •  total  flow  area 

A„ 

a  a  -21  (void  fraction) 

A 

Pg  ■  MSS  density  of  steaa  at  test  tcwperature 
•  velocity  of  steaw 

Cg  "  drag  coefficient  for  turbine  blade  In  steam  flow 


■  MM  iMilly  vf  IlfiM  %t  tMt  tMurtfri 


if  *  vtloelty  «f  IImM 

Cf  •  4rt9  cecffIciMt  for  turtiM  blodo  In 

S  •  (slip). 

V, 


Dividing  both  sIHts  of  tho  cquatlen  by  if  and  rearranging  yields 


This  Is  one  equation  with  six  paraMters  so  that  further  relations  are 
required.  Once  five  of  these  paraMters  are  deteralned.  Equation  (1) 
can  be  solved  for  V^.  the  fluid  velocity. 

The  drag  coefficients  depend  on  Reynolds  nuober  and  the  cavitation 
nMber  and  not  the  type  of  fluid.  Assiaalng  the  ratio  of  drag  coefficients 
Is  one.  Introducing  the  slip  ratio  S  ■  i^if»  and  solving  the  quadratic 
equation  yields 


Za _ 

« 

)  -C. 

^l-o)  -  oS 

♦<  |»S  -  (1-a)  — 

\l  "gJ 

(l«o)  -  aS^ 

L  ‘*9  J 

|l/2 

Knowing  the  fluid  velocity*  the  gas  velocity  can  be  calculated  from 
the  slip  ratio.  The  mss  flow  rate  can  then  oe  computed  as 


(3) 


\  •  ♦  (l-)p^Vj  A  (4) 

«litr«  A  Is  tht  total  flow  oroa.  It  list  bato  astuatd  that  tha  valocity 
prof 11a  Is  flat  ovar  tha  araa. 


3.  WOUHANPS  won. 


Rouhanl^^^  davalopad  a  aodal  basad  upon  physical  masoning  that  tha 
highar  valocity  staaa  Is  driving  tha  turblna  Into  tha  slowar  valocity 
liquid.  Qy  application  of  tha  llnaar  ■oosntuw  aquation  with  no  body 
forces*  tha  following  governing  aquation  was  obtalnad^*^. 


».  •  *t » ‘rrWi 

wham  X  Is  defined  as  tha  flow  quality 


(5) 


(6) 


and 

Ng  ■  MSS  flow  rata  of  gas 
N^«  total  MSS  flow  rata 
By  awploylng  the  ralatlon^^^ 


^  (1-x)  «  Pf 

[a]  The  aquation  pmsented  In  Rafamnca  3  can  be  recast  Into  tha 
a<piat1on  pmsantad. 


(7) 
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Hu  fwtniltf  •qpMtton  cm  ta  rcurl^M  m 


V 


f 


•$  ♦  (1-a)  ^ 

_ !i 

.S*  ♦  (l.a)  ^ 


rgi 


(8) 


4.  VOLUMETOIC  FLOW  MODEL 


This  Model  assunes  that  the  turbine  measures  volumetric  flow  rate 
regardless  of  whether  the  flow  Is  single*  or  two-phase.  Thus,  the 
volumetric  flow  rate  Is  given  as 

Q  •  AgVj  ♦  •  <»AVg  ♦  (l-a)  (9) 

where  IR  Is  the  flow  area  through  the  turbine.  Introducing  slip  yields 

Q  •  *  [1  ♦  a  {S-l)l  . 

and  solving  for 

V  - - 9 -  . 

5.  ALL  WATER  CALIBRATION  DATA 


(10) 


(11) 


Prior  to  the  two-phase  tests,  all  water  calibration  tests  were 
conducted.  The  results  of  the  tests  are 


.-0.116+  34.47 
'0.344  ♦  39.19  Iq 


Full  Flow  (12) 

Free  Field 
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Winn  It  Urn  tranit  fiM  vtlocl^y  1w  tlwi  pipt  tni  It  tin  tiirblM 
owtpwt  mtt. 

For  til  witor  flow*  tin  tliroo  attelt  tiaplify  to  (•■0*  S-l)'. 


»t 

•*« 

Rya*t  and  Rouhanl's  aodel 

(13) 

ft 

Voliantrlc  flow  aodel 

JT 

av 

6.  MOB.  CCMPARISOII  MTM  EXPERlHBiT 

Ctleuittlont  for  tin  tliroo  ndols  won  Mdo  ot  follows: 

pfi  Pg  ■  obtalnod  frtn  ASNE  stoaa  tables^^^  for  tbo  tost 
towpontun 

S  ■  obtoliwd  froa  Roforonco  8*  Flguro  3.  os  a  function 
of  p^Pg 

X  ^  flow  quality  known  froa  tbo  tosts 

« 

«  •  obtalnod  froa  Equation  (7)  onco  S  and  x  aro  known 
«  froa  all  wator  calibration*  Equation  (13). 

Tbo  rangos  of  tbo  tost  varlablos  asy  bo  scon  In  Table  l. 

7.  RESULTS 


If  no  aodcling  wore  conducted  tbon  tbo  basic  calibration  data  for 
tbo  turbine  would  bo  as  shown  In  Flguro  3.  The  UCL  velocity  was  ca1> 
culated  froa  the  relationship 
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Fig.  3  Turbine  calibration  with  no  wodellng. 


where  Is  the  average  fluid  density  and  was  obtained  fron  measured 
test  parameters  and  A  Is  the  net  flow  area  at  the  Instrunent. 

The  effect  of  modeling  on  the  calibration  may  be  seen  In  Figures  4 
and  5.  In  Figure  4  Is  shown  the  liquid  velocity  versus  turbine  volts 
for  Rouhani's  model.  Steam  velocity  versus  turbine  volts  Is  shown  for 
slips  of  1.6  In  Figure  5a  and  Z.3  In  Figure  5b.  Similar  results  are 
obtained  for  the  other  two  models. 

The  best  test  for  the  model  Is  a  comparison  of  predicted  versus 
measured  mass  flow  rate.  This  Is  done  for  the  free  field  tests  In 
Figures  6.  7,  and  8  for  the  Rouhani,  volumetric  flow  and  Aya  models. 


I3M/HI  UIWIS 
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Fig.  4  Liquid  velocity  calibration  froa  the  Rouhanl  aodel. 
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Fig.  6  Predicted  versus  experlwntel  ntss  flow  rate  frow  the  Rouhenl  aodel 


rttptctivaly.  lh»  4tU  frta  ttmam  flfiirM  aai  froa  tlM  full  flow  tMts 
V  vmmriiM  In  Takit  11. 


TWLE  n 

RESULTS  FROM  THE  THREE  TURBINE  NOOELS 
FORIMRO  FUM  KITH  0RA&-DISC  UPSTREAM 


Model 

c  M 
‘•o 

c 

*•1 

Test 

Configuration 

0  [‘>3 
yx 

2o_  X  100 

e.  rci 

30. 

Rouhani 

-  0.186 

1.215 

Free  Field 

3.02 

9.9 

Volumetric  Flow 

-  0.648 

1.033 

Free  Field 

1.93 

6.3 

Aya 

-  0.369 

1.169 

Free  Field 

2.62 

8.6 

Rouhani 

0.033 

0.913 

Full  Flow 

0.216 

— 

Volumetric  Flow 

-  0.028 

0.761 

Full  flow 

0.311 

— 

Aya 

0.036 

0.870 

Full  Flow 

0.198 

•  •• 

w  Vld.1* ‘o  ♦ ‘l**!' 

[b]  The  estlMted  standard  deviation  of  y  on  x  Is  c  . 

y* 

[c]  The  quantity  30.6  Is  the  naximan  flow  rate  for  the  free  field  tests. 


For  the  free  field  tests t  the  agreeaent  between  the  predicted  and 
Measured  aass  flow  rates  Is  surprisingly  good  for  all  three  models. 
However,  the  mass  flow  rate  predicted  from  the  volumetric  flow  model 
agrees  best  with  the  experimentally  measured  mass  flow  rate.  This  Is 
evidenced  by  the  closeness  of  the  coefficient  to  1.0  and  the  small > 
ness  of  a  . 

For  the  full  flow  tests.  Aya's  and  Rouhanl's  models  agree  with  each 
other  very  closely  and  both  agree  better  with  the  experimental  data  than 
does  the  volumetric  flow  model. 
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Hn  •KMlMtlOfi  of  tho  rosuUs  sIiom  that  for  the  free  field  the 
errors  In  mss  floe  rate  are  ^9.9.  ^.3.  and  ^.6X  of  test  range  for  the 
Moiihanl.  voluMtrIc  flow,  and  Aya  wdels,  respectively 


8.  CONCLOSIOWS 


The  three  turbine  models  all  predict  mss  flow  rates  which  are  In 
surprisingly  good  agreement  with  the  measured  rates.  For  the  purposes 
of  the  LOFT  experiments,  the  free  field  tests  are  more  representative. 
For  these  tests,  best  results  were  obtained  using  the  volumetric  flow 
model. 


Presently,  work  Is  continuing  with  the  models.  From  preliminary 
data  It  appears  that  the  Rouhani  and  Aya  models  correlate  the  liquid 
velocity  and  steam  velocity  better  than  the  voluMtrIc  flow  model. 

Linear  extrapolation  of  the  predicted  mss  flow  data  from  the  UCL 
to  LOFT  size  piping  Indicates  that  large  errors  My  occur. 

From  the  above  It  Is  seen  that  none  of  the  models  appear  to  be 
significantly  more  accurate  than  the  others  and  that  caution  should  be 
utilized  In  applying  the  turbine  models  to  large  pipes. 
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SURVEY  OF  FACILITIES  FOR  TWO-PHASE 
FLOW  INSTRUMENT  TESTING 


I.  C»togorie«  of  Two-Phase  Flow  Initrumont  Taat» 

A  limited  review  is  being  conducted  of  facilities,  in  the  United 
States  and  abroad,  which  have  capabilities  for  testing  the  flow  instruments 
used  in  water  reactor  safety  experiments.  Such  experiments  require  flow 
instruments  which  determine  the  mass  flows  of  each  fluid  phase  (steam  and 
water)  and  a  flow  regime  characterising  variable,  such  as  void  fraction. 

Because  mass  storage  effects  are  significant,  flow  instruments  are  required  at 
several  locations  in  a  typical  safety  experiment.  At  each  location,  a  num'ier  of 
measurements  of  differing  principle  must  be  made  to  define  the  mass  flows  and 
void  fraction:  the  number  corresponds  to  the  number  of  unknowns.  Thus,  a 
flow  "instrument"  which  is  to  be  tested,  may  actually  consist  of  several  measur¬ 
ing  devices,  from  which  individual  phase  mass  flows  and  regi**'e  can  be  inferred. 

Testing  such  instruments  requires  facilities  have  a  vairiety  of 
characteristics.  The  variety  derives  from  the  scope  of  the  tests  needed  to  de¬ 
velop  and  calibrate  a  flow  instrument.  The  tests  are  as  follows: 

a.  Conventional  Calibration  Tests 


All  of  the  individual  measurement  devices  making  up  a  flow  instru¬ 
ment,  whedxer  the  instrument  is  the  first  of  a  kind  or  a  follow-on 
unit,  must  be  subjected  to  limited  (conventional)  calibration  testing. 
These  "production  calibration"  tests  establish  the  "meter  factors" 
of  the  measurement  devices  by  observing  the  response  under  flow 
conditions  where  the  nature  of  the  response  is  well  understood 
(usually  single  phase  flow).  Separate  teats  in  both  gas  and  liquid 
phases  may  be  necessary  to  characterize  a  device  fully.  For  ex¬ 
ample,  large  liquid  phase  flows  may  be  needed  to  achieve  rated 
momentum  flux  while  small  gas  phase  flows  may  be  needed  to 
characterize  minimum  velocity  measuring  capability  (turn-down). 
"Meter  factors"  established  in  this  manner  are  used  as  parameters 
in  the  equations  which  model  the  responses  in  two-phase  flow. 


In  addition  to  measuring  meter  factors,  conventional  calibration  tests 
are  needed  to  check  the  repeatability,  linearity,  and  accuracy  of  each 
sensing  device  manufactured. 


b  order  to  combino  output*  of  mooouring  dovteoe  to  cotculoto  tko 
poromotoro  of  two«phoae  flow  in  onglnoorlng  units,  on  onolyticol  or, 
possibly,  on  sinpiricol  model  of  the  response  of  eocb  messuring  de« 
vice  in  two'phase  flow  is  essential.  Usually,  a  plausible  model  can 
be  postulated  on  theoretical  bases.  The  accuracy  of  the  model  must 
then  be  examined  by  an  extensive  series  of  tests  with  a  prototype  flew 
instrument.  In  outline,  the  test  procedure  consists  of: 

(1)  Observing  the  response  of  the  measuring  devices  of  the  prototype 
flow  instrument  over  the  full  range  of  hydraulic  conditions  ex> 
pected  in  service. 

(2)  Calculating  phase  mass  flows  and  other  flow  parameter*  o1 
interest  from  these  responses  again,  over  the  full  hydrs'^.ic 

range. 

(1)  Comparing  the  flow  parameters  calculated  from  the  v.  ]  of 
the  instrument  under  test  with  the  same  flow  paran  ie- 

termined  by  independent  means.  The  independent  tm  aS^  are 
flow  measuring  standards  of  known  precision  for  the  conditions 
of  measurement.  Typically  the  standards  will  include  meas* 
urements  of  the  individual  phase  flows  upstream  -  (  a  mixing 
device,  since  no  two-’phase  flow  measurement  standard  exists. 

Correlation  tests  are  of  the  nature  of  a  calibration  test  !  r  a  proto* 
type  instrument  under  two-phase  flow  conditions.  Hou  r-v  er,  not 
all  of  the  conditions  of  the  ultimate  instrument  inatallai.on  may  be 
duplicated  simultaneously  (i.e. .  in  a  single  test).  On  the  contrary, 
the  calibration  testing  of  the  prototype  instrument  will  in  ).eneral 
require  the  separation  of  effects  which  may  alter  instruiiient  re¬ 
sponse.  e.g. ,  it  is  desirable  to  examine  the  effect  of  temperature 
and  hydraulic  conditions  separately. 

Correlation  testing  may  also  be  used  as  an  adjunct  to  the  design  of 
the  instrument  itself.  For  example,  .\n  analytical  model  of  an  in¬ 
strument  may  assume  a  spatially  uniform  distribution  of  phases  over 
the  flow  cross  section  (though  net  necessarily  equal  velocities).  The 
correlation  testing  might  be  used  to  select  a  means  for  homogeniz¬ 
ing  the  flow  to  obtain  spatial  uniformity. 

The  products  of  correlation  testing  of  prototype  instruments  are: 

(1)  An  experimentally  confirmed  matneR.aticai  correlation  whereby 
the  flow  parameters  of  interest  can  be  calculated  from  measured 
instrument  responses,  for  all  thermal  and  hydraulic  conditions 
of  the  experiment  in  which  the  instrument  is  ulti.nately  to  be  used 

(2)  Quantitative  estimates  of  the  uncertainty  in  the  flow  parameters, 
as  established  via  the  correlation,  for  the  range  of  thermal  and 
hydraulic  conditions  of  the  experiment. 
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bratioM  (or  "motor  loctoro")  lor  oao  to  towphooo  oiqtortoMoto. 
TImoo  mo4iAbro  altar  callbratioaa  team  llioao  maaourod  tta4or 
toa  '^rodaetiaa  calibratioa"  caoditiona  of  (a)  abeva  and  ara 
datormtoad  by  too  carralatiea  toatfag  to  ba  ^praprtoto  far  tba 
taro'phaao  flow  ceoditioM* 


€•  blachaatcal  Daaaleemaat  aad  Proof  Taatiaa 

la  additioa  to  pravidtog  totarpratabla  outputo  oadar  tho  toarmal  aad 
bydraolic  coaditioiia  of  tha  asparimaato  to  which  thay  ara  uaad*  flow 
iMtruinaato  must  withatana  tha  machaaical  loadinga  impoaad  by  tha 
ai^arimaot  and  moat  andura*  without  malfunction,  for  a  tima  period 
datormtoad  by  tha  langth  of  tha  axparimant  and  tha  aceaaaibillty  of  tba 
toatrumanta  for  matoiananca  following  it.  Tha  capability  of  tha  to* 
atrumaota  to  withatand  tho  aiqparimantal  anvironmant  muat  ba  dam* 
onatratod  fay  taating  to  facilitiaa  which  duplicate  that  anvironmant. 


d.  Scalina  Taata 

Calibratioa  and  corratation  taata  normally  raquira  that  all  of  tha 
flow  ba  dlracted  through  tha  maaauring  iaatrumant.  Thia  proca* 
dura  mtoinusaa  uncartaintiea  whan  tha  flow  paramatara  computed 
by  tha  inatrumant  ara  compared  agatoat  thoaa  detarmi^^  1  via  tha 
atandarda.  Tha  intandad  uaa  of  tha  flow  inatrumant.  h.‘'>avar.  may 
tovolva  only  the  aampling  of  a  large  flow  field  with  a  amall  inatru* 
mant  (aa  with  tha  LOFT  drag  diac  turbine  tranaducar).  Taating 
may  be  required  to  eatabliah  proceduraa  whereby  tha  phaaa  (Iowa 
aamplcd  by  the  matrument  can  ba  extrapolated  to  infer  phaaa  flows 
for  the  pipe  aa  a  whole.  Tha  product  of  auch  taata  ia  a  aat  of  acal* 
tog  factors  for  the  instrument*computed  local  flows,  and/or  a 
quantiuttve  estimate  of  tha  uncartatotiaa  to  tha  total  flows  infer* 
rad  from  amall  samples,  ^additionally,  tha  physical  location  of  a 
flow-sainpllug  instrument  C'ntld  ba  aalected  in  a  scaltog  teat. 


Another  form  of  scaling  test  could  arise,  if  a  amall  model  flowmeter 
were  built  (or  correlation  testing,  but  a  larger  version  of  the  same 
mater  ware  intended  for  the  safety  axparimant.  Xn  this  situation, 
full  scale  single  *phasa  (production)  calibration  taating  would  un* 
doubtedly  ba  necasaary. 

a.  Transient  Taata 


Many  of  the  water  reactor  safety  experiments  which  require  tero* 
phase  flow  instruments  ara  of  a  transient  nature  **  flow  changes, 
sometimes  rapid,  ara  inherent  in  tha  c?q>eriment.  For  this  reason. 


1.6-3 


•m  vndcrtlMiilaf  is  Msdsd  of  tko  foctoro  offoctlag  tfio  troMlool 
rsspeosoo  of  tho  dovleos  moklog  tho  flow  taoUnmoot.  This 
koowlodfo  COM*  in  cortoiB  insUneos*  bo  osod  to  iafor  octuol  flows 
■soro  sceurstsly  from  instrumoat  iadlesUons;  ia  otbors*  tho  uo« 
eortolatios  sssociotod  with  tho  trsaslsot  mossuroaMat  esa  bo 
booadsd. 

Aa  uadsrstsndinf  of  tastrumaat  trsMior>t  rospooso  roquiroo  s  dy- 
nsmic  model  for  each  elomsat.  Tho  objectives  of  treasient  tests 

ores 

(1)  To  coafirm  the  form  of  the  model. 

(2)  To  determine  the  numerical  value  of  certain  conetaats  of  the 
model  which  are  operative  only  ia  tranaient  (e.  g. ,  effective 
rotary  inertia  of  a  flow  turbine). 

2.  Facility  Reouiremente 

Facility  requirements  depend  on  the  specific  category  of  a  test 
to  be  performed:  no  single  facility  ia  likely  to  satisfy  the  requirements  of  all 
categories.  Also  it  does  not  appear  that  one*  unique  unalterable  set  of  fa> 
cilities  must  be  found;  there  will  generally  be  several  combinations  of  facil¬ 
ities  capable  of  filling  testing  requirements.  For  example,  it  appears  de¬ 
sirable  to  do  much  of  the  correlation  and  prototype  calibration  testing  in 
gas/water  loops,  but  some  of  this  tasting  might  ^  done  (albeit  with  more 
difficulty)  in  a  well  designed  steam  water  loop. 

The  types  and  characteristics  of  facilitites  which  can  satisfy 
the  requirements  of  etch  testing  category  are  outlined  below: 

a.  Conventional  Calibration  Test  Facilities 


A  water  loop  and  a  separate  gas  loop  are  needed  to  exercise  an 
instrument  over  the  full  range  of  velocities,  nnixture  densities,  and 
momentum  fluxes  it  may  see  in  service.  The  experimental  set  up 
should  be  full  flow;  that  is,  the  Quid  flowing  through  the  instrurt-nt 
should  have  passed  entirely  through  the  standards  against  which  the 
instrument  is  compared.  The  capability  to  determine  quantitatively 
the  effects  of  temperature  on  instrument  calibration  is  also  needed; 
this  can  be  obtained  if  the  water  loop  can  be  operated  at  tempera¬ 
tures  up  to  the  maximum  expected  in  service.  Standards  against 
which  the  measuring  devices  are  compared  should  be  of  high  qual¬ 
ity  (absolute  accuracy  better  than  1%)  and  traceable  to  NB5,  to 
reduce  insofar  as  possible  the  uncertainties  due  to  this  source. 
Some  or  all  of  the  production  test  facilities  might  be  located  at  the 
instrument  manufacturer. 


k.  Cfyl»tioii  {PT»tot»D»  C>Ubr%tioal  T»«t  F*ctUtl»« 

As  liss  bsss  iodlestsd*  ths  tspsrsrion  of  tsmpsrstairo  from  hjrdrsuHe 
offsets  mokos  dssirsbio  ths  oso  of  tss/wotsr  loops  ia  dovoloping  colt- 
tootions  of  o  sot  of  moosurlag  dsvlcso  la  low-phoso  flow.  To  dupllcoto 
kpdroallc  phoaemoso  fttlly*  ths  leUoorlag  choroctoristico  oro  aoodod 
la  fos/ofotor  eoUbrotioB  loops: 

(1)  Ths  ability  to  gsaorats.  •imuttaaaoualy,  phass  vstocitioo  and 
qualltloa  which  dupUcato  those  oxpsctsd  ia  sorvica.  Ths  locus 
of  velocltias  and  qualities  will  dspsad.  of  course,  oa  the  expected 
lastrumeat  locatioa.  A  limitinf  example  for  a  LOFT  broken  loop 

is  shown  ia  Figure  1.  In  this  figure,  the  ordinate  is  the  approximate 
velocity  of  each  phase.  (In  the  broken  loop,  phase  velocities  are 
ospectad  to  remain  nearly  equal  throughout  the  blowdown,  because 
of  tbs  high  velocity  of  the  dominant  phase.) 

(2)  The  ability  to  generate,  simultaoeoualy,  momentum  fluxes  and 
qualities  simitar  to  those  expected  in  service.  As  with  velocity, 
ttis  locus  will  depend  oa  specifics  of  instrument  location  and  ex¬ 
periment;  a  limiting  example  for  a  LOFT  broken  loop  is  shown 
la  Figure  2.  Ia  tits  figure,  one  ordinate  is  void  fraction- 
weighted  momentum  flux 

This  parameter  is  believed  to  govern  the  response  of  drag- 
measuring  devices  in  two-phase  flow.  A  second  ordinate  is  the 
void  fraction,  which  with  quality,  is  necessary  to  characterize 
the  two -phase  mixture  spatially. 

These  capabilities  infer  the  selection  of  a  gas  with  a  density  r».t;io 
ral^iHve  to  water  similar  to  that  which  will  be  encountered  in  service. 
While  useful  testing  can  be  accomplished  with  air  water  mixtures  at 
pressures  near  atmospheric,  fr'l  range  coverage  requires  pressuri- 
aation  to  about  30  atmospheres,  or,  alternatively,  use  of  a  heavier 
gae. 

The  use  of  a  gas  ’  •  ch  as  sir,  rather  than  steam  for  this  phase  of  test¬ 
ing  does  not  only  <ierive  from  the  need  to  separate  temperature -induced 
effects.  It  also  permits  a  visual  examination  of  the  flow  regime  seen 
by  the  instrument  since  pipe  walls  can  be  made  transparent.  This  ex¬ 
amination  may  be  useful  not  only  in  characterizing  instrument  response 
as  a  function  of  regime,  but  also  in  determining  the  effectiveness  of 
homogenising  devices  in  obtaining  a  spatially  uniform  phase  distribution. 
A  homogenizer  is  often  necessary  to  a  repeatable  and  predictable  in¬ 
strument  response,  and  ia  an  important  part  of  instrument  development. 
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Dvpsadiag  m  tbs  ulttmal*  appllcAtiea  of  tbo  lostJfwmoat.  copobilltf 
to  mock-up  hydraulic  coadikloac  may  bo  aocoacavy.  If  tho  iacirumoal 
aador  tost  is  iatondod  for  aa  caportmoat  ia  which  it  measures  full  flew, 
Ihoa  duplication  of  bends  and  otter  hydraulic  fcatuxos»  upstream  and 
inmiediatoly  downstream,  is  aeceasary  to  ensure  that  the  distortions 
of  velocity  and  void  fraction  profiles  which  these  fsatures  will  bring 
about  ia  service  are  duplicated  in  the  calibration.  Duplication  of 
gravitational  forces  may  also  bo  important  ia  some  two-phase  flow 
situations.  This  requires  that: 

(1)  Bends  and  other  hydraulic  features,  as  wall  as  the  flow  element 
itself,  be  oriented  relative  to  the  earth's  surface  as  they  will  be 
ia  the  experiment. 

(2)  Scaling  of  the  instruments  and  hydraulics  be  full;  extrapolation  of 
the  two  dominant  flow  field  forces,  inertial  and  gravitational,  to 
another  scale  requires  preeervation  of  the  Froude  number,  but 
maintaining  Froude  number  is  not  consistent  with  maintaining 
fluid  velocity  and  momentum  flux  up  to  the  ratings  of  the  sensing 
instruments. 

The  impact  of  temperature  on  the  calibration  of  many  flow  instruments 
is,  generally  speaking,  the  result  of  changes  in  instrument  dimensions 
owing  to  thermal  expansion  of  individual  parts  and  of  the  instrument  as 
a  whole.  Calibrations  of  readout  devices,  particularly  strain  gages.can 
also  be  affected  by  temperature.  In  these  situations,  the  evaluation  of 
temperature  effects  can  be  carried  out  in  single-phase  fluid,  operating 
over  the  temperature  range  of  interest. 

If,  however,  the  physical  phenomena  to  which  a  measuring  element 
responds  is  sensitive  to  temperature,  testing  in  two  phases  at  elevated 
temperature  may  be  necessary.  In  any  case,  elevated  temperature 
tests  in  two  phases  can  be  a  desirable  confirmatory  adjunct  to  other 
calibration  testing.  Such  tests  might  be  performed  in  conjunction  with 
mechanical  tests  discussed  in  the  next  paragraph. 

c.  Mechanical  Development  and  Proof  Testine  Facilities 

Obviously,  tests  carried  out  for  purposes  of  calibration  (categories  a 
and  b,  above)  also  serve  to  a  degree,  to  test  the  capability  of  an  in¬ 
strument  to  withstand  the  mechanical  loadings  to  which  it  will  be  sub¬ 
jected.  However,  conclusive  proof  testing  requires  the  simultaneous 
imposition  of  mechanical  loads  a  ■  i  /.ermal  effects.  This  leads  to 
tiie  requirement  that  final  proof  be  carried  out  in  a  facility 

capable  of  supplying  steam  water  -tures  at  pressures,  temperatures, 
flows  and  qualities  spanning  the  expected  service  range.  The  capability 
for  mixed  phases  is  particularly  important,  since  the  random,  rapid 
variations  in  void/Uquid  distributions  characteristic  of  steam  water 
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rnfartiiraa  caa  raault  ia  vary  aavara  vibratory  loads  oa  boarief  •  aad 
•tbar  roaebaalcal  parts.  As  adtb  ealibratloa  tssts*  mochaaical  proof 
tasts  rs<|uira  tho  facility  to  supply  phase  valocitios  aad  momaatum 
flaaas  up  ta  tha  maaimum  capacity  of  tha  iastrumaat,  but  aot  mora. 

If  aa  iastrumaat  sanplas  a  Aow  fiald  (as  doas  tha  LOFT  OTT),  this 
maaas  that  proof  tasting  could  ba  aecomptishad  ia  a  facility  with  tha 
capability  to  supply  that  fraction  of  tha  total  mass  aad  vohuna  fluxes 
sasa  by  tha  instrumant  itself.  It  is  coasidarad  that  tha  machaaical 
daslga  of  tha  supporting  strut  for  such  aa  instrumeat  can  ba  made  suf¬ 
ficiently  conservative  to  eliminate  tha  nead  for  testing  in  a  two-phase 
facility  (though  such  tasts  would  remain  desirable). 

d.  Scalins  Test  Facilitiae 

If  a  small  flow  instrument  is  used  to  measure  a  large  flow  flald  (as 
la  LOFT),  scaling  teats  may  be  dasirabla.  Tha  principle  quastione 
to  ba  answered  by  such  sealing  testa  relate  to  the  interplay  of  inertial 
aad  gravitational  forces  on  the  phase  velocity  profiles  aad  oa  tho  dis¬ 
tribution  of  voids  [for  most  water  reactor  experiments  viscous  forces 
are  believed  to  play  a  minor  role).  Many  of  these  questions  are  best 
answered  visually  -•  by  high  speed  photography  or  similar  techniques. 
The  force  field  and  tha  ability  to  measure  it  can  best  be  achieved  by 
a  large  scale  gas /water  facility  (as  opposed  to  a  steam  water  facility), 
la  addition,  detailed  void  maps  and  other  flow  field  information  may 
be  obtainable  by  laser  or  hot  wire  anemometry.  minature  conductivity 
probes  and  similar  instruments,  all  of  which  can  be  readily  applied 
la  gas/water  experiments  but  are  difficult  to  use  in  steam  and  water. 

The  sealing  test  facility  must  have  the  capability  of  operating  with 
gas/water  density  ratios  similar  to  those  expected  in  the  safety  ex¬ 
periment.  It  must  also  be  adaptable  with  respect  to  hydraulic  geometry; 
Individual  tests  duplicating  the  hydraulics  of  each  instrument  location 
may  be  needed. 

Although  some  measurements  of  individual  phase  flows  prior  to  mixing 
are  aecessary  to  ensure  the  conditions  of  the  experiment  hav«.  been 
duplicated,  great  precision  in  these  measurements  is  probably  un¬ 
necessary.  Operation  of  the  flow  instrument  itself  in  the  scaling 
experiment  while  desirable  may  not  be  necessary. 


e.  Transient  Test  Facilities 

Because  of  the  time  response  of  die  measuring  devices  is  likely  to  be 
a  function  of  void  fraction,  the  transient  test  facility  must  be  capable 
of  operating  with  mixed  (two-phase)  flows.  Again  a  gas/water  facility 
appears  better  able  to  satisfy  requirements  than  a  steam-water  facility 
llie  capacity  of  the  facility  should  be  sufficient  to  supply  maximum 
gas  and  liquid  phase  flows  to  the  flow  instrument,  proper,  (Transient 
testing  of  a  sampling  type  of  Instrument  in  a  full  size  pips  is  not  con¬ 
sidered  necessary. )  A  key  feature  of  a  transient  test  facility  is  the 


to  g«Mrato  vapii  cluuif  ••  to  Iha  vato  flow  of  aitoar  plMM#. 
Tkaaratleat  datarmtoattoa  af  toatfamaad  raapaaaas  to  atop  Horn 
alMagaa  la  mauallp  atralghtfarwanl;  It  la  aceoHlafly  daalraUa  to 
maaaara  raapeaaaa  to  aaparlanaala  arhieli  coma  cloaa  to  daplicattof 
atopa.  Favtkarmora,  acearato  datarmtoattoa  of  toatnunaat  raapeaaa 
maaaa  UmI  ana  muat  maasura  tha  flaw  dlaturbaaca  arlkh  a  ataadard 
maeli  mera  rapid  thaa  tha  iastramaat.  Bat  tha  raapaaaaa  af  tha  maaa 
artof  davlcaa  ma'dag  ap  tha  tostrumaat  may  ba  cjctranMly  rapid  —  a 
law  miUiaacaada,  aad  ftodiag  staadarda  aatfielaatly  fact  will  ba  dlffi- 
calt.  Caaaaquantly,  tha  moat  faacibla  way  ta  produce  a  kaowa,  rapid 
flaw  chaaga  ia  by  maaaa  of  a  quick  acting  «'al  /a.  Tha  fluid  ayatam  ar 
raagamaat  muat  ba  auch  that  tha  chaaga  ia  phaaa  flow  doaa  aot  aacaa- 
aarily  procaad  to  a  care  fiaal  valua:  toatnunaat  raapoaaa  charac- 
tariatica  may  ba  aanaitiva  to  initial  and  fiaal  phaaa  flow  valuaa. 

Facility  &>rvav  fPraliminary  Statual 

Tka  approach  to  tha  aurvay  of  facilitiaa  la  to: 

(a)  Identify  by  litaratura  aaarch  aad  other  maaaa,  candidate  facilitiaa  for 
yartoua  taat  eatagoriaa. 

(b)  By  talaphoaa  coayaraation,  aatabliah  a  preliminary  eat  of  facility 
charactariatica  including: 

(1)  Williagaaaa  to  do  ’Contract  calibration'*  work,  contractural  re> 
quiramanta,  rad  coating  rataa. 

(2)  Fluida  amployad  [a.  g. ,  gaa,  water,  gaa/watar,  ataam,  ataam/ 
water] 

(3)  hfaximum  operating  praaaura  and  tamparatura 

(4)  Nominal  aiaa  of  loop  and  matariala  of  conatructioa 

(5)  Maximum  capacitiea  of  i.oop  ia  tarma  of  phaaa  maaa  flow  rataa, 
▼alocitiaa,  and  apparent  momaatum  flux,  veraua  quality 

(4)  Type  aad  accuracy  of  ataudarda  including  torn>down  range 

(7)  latarfaca  roquiremanta  between  tha  device  under  taat  and  the 
loop  (a.g. ,  acceptable  methoda  of  cloaura,  aupport  facilitiaa 
available) 

(8)  Schedule  lead  time  aad  notifleation  raquiramanta 

(9)  Capability  to  change  hydraulic  geometry;  coat  aixl  achedular 
impact  of  auch  changoa 


<I0)  C«pabUlt)r  to  eeiiatnict  ad  hoc  hfdraallc  medala  («.  g. ,  reactor 
ooaaal  o^at  ploaom  aad  ooaalas) 

(c)  Is  aalactad  caaas  coafirm  facUlly  charaetarlatles  aad  capabUltlas  hf 
vlalt  or  other  contact. 

Moat  of  the  work  nadar  item  (a),  above  haa  bean  completed,  aad 
work  aadar  item  (b)  haa  bagua.  Pralimiaaiy  raaolta  cf  tea  aurvey  are  ahown 
ia  tea  matrix  of  flgara  3,  which  dallnoa.  in  qoalitiatlva  terma.  the  typea  of 
facility  boat  aoited  for  each  category  of  inatramant  teat.  Alao  ahown  are  the 
aamaa  of  facUitiaa  which  have  baaa  idaatiflad  aa  poaaaaaing  aoma  or  all  of  tea 
characteriatiea  raquirad.  All  of  the  aaada  of  a  particular  category  may  not  be 
filled  by  tea  indicated  facility:  if  bowaaar.  it  appaara  that  a  aignificant  fraction 
of  tea  designated  taating  can  be  aecompliahad.  tea  facility  has  bean  iocludad. 

For  example,  tea  ability  of  tea  ORNL  and  ZNEL  air  water  loopa  to  operate  at 
elevated  pressure  is  believed  limited;  nonatealesa.  these  loops  may  be  useful 
(aad  ia  tee  ease  of  tee  ORNL  ioop.  have  been  useful)  ia  development  of  inatru- 
naeat  correlatioas  aad  two«phasa  calibrations.  From  the  preliminary  information 
coacamiag  tee  facUitiaa  on  the  matrix  of  requiremaats.  t^  foUowing  tentative 
ceacluaiona  may  be  drawn: 

(1)  There  have  not  been  a  large  number  of  (single  phase)  gas  loops  identified. 

This  need  is  however  net  considered  crucial,  since  generally  the  gas/water 
loopa  of  the  4te  row  cait  be  operated  ia  a  ategle  phase  mode.  Also,  it  is 
believed  that  additional  faeUitlea  not  yet  identified  may  exist  (used  in  cali- 
bratioa  of  meters  for  gas  pipeliiMs). 

(2)  The  ability  to  operate  existing  small  gas/water  loops  at  high  pressure  appears 
limited:  teis  is  a  significant  gap  in  two>phase  calibration  capabiUties. 

(3)  No  large  sixe  gas/water  loops,  suitable  (or  scaling  studies  have  been 
identified.  This  may  be  a  significant  deficiency  for  programs  such  as 
LOFT. 

(4f  There  are  no  obvious  requirements  for  large  scale  steam/water  loops, 
teough  such  loops  might  in  some  instances  be  used  in  scaling  tests.  (If 
this  is  dons,  however,  instrumenting  wiU  be  difficult,  expenses  high,  and 
results  more  limited  than  a  comparaUe  test  in  air /water. ) 

Aa  argument  miitht  be  advaiwed  ia  favor  of  large  steam/ water 
loops  to  the  effect  that  they  constitute  final  *'proof>of-the-pudding’'  calibration 
of  tee  two'phase  instrument,  because  of  profUe  effects,  however,  such  caIi-> 
brations  are  of  littie  real  use,  uidess  upstream  and  downstream  hydraulic 
geometry  is  duplicated.  In  effect,  the  calibration  facility  wiU  come  to  resemble 
tee  safety  experiment  itself.  A  more  effective  final  test  of  the  instrument  is 
considered  to  be  accomplished  in  the  safety  experiment.  Here  the  instrxunent 
can  be  checked  against  others  of  similar  types  if  redundancy  in  instrumenting  is 
employed,  and  by  mass  balances  using  Tjuid  level  and  other  instruments. 
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(VOID  FRACTION,  PER  UNIT) 


